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ANALYTICAL AND EXPERIMENTAL 
INVESTIGATION OF A l/8-SCALE DYNAMIC 
MODEL OF THE SHUTTLE ORBITER 

Volume H - Technical Report 

By P. W. Mason, H. G. Harris, J. Zalesak, and M. Bernstein 

GRUMMAN AEROSPACE CORPORATION 
Bethpage, New York 11714 

1 - INTRODUCTION 
1.1 OVERVIEW AND SUMMARY 

The Space Shuttle configuration has more complex structural (Jynamic charac- 
teristics than previous launch vehicles, primarily because of the high modal density 
at low frequencies, and the high degree of directional coupling in the pitch plane. An 
accurate analytical representation of these characteristics is a primary means for 
treating structural dynamics problems during the design phase of the Shuttle pro- 
gram. The 1/8 -scale model program was developed to explore the adequacy of 
available analytical modeling technology, and to provide the means for investigat- 
ing problems which are more readily treated experimentally. The basic objectives 
of the l/8-scale model program are to 

• Provide early verification of analytical modeling procedures on a Shuttle- 
like structure 

• Demonstrate important vehicle dynamic characteristics of a typical Shuttle 
design 

• Disclose any previously unanticipated structural dynamic characteristics 

• Provide for development and demonstration of cost-effective prototype 
testing procedures. 

As listed in the FOREWORD, this four-volume report covers the analytical 
and experimental results relating to the Orbiter structure. The analytical studies 
pertaining to the External Tank (ET) and Solid Rocket Booster (SRB) are covered in 
References 9-1 and 9-2, respectively. 
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The mathematical modeling of the Orbiter began in January 1973. Consistent 
with the concepts of testing sub structuring capability in NASTRAN level 15.5 (using 
ALTERS where needed), the structure was divided into five components: fuselage, 
fin, cargo doors, payload, and wing. 

All skin material in the initial NASTRAN model (Model I) was considered to 
be nonbuckled, i. e. , fully effective. This approach is consistent with the design of 
the full-scale Shuttle Orbiter where any wrinkling of the skins would loosen the 
bonded thermal protection system (TPS). 

After a series of static and <tynamic tests of the experimental Orbiter model, 
it was concluded that the NASTRAN model should be modified to account for initial 
imperfections and poor joint designs. These changes in the finite element model 
(NASTRAN Model D) accounted for additional flexibility in six major areas 

• Fin/fuselage supports 

• Forward/mid-fuselage splice 

• Cargo door attachments 

• Wing carry-through structure 

• Effective width of fuselage and wing skins 

• Payload attachment. 

The effects of revising the NASTRAN model were substantial. Hie agreement 
between the second analysis and test is very good for the low frequencies, and as 
generally expected, discrepancies between theory and test become larger for the 
higher modes. Also, the static test data reported herein has a ±10% error in re- 
producibility. It is felt, however, that the agreement of Model II results with test 
data is sufficiently close such that a full coupling analysis may be pursued. 

In NASTRAN Model I the CQDMEM2 element was used to model a fully effec- 
tive structure. The major cause for lack of correlation was traced to the ineffec- 
tiveness of the skin panels due to the presence of initial bows. Model n used bars 
and shear panels to represent the behavior of the physical model. Here, effective 
bar areas and an effective shear modulus were used. It is felt that this type of 
modeling is more descriptive of the actual characteristics of the structure. The 
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effectiveness factors were obtained through the use of a computer program that 
solved the large deflection problem of in-plane loading of plates with initial imper- 
fections. Little information exists in published form that contain the parameter 
range and type of loading of interest to the present study. Work was done by NASA, 
some 20 years ago. It is recommended that this work be updated to cover a broad 
range of aspect ratios, initial imperfections and loadings. Publication of these data 
in chart form would be beneficial to modeling efforts on many structural projects. 

Some remodeling of major joints could have been made in the analysis. How- 
ever, for some joints such as the fin-to-fuselage connection, only static tests could 
give the exact behavior no matter how careful an analysis was made. 

The original analysis plan for the l/8-scale model required that all components 
(Orbiter, ET and two SRB’s) be coupled to determine mated vehicle modes. This 
combined hydro-visco-elastic analysis is theoretically possible within the present 
NASTRAN system, but from a practical stand point it is not. Although modes were 
successfully calculated for the individual components, it was impracticable to couple 
these components and analyze the total (mated) vehicle due to extremely large com- 
puter time requirements. Note that the lack of correlation of the Orbiter analytical 
and test modes for Model I is not a fault of NASTRAN. This lack of agreement did, 
however, cause a shift of emphasis from studying the coupling problem to examining 
the Orbiter correlation by the use of Model II. Coupling is still a major problem. Two 
approaches seem feasible: 

• Modal coupling in lieu of static coupling to reduce the size of the final 
problem 

• Incorporation into NASTRAN of approximate reduction schemes, such as 
FEER (Reference 9-3), that employ an automatic tridiagonal reduction 
algorithm. 

Master Agreement NAS 1-10635, Task 21 will pursue the modal coupling 
approach by making the appropriate ALTERS to the NASTRAN rigid formats. 

The analysis of the l/8-scale model was viewed as a pilot study of the use of 
NASTRAN on realistic aerospace projects. Therefore, some comments extend past 
specific application to the l/8-scale model. Many areas must be considered in 
attempting to judge the NASTRAN system, some of which are not necessarily con- 
cerned with NASTRAN itself. 
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One of the prime reasons for Grumman’s interest in NASTRAN is the work 
that NASA is doing in attempting to establish NASTRAN as the industry standard. 

For multicorporation aerospace projects, standardization is a necessity. Standard- 
ization of computer programs, type of data, etc. , is important from a contractual 
point of view, where subcontractors interface with prime contractors and the prime 
contractors interface with the principal agency. However, with standardization 
there must be sufficient technical flexibility - in this regard NASTRAN has room for 
improvement. Some overall comments concerning NASTRAN (Section 8 contains 
additional detail) are as follows: 

• The learning period required to become proficient with the NASTRAN 
system is lengthy 

• The five Orbiter substructures were coupled using NASTRAN level 15,5. 

No great difficulties were encountered in using NASTRAN to do this. How- 
ever, the column partitioning vectors required in the MERGE instruction 
are somewhat awkward. The equilibrium checks that were incorporated in 
the analysis proved helpful in finding numerical errors in the stiffness and 
mass matrixes and in giving a high level of confidence that there are no 
extraneous constraints. For the l/8-Scale Orbiter Model, five substruc- 
tures were combined to give the total structure. Multilevel sub structuring 
does not appear to present any technical difficulties 

• Existing eigenvalue routines in NASTRAN are inadequate to handle the large 
size problems that are associated with the coupled structure 

• Regardless of many objections, a Government/industry standard like 
NASTRAN is a necessity. 

1. 2 SYMBOLS AND ABBREVIATIONS 


Symbols 

5 deflection, in. 

^ij influence coefficient, i.e. , deflection of point i due to unit load 

at point j. , in. /lb 

f 

m 
1 -4 


frequency of model, Hz 


n 


L 

m 

L 

P 

t 

Abbreviations 

ALARM 

FEER 

ALTER 

CBAR 

CROD 

C SHEAR 

CTRMEM 

CQDMEM2 

CQUAD2 

DOF 

ET 

GLOW 



MPC 


NASTRAN 

SPC 

r. < V • 

SRB 

TPS 


frequency of prototype, Hz 


length of model 
length of prototype 


jnl 


SQS'H 


skin thickness of model, in. 


computer programs to solve eigenvalue problems of very 
large size 

modification of NASTRAN rigid formats 


finite elements in the NASTRAN program used in analysis 
of Orbiter structure 


degrees of freedom 
external Tank 
gross liftoff weight 
liquid Hydrogen 
liquid Oxygen 

multiple point constraint - NASTRAN 
NASA Structural Analysis System 
single point constraint - NASTRAN 
Solid Rocket Boosters 
thermal protection system 
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2 - DESCRIPTION OF PHYSICAL MODEL 

2.1 DESIGN OBJECTIVE 

The l/8-Scale Shuttle Dynamic Model is based on Grumman’ s proposed Design 
619 parallel-bum Space Shuttle, shown schematically in Fig. 2-1. Figure 2-2 shows 
a mockup of the l/8-Scale Shuttle Model basic configuration; Fig. 2-3 shows the 
detailed structural arrangement of the prototype Orbiter. In simplifying the 1/8- 
scale model design, a major objective was to keep model fabrication cost within tar- 
get, while retaining as many of the significant structural dynamic characteristics as 
possible. The alloted funds did not permit fabrication of an exact or ''replica” model 
at the small scale necessary for testing in the existing NASA/Langley facilities. 
Thus, only the general characteristics of the major Oibiter components have been 
simulated; there has been no attempt to model local details. Figure 2-4A is an 
external view of the l/8-Scale Orbiter Model. Figure 2-4B shows the model with 
nonstructural plastic fairings mounted. 

2.2 DESIGN CONSIDERATIONS 
2.2.1 SIMILITUDE REQUIREMENTS 

A first step in the modeling procedure was to determine the major factors that 
influence the response of the prototype being studied. For the Shuttle structure, this 
was accomplished under Task 4 of the Master Agreement Program (Reference 9-4). 

Table 2-1 summarizes pertinent scaling relationships that must exist between 
the model and the prototype. These follow directly from a dimensional analysis of 
the various parameters that influence the dynamic behavior of the structure, and 
from the choice of the model material. Extrapolating prototype behavior from 
model test data is accomplished by directly using these scaling relationships. Note, 
however, that because of design expediency some of the scaling rules have been 
compromised. For example, the local skin stiffness on the model is less than the 
required skin scaled value of the prototype to prevent buckling. Liberty has also 
been taken in modeling the stiffness characteristics by some necessary lumping to 
avoid the large expense of exact scaling of very small dimensions. Thus stiffeners, 
although lumped to some extent, have not been eliminated completely. If they were 
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Fig. 2-3 Prototype Orbiter (Grumman Design 619) Structural Arrangement 
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Fig. 2-4 Views of Finished 1/8-Scale Orbiter Model 
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completely eliminated, then the local stiffness of the skins would not be duplicated*, 
and premature buckling would occur. 

While accurate modeling of the prototype was desirable for extrapolating basic 
Shuttle dynamic characteristics, another prime study object was an evaluation of the 
NASTRAN dynamic analysis and its correlation with model test data. A complete 
static and dynamic analysis was made using NASTRAN with the structure modeled to 
a degree of refinement considered sufficient for preliminary design purposes. Thus, 
the need for direct scaling of the prototype design to obtain an exact model in every 
detail was not considered to be crucial. Note also that the Shuttle design was still in 
a state of flux at the beginning of this study, hence any attempt to model the then 
current vehicle exactly would not be overly beneficial to the Shuttle project. 


Table 2-1 Pertinent Scaling Relationships for the 1/8-Scale Model 


Physical Quantity 

Magnitude 

(a) 

Length and displacements 

L m 1 
L P = 8 


Poisson's ratio, v 

V = V 

m p 


Mass density 

=i 


Modulus of elasticity, E 




e m “ e p 


Area 

A m /IV 
A p = W 

1 

64 

Area moment of inertia, 1 


1 

4096 

Mass moment of inertia, 1' 
TT-1 


1 

32,768 


Physical Quantity 

Stress (a = Ee) 

Force (F = o A) 
Longitudinal stiffness, EA 

Bending stiffness, El 
Torsional stiffness, GJ 
Weight (W = pV) 


Acceleration (F = ma) 




£ 

1 


1 

64 

= 1/64 

= 1/4096 

= 1/4096 
1 

512 


Subscript "m" refers to the model; subscript "p" refers to prototype. 
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2.2.2 LOADS 
LOAD FACTORS 

Table 2-2 summarizes the principal load factors used in the design of the 
l/8-scale Orbiter Model. No attempt was made to limit the lateral and transverse 
loads during the vibration mode surveys. 

LOADS INDUCED BY HANDLING 

The handling load conditions were assumed critical in the design of the 1/8- 
Scale Orbiter Model because of the low overall stiffness of the structure. Table 2-2 
lists the load factors used in the design for the handling conditions. All raising and 
lowering of the model was accomplished with the model in a vertical orientation, 
supported from integral nose fittings. These stipulations were necessary to prevent 
permanent buckling in the lower skins. 


Table 2-2. - Summary of Load Factors Used in 1/8-Scale Model Design 


Condition 

Load Factors 

Limit, g 

Ultimate, g 

Impact during fabrication and assembly 

CM 

II 

X 

c 

3 

Cable and hoist fitting and attachments 

CM 

II 

X 

c 

3 

During raising and lowering model^ 

n x = 1.5 

2.25 

During mode surveys^ 

n x = 1.2 

1.8 

Iq\ 

Lateral forces during handling' ' 

n y = n z = 0.5 

0.75 


a Selected to account for dynamic effects during starting and stopping of the hoist 
system. 

^Selected to provide a 20% envelope on body bending and shears and local inter- 
face loads due to dynamic loading. 

c Value for lateral (Y) and transverse (Z) load factors was selected to account for 
handling forces. 

TT-2 

2.2.3 MASS AND STIFFNESS DISTRIBUTION 

Table 2-3 lists the Orbiter design weight distribution for one-half of the model 
structure. A schematic "stick figure" representation of the fuselage mass distribu- 
tion (for one-half of the structure) is shown in Fig. 2-5. Note that, as stated in 
Table 2-3, the fuselage weight of 100. 91 lb contains the aft OMS ballast of 26. 15 lb 
and the cabin ballast of 29.2 lb per one-half structure. The abort SRB (Fig. 2-1) is 
not included in the model. 
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Table 2-3 Design Weight Distribution for One-half of the 
1 /8-Scale Orbiter Model Structure 


Substructure 

Weight, lb 

Weight, % 

Fuselage 

100.91 ,a) 

48.34 

Wing 

32.86 

15.74 

Cargo doors 

6.52 

3.12 

Fin 

3.85 

1.84 

Payload 

64.62 

30.96 

Total 

208.76 (b > 

100.00 


a The fuselage weight of 100.91 lb contains the aft OMS 
ballast of 26.15 lb and the cabin ballast of 29.2 lb (per 
half structure) 

T"T*3 

°The actual measured weight was 420 lb for the Orbiter model 

Figure 2-6 compares the total axial cross-sectional areas of the l/8-Scale 
Orbiter Model with the scaled prototype Design 619. A similar comparison is made 
in Fig. 2-7 for the area moments of inertia about the Y-axis (vertical bending). 
Curve "A" represents the values used in NASTRAN Model I. Curve "B" represents 
the effect of removing the door longeron. Curve "C" represents NASTRAN Model 
II with partially effective skin and no door longeron. The drop in values from "B” 
to "C" is due to the presence of initial out-of-plane bows in the skin panels. The 
actual value of the inertia in Model n is higher due to the door longeron, but the 
change is not as great as the difference between curves "A" and "B". This is be- 
cause the door longeron is only partially effective, due in some part to the cargo 
door attachments. Curve "D" represents the prototype Design 619 values scaled by 
the ratio of 4096. 

Dead load shear, axial load, and moment diagrams (for the one-half Orbiter 
Structure) supported at the interstage are shown in Fig. 2-8; Fig. 2-9 presents 
axial dead-load diagrams for the half orbiter suspended from nose fittings. 
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"A" 1/8 SCALE MODEL - DOOR LONGERONS 
FULLY EFFECTIVE (NASTRAN MODEL I) 



Comparison of Fuselage Moments of Inertia of Prototype 
and 1 /8-Scale Orbiter Model 




Fig. 2-8 Dead Load Shear and Moment Diagrams for the 1 /8-Scale Orb iter Model 
Supported at the Interstage (Half Orb iter Values) 
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2.3 FABRICATION DETAILS 
2.3.1 FUSELAGE 

The external lines of the l/8-Scale model were simplified to reduce fabrica- 
tion cost, hence more flat side and bottom skins are used in the model than were 
employed in the design prototype. The model fuselage has been fabricated as three 
modules: a forward section, having a tapered non-circular stiffened shell; a U- 
shaped stiffened mid-fuselage section, with straight sides and bottom deck; and an 
aft section, with shell geometry similar to the prototype. (See Fig. 2-10.) 

The forward fuselage (Fig. 2-11) contains attachments for supporting ballast 
weights that simulate the mass of the forward equipment. The hoisting "pick ups' 1 
are attached to longerons which, in turn, are spliced to the mid-fuselage upper 
longeron at Station 64.0. Figure 2-12 shows the forward-to-mid-fuselage interface. 

The fuselage mid-section of the prototype consists of closely spaced frames 

covered by corrugation stiffened skin carrying TPS tiles. The model has similar 

(but simplified) geometry and structural arrangement, consisting of a series of U- 

shaped frames spaced at approximately 10 in. apart. The model skin thickness is 

0.020 in. on the side and 0.025 in. on the bottom. This thickness is enough to 
prevent buckling under static load (if the sheet is perfectly flat), but it results in a 

larger overall section area and moment of inertia than the scaled prototype values 
(as shown in Fig. 2-6 and 2-7). Payload support attachments are provided to permit 
four variations in the payload configuration. The longerons are designed to furnish 
the proper scaled area for the aft end of the fuselage, and are kept constant in area 
for the entire length of the fuselage to limit fabrication costs. To simulate the proto- 
type weight, a distribution of about 0. 8 lb/in. , including the structural weight, was 
required; it was obtained by increasing the thickness of the frame webs where stiff- 
ness is not a factor. 

The aft fuselage, Fig. 2-13 and 2-14, is geometrically similar to the proto- 
type, but contains only two full bulkheads, one aft frame, and two intermediate 
frames. Neither the abort SRB's, nor the cutouts for them and for the OMS in the 
prototype, were modeled. The cross-sectional areas of the struts between the lower 
engines and the aft interstage fittings (Fig. 2-13) are scaled directly from the proto- 
type areas. 
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retail* of Forward 
Fuselage Module 


Flg 2-12 Forward-To-Mid- 

Fuselage Interface 
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T11 Fig. 2-13 Aft Fuselage Details 


The fuselage frames in the cargo bay region are constructed of aluminum 
sheet bent to form a channel section. The tapered side wall channel section and the 
lower portion are attached back-to-back to form a U-shaped frame (Fig. 2-10). The 
material used to form the frames is sufficiently thick to ease fabrication and to avoid 
buckling. In addition, the frames are heavy enough to simulate some of the non- 
structural weight in the prototype. The details of the prototype bulkheads were not 
scaled to avoid excessive costs, and because these details are not significant in 
establishing the primary structural dynamic characteristics of the model. Figure 
2-14, the aft fuselage during assembly, shows the three simulated engine mounts. 

2.3.2 WINGS 

The delta wings, shown under construction in Fig. 2-15 and 2-16, consist of 
six spars and four ribs formed from .032 -in. 2024 aluminum sheet. The covers are 
0.020 in. thick, constant over both upper and lower surfaces. The wing tapers in 
depth from 2.5 in. at the tip (BL 61.58) to 6.0 in. at the root (BL 13.75), with the top 
cover at a constant elevation (WL 51.5). Shear fittings connect the closure rib at 
BL 13.75 to the fuselage frames (the side wall of the fuselage is at BL 12.5). The 
bottom cover is spliced to the fuselage bottom skins and the upper cover is spliced 
to the fuselage carry-through deck extending from station 144. 75 to Station 170. 75. 



T ' 12 Fig. 2-14 Aft Fuselage — Engine Mounts 
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Fig. 2-16 Fuselage Side Wall and Wing Shear Attachments 
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This means that full bending continuity exists only between the last four spars in the 
wing. To some extent, this simulates the behavior of the prototype where the wheel 
well in the forward portion of the wing requires a large cut-out in the lower surface. 
Thus, although the wing is a delta planform, it tends to act more like a high aspect 
ratio structure. The proper weight distribution on the wing has been obtained by 
adding steel ballast plates to the rib and spar webs. 

2.3.3 VERTICAL FIN 

The vertical fin, Fig. 2-17 and 2-18, contains three spars and a closure rib 
that supports the fin ballast weight. The ballast weight is located at the scaled proto- 
type center of gravity, hence the model fin is much smaller in planform than would 
be required if the total fin were modeled. The spar webs are 0.032 in. thick, and 
the covers are 0.020 in. thick. The forward spar is connected to the fuselage bulk- 
head at Station 166.5 by a 0.040 -in. thick angle fitting. The aft spar is connected by 
a similar angle fitting to the fuselage closure frame which is relatively flexible, com- 
pared to the two bulkheads that support the forward and mid-spars. The central 
spar is connected to the engine support bulkhead through gusset panels from the spar 
flanges to clips attached to an angled doubler. 

2.3.4 CARGO DOORS 

The cargo doors (Fig. 2-19) consist of a removable seven-segment semicircular 
cylindrical shell. The semicircular frames at the ends of each segment have a com- 
mon lower flange and independent upper flanges attached to adjacent shell segments. 
These V-shaped frames allow the door segments to breathe in the fore-aft direction, 
yet act in resisting shear forces. The door segments and frames are connected to 
the door longeron which, in turn, is attached to the fuselage upper longeron by four 
keystone shaped fittings (see Fig. 2-20). The door shell segments are 0.016-in. 
thick, compared to a scaled prototype value of 0.00325 in. The heavier gage is re- 
quired to prevent buckling. 

2.3.5 PAYLOAD, BALLAST WEIGHTS, ATTACHMENTS AND OTHER DETAILS 

The payload (Fig. 2-21 and 2-22) consists of two channels intermittently welded to- 
gether at the flanges to form a closed box. Fittings were then attached to this box to sup- 
port it in a statically determinate fashion within the cargo bay area of the fuselage. The 
model is constructed to accommodate four types of payload; however, only the full-up 
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payload case (equivalent to 65,000 lb in the prototype) was analyzed. For this con- 
dition, the forward payload support (FS 78.0) provides vertical and side load re- 
actions. At the aft end, outrigger brackets provide two vertical supports connecting 
the payload to the fuselage sidewall at FS 151.875, BL ± 10.125. Drag and side load 
are taken out at the aft end of the payload on the center line. Typical payload attach- 
ments are shown in Fig. 2-22. 

Major ballast weights (Fig. 2-21) consist of the cabin ballast, of 58.4 lb; two 
OMS ballast weights, 29.2 lb each; and the fin ballast, 6.0 lb. 

A difficult task in designing a scaled model is the simulation of the behavior of 
fittings and connections to the proper scale. In the 1 /8-Scale Orbiter Model, joints 
and attachments have been simplified primarily to reduce fabrication costs. Basi- 
cally, it is believed that the model joints are stiffer than the prototype, except in 
certain locations such as the fin-fuselage connection, the forward-mid-fuselage 
splice, and the forward payload attachment. Analytical modeling of the cargo door 
attachments proved to be rather subtle. The torsional stiffness of the fuselage de- 
pends greatly on the behavior of these attachments. 
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PAYLOAD 


LOCATION OF 
OMS WEIGHT 


Fig. 2-21 Location of Main Ballast Weights 


PAYLOAD] 


NOTE 


This configuration not used in tests to date. Presented to show typical details only, 
Fig. 2-22 Payload Typical Attachment Details 
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3 - FINITE ELEMENT ANALYSIS PROCEDURES 


3.1 BASIC PHILOSOPHY 

The following guide lines were established in setting up the finite element model 
and analysis procedures: 

(1) The entire vehicle would be analyzed using NASTRAN. 

(2) The model would be of sufficient refinement to adequately 
predict overall dynamic behavior. No attempt would be 
made to try to predict local panel motions. 

(3) The modeling detail would be of sufficient refinement to allow 
for the prediction of internal load distributions that would be 
adequate for preliminary design of the structure. Although there 
was no intention of computing internal loads, the analysis was 
considered to be representative of an actual prototype design 
situation. Interest, therefore, existed in evaluating how NASTRAN 
would blend into a design environment. 

(4) The total structure would be analyzed by employing sub structuring 
techniques to again evaluate how well NASTRAN could handle these 
procedures in a design environment. NASTRAN could, in principle, 
handle the entire Shuttle as a single unit; realistically it can not. 

The total vehicle contains the Orbiter, External Tank (ET) and two 
Solid Rocket Boosters (SRB's). This combined hydro-visco-elastic 
problem required far too much computer time to be considered as 

a single unit. The Orbiter by itself could have been considered 

as a single unit, but again, this did not represent a realistic situation. 

(5) Analyses of the ET and SRB's would be performed under separate 

> q, 

contract to investigate the hydro-elastic and visco-elastic capabilities 
of NASTRAN (References 9-1 and 9-2). All four bodies (Orbiter, ET, 
and SRB's) would be coupled using static substructuring techniques to 
obtain an analytical model of the launch configuration. 
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(6) The NASTRAN weight analysis capability would be used to calculate 
the individual component and total weights for the nonfluid portions 
of the model. A supplementary weight check would be conducted 
and the NASTRAN results adjusted where necessary. Structural 
grid points would be used as dynamic mass points using the Guyan 
reduction method as required. Note that this procedure differs from 
Grumman's usual practice, which is to establish a weights model 
independent of the structural model. In this approach, unit loads on the 
weights model mass points are then beamed to appropriate structural 
node points. The dynamic model is thus the same as the weights model 
or a subset of it. This procedure inherently results in a smaller dynamic 
model and additional reduction schemes are not necessary. The 
equivalent reduction takes place in the beaming of unit loads from the 

t 

weights model to the structural model. This method was not used 
because 

- More extensive alters to the NASTRAN rigid formats would have been 
required 

- The NASTRAN weight analysis capability would not have been used 

- Basic mode data would have been produced at nonphysical points 
which might hinder test correlation. 

3.2 OVERALL ANALYSIS FLOW 

The overall analysis flow, Fig. 3-1, represents the originally proposed analysis 
for the combined total vehicle. The Orbiter was divided into five substructures: 
fuselage, cargo doors, fin, wing and payload. The ET was divided into two substruc- 
tures: the LOX tank and the aft portion of the external tank (consisting of the inter- 
tank skirt, LH 2 tank, and aft tank skirt). The SRB was to be handled as a single unit 
(consisting of the forward skirt, propellant cylinder and propellant, and the aft skirt). 

Several technical problems arose during the study which prevented the comple- 
tion of the proposed overall analysis, namely: 

• The Orbiter analysis was completed at the same time that initial 
test results were made available. A rather poor correlation was 
shown to exist for the Orbiter alone 
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• Separate analyses were performed for the following SRB and ET 

structures: 

- LQX tank alone 

- Total ET, comprising the LOX and LH 2 tanks coupled with the 
intertank and aft tank skirts 

- SRB alone. 

The SRB had to be divided into two substructures, and then coupled using an 
"ALTERed" Rigid Format 7. Storage problems were encountered in attempting to ana- 
lyze the total SRB. Analysis of the total tank structure using the ALTERed Rigid 
Format 7 uncovered the practical inadequacy of NASTRAN in solving large hydroelas- 
tic problems. One mode was obtained after using what was considered to be an ex- 
tremely large amount of computer time. 

The two major problems encountered (lack of correlation of analysis and test 
data for the orbiter; excessive computer time requirements for coupling the total ve- 
hicle) forced a decision to abandon the original overall analysis flow. Consequently, 
basic effort was redirected to rectifying the Orbiter analysis to obtain correlation with 
test results; also a separate study (Contract NAS 1-10685-21 with Grumman) was in- 
itiated to find a means for improving the efficiency of the hydro-elastic analysis. Up- 
. 

on completion of this parallel study, current plans are to couple the components of 
the total vehicle using modal coupling procedures, rather than the standard type static 
procedures indicated in Fig. 3-1. 

J 

In the context of this report, the original Orbiter analytical model is designated 
"Model I"; the modified analytical model, which reflects test behavior, is referred to 
as ''Model II. " The analytical procedures used for both models were the same, ex- 
cept where noted. 

Referring to the Orbiter portion of Fig. 3-1, observe that each of the five Or- 
biter substructures was analyzed to produce reduced mass and stiffness matrices for 
selected dynamic degrees of freedom (DOF's) and interface attachment points. Modes 

for these substructures were then obtained with the interfaces held - an exception is 
the fuselage, which was analyzed in a free-free condition. This approach aided in 

checking and understanding the behavior of the combined Orbiter vehicle. Next, the 
five substructure stiffness and mass matrices were merged to form the total Orbiter 
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mass and stiffness matrices. These matrices were again reduced to yield final stiff- 
ness and mass matrices that were used in the modal analysis. 

3.3 SUBSTRUCTURING PROCEDURE 

The basic sub structuring procedure for combining substructures as presented 
in the NASTRAN User's Manual were followed, but with some minor changes in the 
assumptions used, and with more extensive DMAP ALTERS. These ALTERS were 
written for Rigid Format 3 to provide for a more reliable substructuring procedure. 

The analysis was performed in two phases, as shown in Fig. 3-2. In NASTRAN 
Phase 1, each substructure was analyzed and checked separately. The output from 
this phase was then assembled onto a copy tape for the symmetric and antisymmetric 
cases, and coupled in NASTRAN Phase 2. The same set of ALTERS are used for 
both phases, the particular phase being determined by an input parameter. 

The following assumptions were used in formulating this substructuring pro- 
cedure: 

• Any external supports present are included in the Analysis Set (a-set) 

• Any zero stiffness DOF's and symmetric or antisymmetric boundary 
constraints at the model plane of symmetry are included in the Single 
Point Constraint Set (SPC-set). No other DOF's are included in this 
set 

• Masses which are associated with zero stiffness DOF's will be lost 
unless these DOF's are "beamed" to adjacent points using Multipoint 
Constraints (MPC's) 

• The interface DOF's may be sequenced differently and in different 
coordinate systems in any two substructures to be coupled. The MPC's 
are used to relate the appropriate DOF's irrespective of local coordinate 
systems or initial sequencing 

• Grid points, rather than scalar points, were used in Phase 2. This 
allowed total Orbiter mode plots to be obtained in this phase, rather 
than obtaining substructure mode plots in a NASTRAN Phase 3 run. 

PLOTEL elements were used to define the Phase 2 Orbiter pseudo 
structure 
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• The grid cards used in Phase 2 were the same as the grid cards used in 
Phase 1 that contained any a-set DOF. Any DOF's not in the Phase 1 a-set 
were removed by SPC cards. Each substructure had a unique set of grid 
identification numbers, and referenced a common basic coordinate system. 

Although the general theory presented in the NASTRAN User’s Manual for sub- 
structuring is correct, it does not provide analysis checks at various critical points 
in the procedure. Consequently, the following checks were incorporated in the Model 
I analysis by means of DMAP ALTERS: 

• A rigid body check was made in Phase 1 after the generation of the 
reduced stiffness and mass matrices. Temporary rigid body supports were 
included in the deck as "support" cards for this purpose 

• The structural transformation matrixes G_, G , and D were used to gen- 

m’ o 

erate equilibrium matrices for the various constraint sets, except SPC's. 
These equilibrium matrices represent resultants about a chosen origin due 
to unit applied loads at the appropriate DOF's 

• Provision was made to compute either free-free modes, or modes with 
the substructure held at the interface. This is necessary if each sub- 
structure is to be checked independently in Phase 1 

• A rigid body mass matrix relative to the basic origin was computed. This 
was compared with the general mass matrix calculated by the Grid Point 
Weight Generator (GPWG). This check verified that no mass had been 
lost in the reduction process. 


The DMAP statements to perform these functions for Rigid Format 3 were pre- 
sented in Volume IIEB, Appendix Bl, of this report. 

In the Model n analysis, a revised ALTERED Rigid Format 3 was used. The 
new alters provide more convenient checks including an SPC check. These ALTERS 
require a revision in the GPWG module (refer to Volume MA, Appendix A4) to obtain 
a transformation matrix that is generated from geometric data. This matrix ex- 
presses the resultant loads in terms of unit G-set loads about a chosen origin that is 
defined by the GRDPNT parameter. 
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Fig. 3-2 Flow Diagram for NASTRAN Substructuring to Obtain Normal Modes 
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The rigid body error ratio was removed since the support card could be used to 
define temporary supports to obtain restrained modes. This eliminated the necessity 
of inputing a column partitioning vector to remove the temporary supports from the 
a-set DOF's. It was felt that the reduced stiffness check and the incorporated rigid 
body transformation checks were more helpful in pin-pointing errors. The rigid body 
transformation matrices G m , G q , and D were related to the rigid body motions of the 
origin and similar independently derived matrices were subtracted to give rigid body 
transformation checks. 

A rigid body mass matrix was generated from the reduced mass matrix M and 
compared with the M q matrix that is output from the GPWG. This approach ensured 
that no weight had been lost in the DOF reduction process. 

The revised ALTERS to Rigid Format 3 are given in Volume mA, Appendix A4. 
3.4 FINITE ELEMENT MODEL I 

Table 3-1 shows the number of grid points and types of finite elements used in 
Model I for the five Orbiter substructures. Sketches of the finite element model are 
presented in Appendix A of this volume. These drawings were prepared to establish 
the finite element model geometry and member data. In the discussion that follows, 
reference to these sketches is implied. 

NOTE 

Appendix A includes a list of the sketches it contains. The 
summary list included here augments the general drawing 
titles in terms of discussion specifics. Readers should con- 
sult these "background" sketches as they may desire, even 
though the text does not, in every case, cite specific Appendix 
A illustrations. 

• Orbiter general arrangement, Fig. A-l 

• Fuselage frames (total of 21), Fig. A -2 through A-4 

• Fuselage sidewall (skin), bottom deck, keel, and wing carry-through; 
also, cargo door frames, and shell segments, Fig. A-5 

• Orbiter wing. Fig. A -6 

• Fin and payload idealizations. Fig. A-7. 
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The fuselage shell structure was modeled using CQDMEM2 elements - a new 
element in NASTRAN, but one that has been used widely at Grumman. As Fig. 3-3 
shows, it is a warped quadrilateral composed of four constant strain triangles having 
a common central node defined by the intersection of lines that connect the midpoints 
of the opposite sides of the quadrilateral. Figure 3-4 illustrates the suitability of 
this element to represent the fuselage side wall. The curves in this figure compare 
the effective stiffness of a composite beam segment subjected to pure bending with the 
conventional beam formula, El = In these calculations the node displacements 

along the sides have been constrained to remain linear by MPC's. With no caps, the 
error in stiffness becomes larger with increasing aspect ratio. When caps are added, 
they tend to dominate the bending behavior and reduce the error. The dotted region 
represents the locus of beam segments for the fuselage side wall. (A single 
CQDMEM2 element of aspect ratio 1 is 70% too stiff in pure bending, and hence 
should not be used to idealize the web of a spar or rib.) Note that Fig. 3-4 implies 
that the expected error for the fuselage side wall should be in the neighborhood of 
10% or less. But the curve in Fig. 3-4 is only indicative of the behavior of a com- 
posite beam segment. Calculations performed on the fuselage section at Station 116 
indicate that the finite element model is only 1.2% stiffer than the actual El value. 
(Note also that the finite element model contains more segments than used in develop- 
ing Fig. 3-4, and the neutral axis is not at the midheight of the section. ) 


4 CONSTANT 

STRAIN 

TRIANGLES 



1. PARALLELOGRAM FORMED 
BY CONNECTING MIDPOINTS OF 
ADJACENT SIDES. DEFINES THE 
REFERENCE PLANE FOR THE 
WARPED QUADRILATERAL 

2. LOCAL X AXIS IS THE 
PROJECTION OF LINE 1-2 ON THE 
REFERENCE PLANE. THE LOCAL Y 
AXIS IS PERPENDICULAR TO X 
AND LIES IN THE REFERENCE 
PLANE. THE LOCAL Z AXIS IS 
PERPENDICULAR TO THE 
REFERENCE PLANE. 

3. POINT 5 IS THE 
INTERSECTION OF THE LINES 
CONNECTING THE MIDPOINTS OF 
OPPOSITE SIDES OF THE WARPED 
QUADRILATERAL AND LIES IN 
THE REFERENCE PLANE. 

4. THE LOCAL X AND Y 
CONSTRAINTS AT 5 ARE RELAXED 
ELASTICALLY 

5. THE LOCAL Z DOF AT 5 IS 
SLAVED TO THE LOCAL Z DOF'S AT 
POINTS 1,2, 3 AND 4. FOR 
EXAMPLE, A UNIT Z MOTION AT 1 
PRODUCES A Z MOTION AT 5 
EQUAL TO 1/4. UNIT Z MOTIONS AT 
2, 3 AND 4 LIKEWISE PRODUCE A Z 
MOTION AT 5 EQUAL TO 1/4. 


Fig. 3-3 Basic Assumptions Used in Developing the Warped Quadrilateral "CQDMEM 2" Finite Element 
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Figure 3-5 is an alternate rod/ shear panel idealization. Here, the side wall 
has been modeled with five rods and four shear panels, wherein the areas of the rods 
represent the "lumped" adjacent material in the web. A simple calculation indicates 
that this model overestimates the bending inertia by 12.5%. This error is greater 



1.0 2.0 3.0 4.0 5.0 

ASPECT RATIO, L/a 

T-25 Fig. 3-4 Behavior of a Beam Segment Composed of Four “CODMEN 2 " Elements and Two Rods 



Fig. 3-5 Alternate Rod and Shear Panel Idealization For Fuselage Side Wall 
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than the error associated with 4 CQDMEM2 elements of aspect ratio 1 with no cap 
material. One might be tempted to adjust the areas of the rod and panel idealization 
to obtain the correct inertia; however, care must be taken in doing this or else the 
axial stiffness and lateral inertia will be incorrect. With these guide lines in mind, 
the CQDMEM2 elements were used for modeling the sidewalls. Note that in the fore- 
going discussion, there is a tacit assumption that the material is unbuckled and fully 
effective. As pointed out earlier, this assumption follows from the criteria used for 
the prototype design. 


The U-shaped fuselage frames and the keel were idealized using CROD and 
CSHEAR elements. Here, effective cap areas were calculated for the CROD elements 
to represent the appropriate bending behavior. 


In idealizing bending members with rods and panels, some thought must be 
given to the refinement of the model. As an example, consider the simple cantilever 
beam illustrated in Fig. 3-6 where one-sixth of the cross-sectional area (th) of the 
bending material has been "lumped" into each of two equivalent caps (CROD elements) 
to preserve the inertia, and the shear material has been idealized as a shear panel. 
The CROD element uses the average shear load in evaluating the stiffness rather than 
a linear varying distribution. Figure 3-7 shows the error associated with this ap- 
proach. This figure compares the strain energy of a constant stress bar with the 
strain energy for a linearly varying stressed bar (an element usually associated with 
a force method bar and panel idealization). For this simple case the ratio, R, is 
zero; from the curve the error is 25 %. This means that the bending deflection is 
actually 33 % greater ( * ^ where € is the percent error from the curve). The only 
means for improvement, without using a different element, is to further refine the 
model as shown in Fig. 3-6. The ratio, R, for each of the four top cap rods is 0, 

1/2, 2/3 and 3/4 for rod numbers 1, 2, 3 and 4, respectively; the errors are 25, 3.6, 
1.3 and 0. 68 %, respectively. The combined error is only 1.57 % because the fourth 
element contributes the most to the strain energy. Note that the shear deformation 
is the same for either one large panel or the summation of four smaller panels. Thus, 
in terms of idealizing beam bending behavior, the limitation placed on the aspect ratio 
of a panel is predominately in terms of a control on the number of rods, and not so 
much on the shape of a panel. 
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As indicated previously, CQDMEM2 elements should not be used to idealize 
beam webs - they are too stiff in undergoing this type of behavior. Since CQDMEM 
elements are even stiffer, the same comment applies. It is advisable not to use a 
CQDMEM1 (isoparametric element) to represent the web of a beam, even if the direct 
stress properties are set to zero in attempting to simulate a shear panel. The iso- 
parametric element under these conditions does not behave in the same fashion as the 
conventional Garvey shear panel (element CSHEAR). 

CBAR elements were used to represent thin-ring type of frames, such as the 
engine compartment closure frame (Appendix A, Fig. A -4). These elements, 
with offsets, were also used to represent the fuselage upper longeron. 

The rib and spar webs of the wing and fin (fabricated from bent sheet) were 
idealized with CSHEAR elements. Again the effective rib and spar bending material 
was incorporated into CROD elements located in the covers. The covers themselves 
were represented by CQDMEM2 elements and, in the case of the wing, some CTREM 
elements located at the leading edge. Intermediate node lines were established for 
the wing and fin to further refine the grid. For the wing, the geometry of these node 
lines was established by the location of fuselage frames. 

The shell portion of the cargo bay doors was idealized using CQDMEM2 ele- 
ments, except for a few CQUAD2 elements that were required for local stability at 
Stations 102.12 and 129.0. The door frames were idealized with CSHEAR and CROD 
elements. Note that these frames contain two webs (Fig. 3-8), one common lower 
cap, and two upper caps that connect to the forward and aft shell segments. This al- 
lows the doors to "breathe" in the longitudinal direction. The door frames and shell 
segments were connected to the door longeron represented by CBAR elements. The 
longeron was, in turn, connected to the fuselage at four discrete points with MPCs. 

In Model n, the details of this region were changed: the door longeron was set at 
WL 62. 81, while the fuselage longeron remained at WL 62.5. 

The payload was idealized as a series of CBAR elements extending from Sta- 
tion 76. 0 to 151. 875. The fore and aft fittings that support the payload were assumed 
to be rigid. (Data for the MPC's is given in Appendix A, Fig. A-7.) 

NASTRAN plots of the finite element model are shown in Fig. 3-9 through 3-14. 
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Fig. 3-8 Schematic of Cargo Door Frames 



T -30 Fig. 3-9 NASTRAN Plot of Orbiter Fuselage 


3-15 













n 



3.5 MASS DATA 

A mass matrix was generated by NASTRAN for each of the substructures 
using weight data in three forms 

• Assigned as structural weight density on MAT1 cards 

• Assigned as nonstructural weight on member property cards 

• Lumped into concentrated weights. 

The sketches in Fig. 3-15 illustrate how the weight was distributed to the various 
structural grid points. Table 3-2 gives the half-structure Orbiter weights and centers 
of gravity of the various substructures. Table 3-1 indicates the number of DOF's for 
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FLANGE-ROD (p^' 0.0016 LB/IN. 3 ) 
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COVERS 



COVERS (p g - 0.1 LB/IN. 3 ) 





PAYLOAD BEAM 


BALLAST 


WEB FLANGE 


WEBS (p g “ 0.1 LB/IN. 3 ) 


(p N - 0.0135 LB/IN.* FOR BALLAST) 


•NOTE: NO BALLAST ON WEBS AT 
BL 29.7 AND BL 46.5 


WEIGHT 1 WING - 32.3 LB 

( STA 151.04 
BL 29.80 
WL 49.08 


FRAME WT. ASSIGNED TO 
OUTSIDE RODS (p N - 0.0088 LB/IN. 3 ) 


CARGO DOOR 


COVER 


FRAMES 


PAY LO AD-C8 AR (p s - 0.1 LB/IN. 3 ) 
(p., - 0.245 LB/IN. FOR BALLAST) 


3NCENTRATED WEIGHTS (BRACKETS & OVERHANG) 1/2 STRUCTURE 


* NOTE: THE TORSIONAL WEIGHT 

l ' INERTIA OF PAYLOAD BEAM 8i 

BALLAST WAS INCLUDED AS 
!4 CONCENTRATED WEIGHT AT GRID 

POINTS ALONG BEAM. BASED ON 
LB-IN. 1 

4.3 

IN. 

1/2 PAYLOAD WEIGHT - 64.6 LB 
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WL 62.23 
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.GRID 4400 l x (ROLL) = 1.56 LB-IN 3 
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COVER (p s « 0.1 LB/IN. 3 ) 
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SKIN LINE AS NON-STRUCTURAL WEIGHT 
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OF MEMBER DATA IN VOL III 
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301 .39 ADDITIONAL FWD INTERSTAGE ATTACHEMENT MATERIAL 
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2.25 AFT C INTERSTAGE FITTING 
.13 

.33 ADDITIONAL ENGINE THRUST BACKUP MATERIAL 
.22 
.12 
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• P N - ASSIGNED AS NON-STRUCTURAL WEIGHT ON MEMBER PROPERTY CARD. 

C. FUSELAGE 


Fig. 3-15 Model Weight Data 
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Table 3-2 NASTRAN Model Weights and Center -of-Gravity Locations 


Substructure 

Weight, lb 
(a) 

Center of Gravity 

X 

Y 

Z 

Fuselage 


116.36 

8.43 

58.37 

Wing 


151.04 

29.80 

49.08 

Cargo door 

6.4 

115.30 

8.22 

69.78 

Payload 

64.6 

116.02 

- 

62.23 

Fin 

3.8 

182.73 

- 

86.54 

Total (1/2 Orbiter) 

207.9 

122.94 

- 

59.05 


TT-6 Substructure weights are for half structure (including ballast) 


the various substructures after SPC's and MPC's have been applied, and after the 
Guyan reduction. The final number of DOF's for the symmetric and antisymmetric 
cases were 339 and 324, respectively. Figure 3-16 illustrates the location of the 
retained DOF's that were used in the dynamic analysis (Model I). 

3.6 COMPONENT RESULTS FOR MODEL I 

The results for the individual substructures were obtained from the Phase 1 
analyses. This output was generated to facilitate understanding the behavior of the 
substructures and the total vehicle. Subsection 5.5 discusses the use of the fuselage 
modes in conjunction with a simple beam model to aide in rationalizing the behavior 
of the total vehicle. Modes were calculated for the substructures either restrained 
at their interconnection, or in a free-free state, whichever seemed more appropriate. 
The substructure modes also helped in uncovering input errors, such as an incorrect 
constraint specification. An additional output from Phase 1 consisted of the weight 
and center-of-gravity location. This helped to ensure that all mass data was included 
prior to analyzing the total vehicle. In this regard, it should be noted that although 
an SPC can properly constrain a nonstrainable DOF in the stiffness matrix (such as 
motion normal to a bulkhead), any mass located at this DOF will be lost. Mass 
should not be placed at these DOF's and MPC's should be used to constrain the mo- 
tion. Obtaining the substructure weight and center-of-gravity data from the sub- 
structure mass matrix, using rigid body transformations, thus provides a valuable 
check. 

The lowest frequencies for the various substructures are indicated in Table 3-3 
with the corresponding support and boundary conditions. Of particular interest are 
the values for the fuselage. These were checked by alternate means (described in 
Section 5) to ensure their validity. 
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Orbiter Phase 2 Degrees-of-Freedom 
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Fig. 3-16 Orbiter Phase 2 Degrees-of-Freedom - Model I (Sheet 2 of 2) 
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Table 3-3 Lowest Frequencies for Orbiter Substructures — Model I 



Frequency, Hz 

Substructures 

1st 

2nd 

3rd 

• Fuselage 




— Free, Symmetric 

62.2 

129.9 


— Free, Antisymmetric 

28.7 

89.1 

128.3 

• Payload 




— Restrained, Symmetric 

81.2 

268.5 

627.7 

— Restrained, Antisymmetric 

68.6 

175.4 

462.8 

• Cargo doors 




— Free, Symmetric 

4.6 

10.7 

17.6 

— Restrained, Antisymmetric 

156.4 

622.2 

1054.6 

• Wing 




— Restrained 

77.6 

158 3 

259.9 

• Fin 




— Restrained, Symmetric 

264.2 

841.3 

1263.3 

— Restrained, Antisymmetric 

107.8 

407.2 

1018.7 


TT-6 


3.7 ORBITER RESULTS FOR MODEL I 

In Phase n of the coupling analysis, the five Orbiter substructures were coupled 
and modes were calculated. PLOTEL elements were used to connect the reduced dy- 
namic grid points for plotting purposes. Figures 3-17 and 3-18 are NASTRAN plots 
of the first two symmetric modes. These plots show the deformed structure together 
with the X, Y and Z displacement components measured from the undeformed posi- 
tion. The first mode (at 53. 2 Hz) exhibits fuselage vertical bending, fin pitching, 
and wing motion. Most of the wing motion is due to bending in the root region; the 
outer portion remained almost linear. Note that the point of maximum deflection is 
at the fin tip. 

The second mode (at 62. 6 Hz) consists principally of wing bending together with 
some payload and fuselage vertioal bending. (Additional plots of higher modes are 
shown in Volume IIIB, Appendix B19. ) Frequencies for the symmetric and anti- 
symmetric cases are given in Table 3-4. 
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A. f , = 53.2 Hz Fuselage Side Wall. Wing and Fin 



B. f , = 53.2 Hz Fuselage Bottom Deck, Wing Lower Cover, and Payload 

Fig. 3-17 First Symmetric Mode 
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A. f a - 62.6 Hz Fuselage Side Wall, Wing and Fin 



B. f j = 62.6 Fuselage Bottom Deck, Wing Lower Cover, and Payload 


Fig. 3-18 Second Symmetric Mode 



Table 3-4 Frequencies for Symmetric and Antisymmetric 
Modes — Model I 


Mode 

Frequency, Hz 
(Symmetric Case) 

Frequency, Hz 
(Antisymmetric Case) 

1 

53.2 

52.9 

2 

62.6 

72.6 

3 

75.2 

85.1 

4 

108.5 

92.0 

5 

133.7 

101.5 

6 

156.3 

135.3 

7 

162.2 

- 

8 

175.2 

- 

9 

216.4 

- 
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The coupled Phase 2 stiffness matrix formed for dynamic analysis was also 
used to calculate static deflections for comparison with data measured during the 
horizontal static test. This was accomplished by inverting the stiffness matrix re- 
strained at the test support points and calculating the deflections for loads applied at 
the test locations. The DMAP ALTER statements to do this are presented in Vol- 
ume in, pages A21-38 and A21-39. 
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4 - EXPERIMENTAL INVESTIGATION 

4. 1 TESTING PROCEDURES 

All static and dynamic testing has been conducted at the NASA/Langley Research 
Center. Three series of tests have been performed to date on the l/8-Scale Orbiter 
Model, and this series forms the basis for the correlation discussed in this report. 
Additional tests of the l/8-Scale model are planned on a time-shared basis with on- 
going tests of other components (ET and SRB). 

The main objective of the first series of tests was to determine the symmetric 
free-free modes and frequencies. For these tests, the model was suspended horizon- 
tally and inverted on soft springs attached to the interstage fittings (Fig. 4-1 and 4-2. ) 

In the second series of tests, the model was suspended vertically from two inte- 
gral nose fittings. A. "wiffle tree" arrangement was used to apply an axial preload 
to the model through the lower engine thrust structure and the fuselage bottom skin 
(Fig. 4-3 through 4-5.) Strain measurements were made to determine the behavior 
of selected typical panels under the effect of the axial preload (Fig. 4-6). The second 
series of mode surveys was conducted to determine the effect of the model dead weight 
and preload on the natural frequencies. 

The last of the three series of tests was performed with the model in a horizon- 
tal position, right-side up, and supported at the interstage fittings. The model was 
loaded iocrementally with statically applied loads at a number of locations. Deflec- 
tion measurements were made over the whole model (Fig. 4-7) for each applied load 
during loading and unloading. This data was used to obtain influence coefficients for 
comparison with the analytical results. The analytical results are presented in Appen- 
dix B of this volume. 

Details on the testing procedures are given in Reference 9-5. 

4.2 MODE SURVEY WITH ORBITER MODEL SUSPENDED HORIZONTALLY AT IN- 
TERSTAGE POINTS 

Symmetric free-free modes and frequencies were obtained while the vehicle was 
suspended in a horizontal and inverted position, with soft springs attached to the inter- 
stage points (Fig. 4-1). The shaker was mounted at the nose bulkhead and oriented 
to produce an exciting force in the transverse (±Z) direction. Thirty nine components 
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T-39 Fig. 4-1 Orbiter Model Suspended in Horizontal Inverted Attitude for 

Vibration Tests — Forward View 


T-40 Fig. 4-2 Orbiter Model Suspended in Horizontal Inverted Attitude for 

Vibration Tests — Aft View 
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T-41 Fig. 4-3 Orbiter Model Suspended in Vertical Attitude for 

Vibration Tests With Axial Preload 
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T-43 Fig. 4-5 Linkage for Applying Axial Preload to Fuselage During Vibration Tests 

of Orbiter Model Suspended in Vertical Attitude 
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T- 44 Fig. 4-6 Typical Strain Gage Installation on Orbiter Model Fuselage Side Wall 

Panels and Longerons 
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of Strain and Dial Gages for Static Tests 









Table 4-1 Summary of Measured Symmetric Free-Free Modes (Horizontal Inverted Orientation) 


Mode 

No. 

Frequency, Hz 

Mode Description 

1 

43.6 

Fuselage 1st Bending 

2 

52.7 (avg of 51.2 & 54.2) 

Wing 1st Bending (vs Payload Vertical) 

3 

58.2 

Wing 1st Bending (vs Aft Fuselage Vertical) 

4 

80.1 

Fin Fore-Aft 

5 

95.7 

Payload Fore-Aft (vs Fwd Fuselage- Fore-Aft) 

6 

104.1 

Wing 1st Torsion 

7 

107.0 

Payload Aft Vertical 
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Table 4-2 Variation of Fundamental Frequency with Axial Preload 


Axial Load, lb 

Measured Frequency, 

(a) 

Hz 

0 

43.67 

200 

43.86 

400 

44.05 

600 

44.25 


TT-9 a The axial preload is in addition to the model dead load 



Fig. 4-8. Location of Strain Gages - Section View 
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of acceleration were monitored during the test to obtain data for plotting mode shapes. 
Figure 4-2 shows the aft portion of the vehicle and some of the accelerometers. (Note 
the accelerometers bonded to the wing and fin ballast weight. ) Seven symmetric free- 
free modes and frequencies were measured as listed in Table 4-1. Plots of the 
measured mode shapes are given in Appendix C of this volume. 


4.3 MODE SURVEY WITH ORBITER MODEL SUSPENDED AT FUSELAGE NOSE 

For this series of tests the model was suspended from cables attached to the 
nose lugs (Fig. 4-3 and 4-4). Both static and vibration tests were performed with the 
Orbiter in this position. This test was initiated after the horizontal test at a point 
where it was perceived that the results did not correlate with the NASTRAN Model I 
analysis. At this point, the ineffectiveness of the skins was believed to be a major 
contributor to the lack of correlation. Consensus indicated that by preloading the 
vehicle a quantative value could be obtained for the skin effectiveness factors that 
should be applied in the analysis. Consequently, a simple two-tier "whiffle tree" 
arrangement was attached to the aft end of the model for applying axial load (see Fig. 
4-5). Loading in increments of 200 lb up to a total of 600 lb was applied with the fre- 
quencies being measured for each load level. These values are shown in Table 4-2. 
The small change in frequency due to the applied preload indicated that the 600-lb 
load was not sufficient to eliminate the initial imperfections. 

Static strains were also measured with the vehicle in the nose-up position. Strain 
gages were attached at 20 locations (Fig. 4-6 through 4-8.) (Note that Fig. 4-7 also 
shows the location of dial gages used in the horizontal static test. ) Strains were 
measured for increments of axial load, some of which are shown plotted in Figures 4-9 
through 4-13. As indicated in Fig. 4-8, gages were mounted on both the inside and 
outside of the fuselage skins at a given station. For example, gages No. 3 and 4, on 
the left side wall at Station 99, are attached to the outside and inside, respectively. 

The difference in the inner and outer gages, divided by two, is the bending strain 
which was then used to estimate the change in the central deflection of the panel for 
the specified load range. With increasing load the curves should become parallel. 

The largest difference between the two curves is shown on Fig. 4-12 for gages No. 7 
and 8. (Methods that were used for analyzing this data are discussed in Section 5. ) 
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T-47 Fig. 4-9 Load-Strain Curves for Side Wall Panels (Sta 99.0) 
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Fig. 4-10 Load-Strain Curves for Side Wall Panels (Sta 110.0) 
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LONGITUDINAL STRAIN, ji-IN./IN. 


Fig. 4-11 Load-Strain Curves for Bottom Deck Panels (Sta 110.0) 
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In Fig. 4-13, the strains are plotted for gages No. 1, 2, 17 and 18, which were 
attached to the fuselage upper longerons and door longerons at Station 116. An impor- 
tant point is that the fuselage and door longerons are not straining the same amount. 
This is because there is an ineffectiveness factor associated with the door longeron 
due to the manner in which the door is connected to the fuselage. 



<L 
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Fig. 4-13 Load-Strain Curves for Longerons (Sta 116.0) 
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4.4 STATIC TESTS 



Static testing was conducted with the model setup shown in Fig. 4-14 through 
4-18. The model was supported at the interstage points by fittings that rested on 
pylons. The pylons were approximately two feet high, and were bolted directly to 
the test floor. Loading was applied in increments by means of dead weights (see Fig. 
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Fig. 4-14 Details of Model Support for Static Tests 
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4-15 through 4-17). Deflections were measured by dial gages, most of which had an 
accuracy of ±0.0001 in. Influence coefficients were obtained from the measured load- 
deflection curves for comparison with calculated values. The tests and type of load- 
ing are indicated in Table 4-3; plotted test data is given in Volume IHA, Appendix 
A3. 

Figure 4-17 shows the aft portion of the model during the application of a longi- 
tudinal (+X) force to the OMS ballast center of gravity. Loads were applied to the 
ballast weights to investigate the local flexibility of the attachments of these relatively 
heavy masses. Figure 4-18 shows a closeup of the cabin ballast weight instrumenta- 



Fig. 4-15 Static Test Setup (Side View) 
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tion during static test. From a physical point of view, all static testing was kept as 
simple as possible using only available instrumentation. This was done to control both 
the cost and the time schedule. 

Influence coefficients as shown in Table 4-4 were obtained by averaging the 
loading and unloading plots for both downloads (-Z) and uploads (+Z) (Fig. 4-19 and 
4-20). Note that Figure 4-19 shows the deflections measured at two symmetrically 
placed gages (dial gages No. 1 and 2) during tests with the cargo doors off and load 
applied at the nose of the model. A similar set of load deflection curves is shown in 
Figure 4-20 for the mid-fuselage region. For the latter case, the dial gages for the 
+Z loading were moved aft because of interference of the gages with the loading wires. 

The average of gages No. 1 and 2 for loading and unloading in the -Z direction 

—6 —6 

at the nose with the cargo doors off (test No. 7) is 74.0 x 10 and 87.5 x 10 in. /lb. 

—6 

The average of these two values is 80. 75 x 10 . Similarly, for test No« 8 (+Z 

“6 -a 

loading at the nose), these values are 84.5 x 10 _ and 86.2 x 10 , and the average is 

—6 —6 
85.35 x 10 . The NASTRAN Model I analysis gave 55.4 x 10 . All seven numbers 

appear in Table 4-4, columns 16 through 22 for the 6 ^ coefficient. The average of 



T-54 Fig. 4-16 Static Test Setup (Front View) 
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T-55 Fig. 4-17 Apparatus for Loading and Measurement at Aft End of Fuselage During 

Static Tests 
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T-56 Fig. 4-18 Dial Gages Used for Measurement of Forward Cabin Ballast Deflection 

During Static Tests 
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Table 4-3 Static Tests 


Test 


Payload 


Load 

No. 

Loading Condition 

Doors 

Direction 

Application 

4 

• Mid Fuselage 

— Normal Symmetric 

On 

+Z 

FS 116 

3 

— Normal Symmetric 

On 

-Z 

FS 117.5 

14 

— Normal Symmetric 

On 

-z 

FS 117.6 

5 

— Normal Symmetric 

Off 

+z 

FS 116 

6 

— Normal Symmetric 

Off 

-z 

FS 117.5 

15 

— Normal Symmetric 

Off 

-z 

FS 117.5 

35 

— Lateral Antisymmetric 

Off 

-Y 

FS 117.5 

17 

— Lateral Antisymmetric 

Off 

-z 

FS 117.5 

7 

• Fuselage Nose 

— Normal Symmetric 

Off 

-z 

FS 46.8 

8 

— Normal Symmetric 

Off 

+z 

FS 46.8 

30 

— Torsion Antisymmetric 

On 


BL± 11.5 

29 

— Torsion Antisymmetric 

Off 

< 

BL± 11.5 

9 

• Wing Tip 

— Normal Symmetric 

On 

+z 

FS 162 

13 

— Normal Symmetric 

On 

-z 

FS 162 

10 

- Normal Symmetric 

Off 

-z 

FS 162 

11 

— Torsion Antisymmetric 

On 

-K?X 

FS 162 

12 

— Torsion Antisymmetric 

Off 


FS 162 

18 

• Fin Ballast 

— Long. Symmetric 

Off 

+x 

Fin eg 

19 

— Long. Symmetric 

Off 

-X 

Fin eg 

KH 

— Long. Symmetric 

Off 

-X 

Fin eg 

K 

— Long. Symmetric 

On 

+x 

Fin eg 

mm 

— Long. Symmetric 

On 

-X 

Fin eg 

Hi 

— Lateral Antisymmetric 

Off 

-Y 

Fin eg 

22 

• OMS Ballast 

— Long. Symmetric 

Off 

+X 


23 

— Long. Symmetric 

Off 

-X 


25 

• Payload 

- Normal Symmetric 

Off 

+z 

FS 117.5 

26 

— Normal Symmetric 

Off 

-z 

FS 117.5 

28 

— Torsion Antisymmetric 

Off 


+Z@ 117.5, -Z @ 78.0 

32 

• Cabin Ballast 

— Normal Symmetric 

Off 

-z 

FS 55.4 
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Comparing 


all the test values, both download (-Z) and upload (+Z), is 83.05 x 10 

the analysis with this averaged test value gives 

6 n (Model I) „ 4 

= 22ii_ = 66 . 6 % 

11 (Test) 83.05 (4-1) 

For the mid-fuselage region, using test No. 5 and 6 the same type comparison 


gives 


« 33 (Model T) 35>61 


= 78 


(4-2) 


33 (Test) 45.87 

Both of the previous values are for the condition with the cargo door off. Comparing 
the analysis values for the door-on condition with the door-off condition 

31.0 


i (Model I - Door On) 
_oo_ 


a 33 (Model I - Door Off) 35. 61 


= 87 


and with comparable test data 


* 33 (Test - Door On) 
i 33 (Test - Door Off) 


40.7 

45.87 


=88.5 % 


(4-3) 


(4-4) 


which indicates about the same percentage effect of the door. In the analysis the skins 
were assumed to be fully effective, while Equation 4-2 indicates that the test model 
fuselage inertia is only 78 % of the analytical value. If 

Ip^ = analytical value of the moment of inertial of the fuselage 

I D ^ = analytical value of the moment of inertia of the doors about the com- 

bined neutral axis 


FT 


DT 


= test value of the fuselage inertia 

= test value of the inertia of the doors about the combined neutral axis 


then, from the above ratios, 


FT 


= 0.78 


FA 


FA 

T FA + T DA 


(4-5) 


(4-6) 


■ = 0. 87 
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i 

FT 

= 0.885 

X FT + I DT 


(4-7) 


Combining Equations 4-5 through 4-7, then 


DT 


■= 0.694 


DA 


(4-8) 


In Model I, the door longeron is not fully effective because it is attached to the 
fuselage at only four locations. In Subsection 5.12, the effectiveness factor asso- 

2 

ciated with Model I is estimated to be 0. 75. This factor being defined by I D = kA^d , 
where k is the effectiveness factor, A^ is the area of the door longeron, and d is the 
distance from the door longeron to the combined neutral axis. This means that the 
contribution of the door longeron to the combined inertia is only 52 % of the conven- 
tional section property value (0. 75 x 0.694 = 0. 516). The loss in effectiveness given 
by Equation 4-8 is due to the fact that Model I does not account for local bending or 
shear deformation in the straps and door longerons. 
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5 - ANALYTICAL MODEL INVESTIGATION 

5.1 BACKGROUND 

Section 5 contains a series of separate analytical model investigations conducted 
at Grumman. The relationship between these separate subtasks is illustrated in Fig. 
5-1. After evaluating the generalized mass and stiffness for the five Orbiter sub- 
structures (Subsection 5.2), a number of possibilities were suggested for explaining 
the differences between the Model I analysis and the test results. Consequently, 
separate studies were initiated to: 

• Investigate the adequacy of the CQDMEM2 element used in a major portion 
of the finite element model (Subsection 5.3) 

• Establish a simplified independent beam model of the fuselage and calculate 
frequencies (Subsections 5.4 and 5.5) 

• Evaluate the Guyan reduction procedure (Subsection 5.6) 

• Investigate the modeling in the forward region of the fuselage (Subsection 5. 7) 

• Evaluate the stress distributions due to 1 g and dynamic loads and check these 
stresses against linear buckling values (Subsections 5.8 and 5.9). 

While these studies were being performed, additional tests were made at 
Langley (Subsections 4. 3 and 4.4). These tests, especially the static test, confirmed 
some of the suspicions concerning the effectiveness of the fuselage skins and cargo 
door longerons. Test results lead to additional analytical tasks at Grumman to: 

• Study the effect of initial imperfections on the stiffness characteristics of 
panels loaded in tension and compression (Subsection 5.10) 

• Evaluate the test data to determine experimental effectiveness factors and 
estimate the amount of initial imperfection that existed in the l/8-scale 
model panels (Subsection 5. 11) 

• Analytically study the effectiveness of the cargo door attachments (Subsection 
5.12). 

While these studies were being performed at Grumman, additional analytical and 
experimental work was performed at Langley. This effort included: 
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• A study of the effects of varying the panel aspect ratios in a finite element 
model on the predicted frequencies 

• Examining the wing carry-through region and studying the local effects of the 
cutouts 

• Investigating the wing/fuselage shear ties 

• Evaluating the effects of removing the direct stress capability in the fuselage 
panels on the frequencies of the Orbiter 

• Studying various methods for modeling the cabin ballast 

• Analytically investigating the fin-fuselage connection 

• Investigating the effects of initial panel imperfections on the stiffness and 
frequency characteristics of a cantilevered beam constructed of bars and 
shear panels. (This investigation included the fabrication and testing of a 
separate built-up cantilever model that resembled the fuselage sidewall. 
Provision was made to "force in" known initial imperfections in the panels. ) 

The results of all these separate investigations led to a set of remodeling guide- 
lines for establishing Model II. These guidelines are discussed in Section 6. 

5.2 SUBSTRUCTURE CONTRIBUTIONS TO THE GENERALIZED MASS AND 
STIFFNESS 

The contributions of each of the five major components to the generalized mass 
and stiffness were calculated and used to evaluate the Rayleigh quotient: 


^2 * IM T [*j] l»l) K g 

Z l*i ) T C M i] IM Mg 


(Fuselage) 

(Wing) 

(Doors) 

(Fin) 

(Payload ) 

6299.78 

0.748 + 

0. 101 + 

0.130 + 

0.015 

+ 0.006 

0. 05629 

0. 548 + 

0. 156 + 

0. 024 + 

0.224 

+ 0. 048 
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{ *i } = portion of modal vector for the ith component 
[K. ] = stiffness for ith component 
[Mj] = mass matrix for ith component. 

(jJ 

The frequency is f = 5 — = 53. 2 Hz. 

The terms within the brackets are the values calculated for the fuselage, wing 
doors, fin and payload, respectively, for the first symmetric mode. These values 
are arranged in tabular form in Table 5-1. Note that the fuselage accounts for ap- 
proximately 75 % of the generalized stiffness. The fuselage weight of 100. 91 lb in- 
cludes the OMS ballast of 26.15 lb and the fuselage cabin ballast of 29. 2 lb per half 
structure. 


Table 5-1 Orbiter Substructure Contributions to First 

Symmetric Modal Mass and Stiffness (Model X) 


Substructure 

Weight 

Modal Stiffness 

Modal Mass 

lb 

% 

Value 

% 

Value 

% 

Fuselage 

100.91 

48.34 

4710.25 

■za 

.03087 

54.8 

Wing 

32.86 

15.74 

635.15 

WSSm 

.00880 

15.6 

Door 

6.52 

3.13 

822.34 

Mil 

.00131 

2.4 

Fin 

3.85 

1.84 

97.03 

1.5 

.01260 

22.4 

Payload 

64.62 

30.95 

35.01 

0.6 

.00271 

4.8 

Total 

208.76 

100.00 

6299.78 

100.00 

.05629 

100.00 
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The fin, which weighs only 2 % of the total vehicle, accounted for 22 % of the 
generalized mass. On the other hand, the payload weighs 30 % of the total and ac- 
counted for only 5 % of the generalized mass. 

Comparisons of the analytical Model I modes with test modes are presented in 
Section 7. For the first mode, both test and analytical values compare rather well 
except for the motion of the fin ballast. For this DOF, the ratio of the test value to 
the analysis value is approximately 


0 Test 
<t> Analysis 


1.25. 


(5- 
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An estimate on the analytical frequency that accounts for the additional flexibil- 
ity of the fin can be obtained by using the Rayleigh quotient with a modified mode 
shape that reflects the test value. That is 


*' = * Analysis + A *. f®- 4 ’ 

Since all points are in fair agreement, except the fin ballast, the only quantity in the 
quotient that needs to be modified is the fin generalized mass 

( M G ' ) Fin = ( °* 224) (1 ° 25)2 = °* 350 * (5-5) 

The generalized stiffness for the fin is small compared to the total and, hence, any 
small additional stiffness due to the modified mode shape can be neglected. An 
estimate on the frequency is thus 


l_ / 6299/78 \ 
27TV0. 05629 / 


1/2 


( 1.0 - 


1/2 


0.224 + 0.350 


= (53.2) (0.89) 


1/2 


= 50.2 Hz. 


(5-6) 


This value is still much higher than the measured value of 43.6 Hz and indicates 
that the "softness" of the fin supports is not the sole contributor to the lack of corre- 
lation. It might be suspected, however, that additional fin "softness", which had not 
been included in Model I, could account for approximately 30 % of the discrepancy. 

If Ul is a mode shape it must satisfy the equation 

[K]{*} -oj 2 [M]{*} =0. (5-7) 

Differentiating this equation with respect to a design parameter, p. (either a stiffness 
or mass parameter) 

ft ] 1 - 1 * M ftl 


-% Ml-' 

[M]{|4j = 0 (5-8) 
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Multiplying Equation 5-9 by { 0 } , and transposing Equation 5-7, it can be seen that 

the last term in Equation 5-9 vanishes. Rearranging the result, 


dcj 

dP { 


M 


dK 

d p i 


{#} 


M [M] U> 


M t |^1 { 0 } 

_</p _ 


(5-10) 


If p. is restricted to stiffness changes only, then; 


doj 2 _ 

apT" 


U} T [m] m . 


(5-11) 


Equation 5-11 provides a means for judging the change in the frequency due to 
changing the stiffness of a given design parameter, p., If the design parameter p. is 
taken to be a representative stiffness quantity of a given component (for example, a 
scalar multiplying the moment of inertia distribution of the fuselage) then the terms in 
the numerator of the Rayleigh quotient represent the relative values of the derivatives 
with respect to the scalar p.. Thus, if the inertia distribution of the fuselage is 
scaled by Pp, then an estimate on the new frequency, is 



where represents the old frequency. Once a change is made, the mode shape 
will of course change, and new values of the derivatives should be calculated. For 
the particular set of design parameters represented by Model I, however, a change 
in the fuselage inertia will produce the largest change in frequency. The next item 
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of most importance is the cargo door stiffness. In accordance with these observa- 
tions, separate investigations were initiated to study the fuselage El distribution, the 
adequacy of the finite elements used in the fuselage, and the finite element represen- 
tation around the cabin and OMS ballast. These investigations will be discussed in 
the next subsection. 

5.3 ADEQUACY OF THE CQDMEM2 ELEMENT 

Correct behavior of a finite element under in-plane bending is a rather severe 
test for adequacy of the element. A comparison of different NASTRAN elements 
under pure bending is given in Section 15 of the NASTRAN Theoretical Manual. For 
convenience, this data is reproduced in Fig. 5-2. Note that the CQDMEM2 element 
is 71 % too stiff ( q gg 5 -1 = 0. 71) using a Poisson's ratio of 0. 3. The CQDMEM 
element is much worse, and the CQDMEMl is not much of an improvement over the 
CQDMEM2 element. The results of "stacking" four CQDMEM2 elements to form a 
beam segment was discussed in Section 3. The effects of adding CROD elements to 
account for flange material was described, and it was noted that the behavior ap- 
proaches beam theory with an increasing ratio of cap material. The neutral axis of 
the l/8-scale model fuselage is not at mid-height, hence these curves do not strictly 
apply. 


POISSON’S RATIO - 0.3 

REF: NASTRAN THEORETICAL 
MANUAL 


A0/M 

Beam Theory 

QDMEM 

QDMEMI 

QDMEM2 

8/El 

1.000 

.224 

.675 

.585 


T-59 

Fig. 5-2 Comparison of Flexibilities for a Pure Bending Couple Applied to Square Membrane Elements, 8 = h 
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Figure 5-3 illustrates a typical fuselage side wall section. It contains two 
CROD elements , the upper representing the upper fuselage longeron, and the lower 
representing the lower longeron and deck. The side wall was modeled with five 
CQDMEM2 elements having a thickness of 0. 020 in. The moment of inertia, cal- 
culated from section properties, is 33. 0. The ratio of cap to web inertia being 


T 4 _ cap _ 22.87 _ „ 

1 U 105 ' ‘ 


(Note that the neutral axis is located 4. 78 in. above the lower cap. ) The heights of 
the five panels agree with the location of node lines at fuselage Station 116. A length 
of 10 in. was used for the composite segment, which is slightly larger than the maxi- 
mum distance between frames in the fuselage. The sides of the element were con- 



I (SECTION PROPERTY) = 33 IN. 4 
I ( EFFECTIVE) * 33.4 IN. 4 

T-60 Fig. 5-3 Fuselage Side Wall in Pure Bending 
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strained by MPC's to produce .pure bending, as shown in the figure. These calcula- 
tions yielded an "effective inertia" value of 33.4. The error is, therefore 


33.4 - 33.0 
33.0 


= 1.21 %. 


For shorter length segments the error would be less (as is indicated by the curves 
in Fig. 3-4). Thus it was felt that the idealization should work rather well, at least 
for predicting the frequencies of low modes. For high fuselage modes the number of 

transverse station cuts might not be sufficient, because of the more abrupt variations 
in the moment diagram for higher modes. However, before these effects become 
important, it was believed that other low mode effects not already included would be- 
come more important. 


In all of this discussion it has been assumed that the section is fully effective - 
that is, non- buckled. 

5.4 FUSELAGE El CURVES 

A beam equivalent model of the fuselage was created and analyzed by NASTRAN 
to check the validity of the finite element model. The simplified model is shown in 
Fig. 5-4 with section properties and lumped weights indicated at each station. The 
moments of inertia were obtained from section cuts and from deflections that were 
calculated using the finite element model. The fuselage vertical deflections (from 
the NASTRAN analysis) for a unit moment applied at the fin ballast are plotted in 
Fig. 5-5. The structure was supported at the interstage points (FS 68.25 and 166.5). 
As shown in the figure, a correction was made to account for the non-zero displace- 
ment at Station 68.25. This is because the forward support is on the centerline; 
hence, the local deformation of the frame produces a linear shift in the side wall dis- 
placements. These corrected deflections were then used to calculate the second de- 
rivatives by the formula 


yi" = ^ ( y i - 1 " 2y i + y i + 1 ) . 


The values of El were then computed from (EI)^ 


M. 


where is the moment at station i due to the unit moment at the fin ballast. The 
inertia values that were obtained are shown plotted as curve "B" in Fig. 5-6. Curves 
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"A" and "C" were obtained from section cut calculations. Curve "A" assumes the 
door longeron to be fully effective, while curve "C" does not include the door longe- 
ron. 

The position of the neutral axis was located by using the horizontal displace- 
ments calculated from the finite element model. These displacements were used to 
calculate average strains from which the neutral axis was located by assuming that 
plane sections remain plane during bending. The position of the neutral axis, plotted 
in Fig. 5-7, compares favorably with that obtained from section cuts. The location of 
the neutral axis does not appear to be influenced greatly by the wing except for a 
small effect near the aft end. 

A similar procedure was used to determine the bending properties of the aft end 
of the fuselage and the fin. Figure 5-8 shows the aft fuselage and fin stub with a unit 
horizontal load applied at the fin ballast (Z = 88.5). The horizontal displacements, 
shown in Fig. 5-9, were used to calculate the equivalent moment of inertia distribu- 
tion (also shown in Fig. 5-9). An elevation of Z = 51.5 was used as a zero reference 
for measuring all displacements. 

5. 5 DYNAMIC ANALYSIS OF BEAM MODELS 

The El curves discussed in Subsection 5.4 were used to create a stick model of 
the fuselage (Fig. 5-3). The model contained 32 nodes and 31 beam elements. The 
value of I that was used corresponded to curve "C" in Fig. 5-6 (section cut data - no 
door longeron). The neutral axis was assumed to be at a constant WL of 51.5. All 
mass data for the fuselage was lumped at the 32 node points. The wing, fin, payload, 
and cargo doors were not included. The MPC's were used to constrain the model 
such that no axial deformation was considered. The results for the first three sym- 
metric free-free modes for the beam model are shown in column "A" of Table 5-2. 
Results are also shown from the Phase I finite element analysis for the Model I 
fuselage substructure. The first mode of the beam model was high by approximately 
11 %. The beam model was modified to include axial deformation and the MPC's 
were removed. The results of this change are shown in column "B". Note that the 
effect is small on the first mode. Next, the location of the neutral axis was modified 
in the cabin region of the fuselage, and a more accurate distribution of the cabin bal- 
last was used. The results for this change are shown in column "C". Finally, shear 
deformation was included by using the fuselage side wall as the shear area. The re- 
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Table 5-2 Results of Fuselage Beam Model Analyses 


Mode 

NASTRAN 
Finite El. 
Model 1 
(Fuselage) 

NASTRAN Beam Models 

A 

Using MPC's 

B 

No MPC's 

C 

Neut. Ax. 8i Mass 

D 

Shear Defl. 

1 

62.2 

69 

68.3 

64.3 

61.2 

2 

129.9 

142.2 

134.6 

138.1 

115.8 

3 

191.3 

265.4 

244.1 

230.5 

189.9 

TT-12 







suits of this change are shown in column "D" of Table 5-2, and in Fig. 5-10. In this 
case, the first mode differed from the finite element model by only 1.6 %. The third 
mode differed by only 0. 7 %, thus showing the effect of including shear deformation 
in calculating higher modes. 

At this point in the analytical investigation, it became clear that no major er- 
rors had been made in the NASTRAN Model I analysis of the fuselage. The discrep- 
ancies between analysis and test were due either to error in the testing, or to a 
fundamentally different behavior between the 1/ 8-scale model and the basic assump- 
tions used in establishing the finite element model. 


STATION 


+2 
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Fig. 5-10 1/8-Scale Shuttle Beam Model - First Three Modes for Fuselage Symmetric Case 
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Reproducibility of the test results for horizontal and vertical dynamic tests 
(Subsection 4.2 and 4.3) eliminated any suspicion of the testing procedures. Unfor- 
tunately, the schedule did not permit dismantling the l/8-scale model and dynam- 
ic test of the individual components. Thus, no test data existed for the fuselage sub- 
structure to compare with the analytical results just discussed. Consideration was 
given to extending the beam model to represent the entire Orbiter. However, time 
required to perform other analytical investigations prevented this from happening. 

The advantage of a simple Orbiter model over the extensive finite element model 
would be that basic parameters could be investigated more easily to tune the model to 
the test results. Derivatives of frequency with respect to various stiffness param- 
eters (such as the moment of inertia of the fuselage) could be easily obtained. This, 
of course, could be obtained by using the finite element model and slaving a large 
number of elements to represent a single parameter. 

5. 6 EVALUATION OF THE GUYAN REDUCTION PROCEDURE 

The Guyan reduction involves an approximation in establishing the transforma- 
tion that is used to reduce the number of DOF's in a problem. Initially, the inertia 
forces at the DOF's to be eliminated are neglected. This allows formulation of a 
transformation, based on statics, that expresses the DOF's to be eliminated in terms 
of those to be retained. The reduction of the mass and stiffness matrices follow from 
virtual work considerations. If harmonic motion is assumed then the equations of mo- 
tion are 

{f} = [k]{a} -CU 2 [m]{aI (5-12) 

where 

{ F } = a vector of external forces 
{ A } = a vector of node displacements 
[k] = stiffness matrix 
[m] = mass matrix 
u = circular frequency. 

Partitioning Equation 5-12 into parts yield 



(5-13) 
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where the subscript "1" refers to DOF's to be eliminated, and "2" refers to DOF's 

2 2 

to be retained,, Neglecting the inertia terms, - cj and - cu M^, leads to 

{F x } “ +[K 12 ]{A 2 }. (5-14) 

If the external forces { F^} are zero, 

( A l) = - [ K 12](M- 

Using this relationship, the transformation, [t] , can be established 

IM W 1 [k 12 ti 

a 2 'J = — y — i A 2rW{*2i <®- 16 > 

or 

{ A) = [t] (A r | 

where the subscript "R" indicates reduced degrees of freedom. 

From virtual work considerations 

(F R ) = [T] T |f) . 

Combining equations 5-12, 5-17 and 5-18 yields 

(F R ) = [t] t [k] [t]{a r | -V 2 [t] t [m] [t]{a h | 

-[[Kr] -" 2 [Mr^rI < 5 - 19 ) 

where [K^ ] and [iVL^] are the reduced stiffness and mass matrices, respectively. 

/ 

The choice of which DOF's to eliminate becomes important because of the as- 
sumption involved in the transformation. To illustrate this point, and to discuss 
general methods for evaluating large systems, consider the simple two-DOF problem 
shown in Fig. 5-11. It consists of two masses and two springs. The values of the 
masses are 1 and 5, and the stiffness of the springs are both 1.0. Assuming har- 
monic motion, the equations of equilibrium are 


( 5 - 17 ) 


( 5 - 18 ) 
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/AE\ 

V~/ = ‘ 10 FOR EACH SPRING 
T-67 Fig. 5-11. Simple Two-DOF Problem 



If DOF No. 2 is eliminated, the transformation matrix, [T ] , is 



[T] IV 


The reduced mass and stiffness matrices for this case are 


[k^] = [t] T [k][t] = [l] 

[MjJ = [T] T [M][T] = [6]. 


Thus, 

( F R I = [1](V -O, 2 [6] |a R ) = 0 


and hence. 


oj 2 =-|= 0.1666. 

b 


IS 


On the other hand, if DOF No. 1 is eliminated, then the transformation, 

/ \ ~i 


H 

-[Kill _1 [* 12 ] 

A2 ) = 

I 



(5-20) 

(5-21) 

(5-22) 

(5-23) 

(5-24) 

(5-25) 

>]• 

(5-26) 

(5-27) 
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from which 

[k r ] = [ 1 / 2 ] 
[m r ] = [ 5 . 25 ] 

and 

oj 2 = 0.0952. 


The exact solution is 

ujj 2 = . 09502 and = 2. 10498 . 


The mode for the lower frequency is 


Ajl (0.5249 | 

*2) 1.0/ 


A method for evaluating the Guyan reduction procedure has been suggested 
(S. Goldenberg of Grumman) and consists of first evaluating the inertia forces asso- 
ciated with a unit value of the mode shape. Returning to the first case (DOF No. 2 
eliminated), 

{F r }= -J [M] [T] {*} (5-28) 



(5-29) 


2 

The scalar cu may be neglected since only the shape of the distribution is of interest. 
Using these forces the static deflection shape is calculated from 


or 


{a}=[k]“ 1 {Fj} 


(5-30) 



(5-31) 
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This static deflection shape is now used in the Rayleigh quotient 


w 2 {a} T [k] {A} 
{A} [M]{A}« 


(5-32) 


For the present case this gives u> 2 = 0. 09516. The static deflection shape | 

lo 6466) 1 11J 

is equivalent to (^) 

Repeating this procedure for the case when DOF No. 1 is eliminated yields 


u 2 = 0. 09501 




0. 5238 

H 


1 i . 


These results are summarized in Table 5-3. Although a rather extreme case, 
this shows that the results that are obtained by the Guyan reduction are heavily depen- 
dent on the DOF's being eliminated. Eliminating No. 2 is a poor choice, as would be 
suspected. Regardless of the choice, however, it is possible to judge the validity of 
the results by first calculating the inertia forces due to the ,, expanded ,, Guyan mode, 
and then using these forces to calculate static deflections. The frequency can then be 
obtained by using these static deflection shapes in the Rayleigh quotient. For the 
present example, the elimination of DOF No. 2 yielded an error in frequency of 75. 3 
%, which was improved by using the static deflection technique to being in error by 
only 1.47 %. The elimination of DOF No. 1 is in error by 1. 89 %, and this was im- 
proved to being in error by only 0. 01 %. For this simple case the exact solution is 


Table 5-3 Summary of Calculations for Simple 
Two- DOF Problem 


Calculation 


% 

Error 

Mode 

Shape 

(1) Guyan Elimination of No. 2 

0.1666 

75.3 

1 

, 1 

(2) Guyan Elimination of No. 1 

0.0952 

1.89 

0.5 

.1 

(3) Using Static Shape From ( 1 ) 

0.09516 

1.47 

0.5455, 1 

(4) Using Static Shape From (2) 

0.09501 

.01 

0.5238, 1 

(5) Exact 
TT-13 

0.09502 

— 

0.5249, 1 
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known, hence an error can be calculated. Of course, for large problems this is not 
the case. Thus it becomes desirable to have procedures that can estimate the error 
or follow the convergence. Rewriting Equation (5-12). 




K LI " u ^1 ; ^2 - w ^*12 

H 

hj 


."Sil -“ 2 “21 i ^2 -“ 2 “ 22 . 

u 


(5-33) 



(5-34) 


Setting the external force { J =0 leads to 

= - [ A u] " X [Ai 2 ] |a 2 ) (5-35) 

from which 

{-^-} = [-~ --r~- ]w * M U»> • < 5 - 36 > 

Using this transformation in conjunction with the full stiffness and mass matrices 
leads to 

[K r ] = [T] T [K] [T] (5-37) 

[M R ]=[T] T [M] [T] (5-38) 

where [K r ] and [M r ] are the reduced stiffness and mass matrices. Here, however, 
the transformation [t] is frequency dependent. Note that if u; is set to zero, 

A.. = K.. and A = K , and the transformation is identical to the transformation 

11 11 1 C* 1 Li 

used in the Guyan reduction. An iteration procedure can be used to evaluate the 
transformation [t] , assuming that only the first mode is of interest. Using this 
procedure on the previous problem where DOF No. 2 has been eliminated leads to the 
values shown in Table 5-4. Observe that convergence is occurring, although not as 
fast as the one-step static deflection procedure. 

A shortcoming of the methods just discussed is that they do not give an error 
bound on the original calculations. For example, in the static deflection procedure, 
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Table 5-4 Convergence of Frequency Using Frequency- Dependent 
Reduction Transformation for Simple Two- DOF Problem 


Cycle 

Assumed 

Mode 

Calculated 

o? 

Change in 

1 

0 

1. 1 

0.16 


2 

0.16 

0.20, 1 

0.135 

0.025 

3 

0.135 

0.333, 1 

0.109 

0.026 

4 

0.109 

0.4545, 1 

0.0968 

0.0122 

5 

0.0968 

0.5163, 1 

0.09504 

0.0017 

Exact 

TT-14 

— 

0.5249, 1 

0.09502 

— 


a new mode and frequency must be calculated to estimate the improvement In the 
original value. There Is no Indication of the closeness of this second value to the 
exact value. 

A procedure for predicting a frequency error bound has been suggested by 
I. U Ojalvo of Grumman (Reference 9-6). This method was used to determine the 
error bounds on the first three symmetric modes from the Model I fuselage Phase 1 
analysis. This was done to determine if enough dynamic DOF's (or the correct ones) 
were retained after the Guyan reduction to ensure that this was not the cause of the 
discrepancy between the analysis and test results. The error bound formula (derived 
in Appendix D of this volume) is 


!u} T [m] {u} >/, _ X + g \ 2 (5-39) 

{?} T [m] i*} \ X + e / 

where 

I u )=[ k m ]' 1 h 

and 

r ; -- ! K m] “(A + «) [m] {*) . 

Hero, { \ = expanded calculated mode obtained from 

(*)“ [ T ]l”al 

where |y>^ j is the mode for the reduced dynamic 
DOF system, and [t] is the reduction transformation. 
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Since •{ ? } is not exact, the full system of equations are not 
satisfied, that is { F} ^ 0. 

— — 2 

X = to = calculated eigenvalue (approximate) 

2 

X = cj = exact eigenvalue 

« = an increment added to the eigenvalue. This 
is used to convert the singular stiffness 
matrix [Kj to by 

[Km] “ [ K ] + « [m] . (5-40) 

The results of the error bound calculations are summarized in Table 5-5. 

Note that the computed first symmetric mode, which is the first fuselage bending 
mode, has an accuracy of greater than 95%. The test and the analysis results, 
however, disagreed by approximately 20% for the first fuselage bending mode. 

Thus, it was concluded that the Guyan reduction was not the cause of the discrepancy. 


Table 5-5 Results of Error Bound Calculation for First 

Three Symmetric Fuselage Substructure Modes 


Mode 

Mode Description 
(Frequency, Hz) 

Accuracy of Frequency 
(Greater than, %) 

1 

1st Bending (62.2) 

95 

2 

2nd Bending (129.9) 

82 

3 

3rd Bending (191.3) 

71 

TT-15 




5.7 FORWARD FUSELAGE INVESTIGATION 

The forward fuselage, from Station 46.75 to 87.5 shown in Fig. 5-12 and 5-13, 
was investigated to determine the sensitivity of the modeling of the ballast deck and 
surrounding structure on the calculated frequency. Symmetric modes were calcu- 
lated for three cases with the models cantilevered at Station 87.5. 

CASE 1 - The Case 1 analysis was used as a basis of comparison for Cases 2 and 3. 
The structure and dynamic DOF were identical to the forward portion of Model I. 
Figure 5-14 shows the first mode results. The following calculations were per- 
formed to verify that these results were representative of the first symmetric Model 
I Orbiter mode. Initially, the idealized ballast center of gravity and weight were 
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STA64 



BALLAST GENERALIZED MASS 


0 T M0 ■ 


3gg-j [28.8 (0.825) 2 + 1213 (0.0301 ) 2 ] 
0.051 + 0.003 = 0.054 


62.50 HZ 
<t> J M<t> • 0.0605 
" 9327 
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Fig. 5-14 Cantilevered Forward Fuselage (First Mode) 


determined (Fig. 5-12). Then the generalized mass of the ballast was calculated by 
scaling the motion of the center of gravity (Fig. 5-14). This was compared with the 
value obtained by using the Model I Orbiter mode after adjusting the mode shape to 
account for the relative motion at Station 87. 5. Figure 5-15 shows the forward por- 
tion of the fuselage and the Model I mode shape with that adjustment included. Figure 
5-16 shows the same mode shape, but with the undeformed fuselage translated and 
rotated so as to align the fuselage Stations at Station 87.5. From Fig. 5-14, the 
ballast generalized mass is 

M G = Wl [ 28 - 8 (0-825) 2 +1213 (0.0301) 2 ] 

= 0.054 

The first term in brackets is the contribution of the linear motion (0. 825) of 
the ballast mass (28. 8) while the second is that for the rotation (0. 0301) of the ballast 
inertia (1213). 
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T-72 Fig. 5-15 Forward Fuselage First Symmetric Mode - Model I Orbiter Analysis 


STA 87.5 



0 T M0 = (0.212) 2 + (0.00783) 2 

= 0.00335 +0.00019 = 0 00354 

SCALE CONVERSION FACTOR = - 15.6 TO CONVERT TO CANTILEVER 

' ' MODE SCALE (TEST RUN 1 ) 

.-. 0 T M0 - 0.00354(15.6) - 0.055 WHICH AGREES WITH TEST RUN 1 RESULTS 
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Fig. 5-16 Forward Fuselage First Symmetric Mode — Model I Orbiter Analysis Relative to Sta 87.5 (53.24 Hz) 
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Similarly, from Fig. 5-16, the ballast generalized mass is 

m g = 3^74 [ 28.8 (. 212 ) 2 +1213 (0.00783) 2 ] 

= 0.00354 


Adjusting this value for the scaling in the mode shape yields 

(0.00354) 1 0 .! 

\ 0 . 

which agrees with the cantilevered case and verifies the use of the cantilevered 
model to adequately represent the behavior of this portion of the fuselage. Results 
of the calculations for Case 1 are tabulated in Table 5-6. 

CASE 2 - This model is the same as Case 1, except that dynamic DOF's were re- 
tained at the corners of the ballast deck. The results (Refer to Table 5-6) indicated 
that additional DOF at the ballast deck will not alter the frequency of the first mode. 
The shape is plotted in Fig. 5-17. 

CASE 3 - The ballast deck was replaced with a rigid deck where the center of gravity 
and weight were determined from actual dimensions of the ballast (see Fig. 5-13). 
The DOF's at the center of gravity were retained with all deck ties described by 
MPC's. The shell structure from Station 46.75 to 64.0 was changed from CQDMEM2 
elements to an equivalent CROD and CSHEAR arrangement. These changes also did 
not produce an appreciable change in the first mode results. (Refer to Table 5-6 and 
Fig. 5-18.) 

The results of this investigation eliminated the ballast deck idealization and 
the selected DOF's around it as being a major cause for the discrepancy between the 
Model I analytical results and the test values. 


■y 


= 0.055 


Table 5-6 Summary of Forward Fuselage Analytical Investigations 


Mode 

Case 

Frequency, Hz 

Generated Mass 

Generated Stiffness 

1 

1 

62.50 

0.0605 

9327 


2 

62.39 

0.0606 

9308 


3 

62.82 

0.0587 

9140 

2 

1 

239.13 

0.0428 

96577 


2 

228.15 

0.0617 

126752 


3 

235.52 

0.0768 

168089 

3 

1 

331.74 

0.0517 

224781 


2 

301.28 

0.0668 

239450 


3 

316.63 

0.0324 

128251 
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5.8 STRESS DISTRIBUTIONS 
5.8.1 1 g STRESSES 

The 1 g stresses were calculated for the fuselage with the Orbiter supported 
horizontally at the interstage attachments. A moment diagram for the half structure 
(Fig. 5-19) was obtained by representing the weights as a system of 'lumped" 
concentrated loads. At Station 119. 0 the value of the 1 g moment is 650 in. -lb (for 
the half structure). 

NASTRAN Model I stresses were obtained by using the reduced mass matrix 
with a 1 g acceleration. These stresses are illustrated in Fig. 5-20 for a fuselage 
section at Station 119. 0. A summation of the cap loads at Station 119. 0 gives an 
internal moment of 630 in. -lb. This compares rather well with the value of 650 
in. -lb from the moment diagram, considering that the moment diagram is only an 
approximation. 


A comparison of the 1 g Model I stresses with simple (Me) values is shown in 
Table 5-7 for the top and bottom longerons. Two values of the moment of inertia 
Iyy have been used; one for a fully effective section which is representative of Model 
I; the second for a partially effective section which is believed to be more repre- 
sentative of the actual 1/8-scale model. The values for the NASTRAN Model I 
analysis are within an acceptable range provided by the Me check. Here, the values 
of M and I are questionable. 


5.8.2 DYNAMIC STRESSES 


Dynamic stresses were estimated by multiplying the analytical amplitude for 

the first mode by that obtained in the test. For harmonic motion the displacement, 

6 , at any station, X, is related to the amplitude, A by 
z z 

5 Z = A z sin cut (5-41) 


differentiating and substituting 


S =_ 

z <u2 


(5-42) 


In the test, the RMS value for the accelerometer at Station 125. 5 was 0. 412 g for the 
first mode which has a frequency of 43. 6 Hz. (The amplitude of the harmonic force 
was 3 lb applied in the transverse direction at the nose fititng). The maximum dis- 
placement S max (at Station 125. 5) is thus 
z 
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Orbiter Model 1-g Moments Diagram (Half Structure Values) 
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& 




MOMENT @ STA 119 ABOUT WL 51 .5 (COMPUTED FROM CAP LOADS) 


13.95 

(11) 

= 

153 

18.16 

(11) 

= 

200 

6.21 

(7.88) 

= 

49 

3.42 

(5.2) 

* 

18 

1.48 

(2.5) 

= 

4 

2.27 

(1.67) 

= 

4 

35.77 

(5.66) 

m 

202 


630 IN.-LB (FIG. 5-19 GIVES 650 IN.-LB) 
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Fig. 5-20 NASTRAN Model I 1-g Stresses and Cap Loads 
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, max w (l-^) P-412M386) , 0.003 in. 

Z (27 t x 43. 6) 

The strains in the top fuselage longeron were obtained by taking the difference 
of the node displacements for the first mode. These values were then normalized to 
a unit value of the mode shape at Station 125. 5 . The normalized strain is relatively 
constant over the mid-portion of the fuselage indicating that deformation has an 
approximately constant radius of curvature for this region. For Model I, the nor- 
malized strain was calculated to be 0. 0102 as shown in Table 5-8; a similar analysis 
performed later for Model II gave a value of 0.0112. Table 5-8 also lists other key 
information for Models I and II and the horizontal dynamic test. 


Table 5-8 Selected Analysis and Test Results for First Free-Free 
Symmetric Mode 



Model 1 

Model II 

Horiz. Dyn. 
Test 

Frequency, Hz 

53.2 

44.2 

43.6 

e x (Top Longeron) 

0.0102 

0.0112 

- 

6 Z (Sta 125.5) 

1.0 

1.0 

1.0 

5 X (Fin Ballast) 

3.36 

3.73 

4.0 

5 Z (Sta 116.0) 

1.04 

1.05 

1.0+ 

6 Z (Wing Tip, Sta 162) 

1.53 

0.396 

1.12 
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The maximum dynamic stress in the longeron can be estimated by using the 
analytical strain for a unit deformation and the test value for the amplitude 

£ = (0.0102) (0.003) = 30.2 X 10 -6 

hence, 

cr= E€ = (10.5 X 10 6 ) (30.2 X 10" 6 ) 

= 317 psi (top longeron) . 

For a 16-in. deep section, the neutral axis is located approximately 6 in. from the 
bottom longeron. Thus the bottom longeron stress is approximately 

cr=j^ (317) = 190 psi (bottom longeron). 
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In the horizontal dynamic test the vehicle was suspended in an inverted position 
from the interstage points. The total maximum stresses between Stations 119 to 
125.5 were thus estimated to be 

Total = 1 g + dynamic 

<x (bottom longeron) = -94 ±190 psi 
<r(top longeron) = +137 ± 317 psi. 

Hence, the maximum compressive stresses are, 

cr (bottom longeron) = -284 psi 
o-(top longeron) = -180 psi. 

5. 9 LINEAR BUCKLING ANALYSIS 

Elastic buckling stresses for the side wall (in both the inverted and right-side 
up conditions) and bottom deck panels are given in Table 5-9. Buckling loads for 
typical panels in other parts of the model are given in Table 5-10. The coefficients, 
k, used in the buckling formulas, are based on the assumption that the panels are 
simply supported. The buckling stresses are plotted in Fig. 5-21 and 5-22 together 
with the 1 g stresses from Subsection 5.8. The lowest buckling stress was found 
to be 167 psi for a 9^-xl2^-in. bottom deck panel. This panel size and shape, which 
is typical of most of the panels in the lower deck, is also the most critical from a 
margin-of-safety point of view (MS = 0. 412 for 1-g plus dynamic stresses). 

The 1 g stress in the central portion of the lower longeron was calculated to 
-94 psi (Subsection 5.8). The stresses within the bottom deck panels are somewhat 
lower than the longeron stresses because of local shear lag effects (see Fig. 5-20), 
however, the longeron stress may be used as a representative value. In Subsection 
5.8 the dynamic stress in the lower longeron was estimated to be 190 psi. This 
maximum amplitude value is almost constant throughout the central region of the 
fuselage. (The middle portion of the fuselage is basically a light beam connecting 
two heavy masses - one is the forward cabin region; the second is the aft thrust 
structure. This simple two-mass beam model has a constant bending moment along 
the beam for the first free-free mode). 

Superimposing the harmonic dynamic stress and the static 1 g stress leads to 
the diagram shown in Fig. 5-23. The dynamic stress of 190 psi amplitude oscillates 
about the static stress of -94 psi with the total value exceeding the elastic buckling 
stress of 167 psi for approximately 4/10 of each cycle. 
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Table 5-9 Elastic Buckling Stresses for Simply Supported Side Wall and Bottom Deck Panels 

m 1 - » - « « ^ k n 2 D ■ - a - » 


r 


“►j No \+ a ►] 



k * (- m + 1 

a Dm 


Side Wall 



Bottom Deck 


Frame 7 Sta 97 

♦ 

8 106.5 

9 116.0 

10 119.0 

11 125.5 

12 135 

13 141.75 

14 144.75 


Side Wall (Normal) 


(a/bl | k |o c . psi 


Bottom Deck 


IE23I 


Side Wall (Upside Down) 


(a/b) 


k 

o c . P» 

35.3 

533 

35.3 

533 


«c- P*' 


9.5 1.6 16 0.594 17.9 270 12.5 0.76 4.33 167 2.67 16 0.594 35.3 533 

9.5 1.6 16 0.594 17.9 270 12.5 0.76 4.33 167 2.67 16 0.594 35.3 533 

3.0 1.6 16 0.187 - - 12.5 0.24 19.4 755 2.67 16 0.187 - - 

6.5 1.6 16 0.406 23.6 356 12.5 0.52 5.95 230 2.67 16 0.406 39.5 596 

9.5 1.6 16 0.594 17.9 270 12.5 0.76 4.33 167 2.67 16 0.594 35.3 533 

6.75 1.6 16 0.421 23 347 12.5 0.54 5.72 221 2.67 16 0.421 38.8 585 


1.6 16 0.187 


12.5 0.24 19.4 755 2.67 16 0.187 


See Reference 9-8 
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Table 5-10 Elastic Buckling Loads 


Wing Cover 


Square Plate 


Fuselage Side 
Wall 


Fuselage Bot. 
Deck 


in. I in. 


0.020 

16 

8.7 

10? 

0.020 

10 

10 

10? 

0.020 

9.5 

17 

107 

0.025 

9.5 

12.5 

107 


1.835 


10? I 0.316 1.0 


CLAMPED 


SIMPLE SUPPORT 


381 66.4 


10.07 367 73.4 


162 55.2 


417 30.5 


4.03 

192 

33.4 

4.00 

146 

29.2 

5.50 

69 

23.4 

4.31 

155 

11.3 
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AVERAGE STRESSES IN BOTTOM DECK PANELS DUE TO DEAD WEIGHT 



Fig. 5-21. Comparison of 1-g and Critical Buckling Stresses for Bottom Deck Panels — Orbiter Suspended 
Inverted from Interstage Fittings 



Fig. 5-22 Comparison of 1-g and Critical Buckling Stresses of Fuselage Side Wall Panels - 
Orbiter Suspended Upside Down from Interstage Fittings 
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TENSION 

+96 
0 

STRESS. PSI 

-94 

-167 

-284 
COMPRESSION 

T _gi Fig. 5-23 Stress Variation in One Cycle for a 9-H x 12-16-Inch Bottom Deck Panel 
Located in the Mid-Fuselage Region 

These results should be viewed qualitatively. The elastic buckling stress is 
based on the assumption that the panels are simply supported. A fully clamped 
panel with an aspect ratio of 0.75 would have a buckling coefficient, k, equal to 11.69 
which is 2.7 times the simply supported value of 4.33. For this case, the buckling 
stress would be -450 psi, which is well above the maximum stress of -284 psi. The 
panels are certainly far from being clamped, and most likely the constraints are 
more represented by the simply supported condition. Any additional constraint will 
of course raise the buckling value. If the buckling value is exceeded, as indicated in 
Fig. 5-23, then theoretically the "effective" modulus of elasticity will immediately 
drop to one-half of the linear value. This assumes that the panel is initially flat. If 
the panel contains an initial imperfection, then the phenomenon is altered and the 
reduction of E to an effective E takes place more gradually and well before the elas- 
tic buckling load. This is discussed in more detail in the next subsection. Since 
there is ample evidence that the panels do contain initial imperfection (refer to Sub- 
section 5.11), any concern about exceeding the elastic buckling value is somewhat 

academic. Of more interest is the fact that lower deck strains are oscillating about 

“6 ""6 

a static strain of -9.0 X 10 in. /in. with an amplitude of approximately 18.0 X 10~ 

in./in. These values set a bound on the range of the "Effective" modulus that will be 
selected in later sections. 
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5.10 EFFECT OF INITIAL IMPERFECTIONS ON PANEL STIFFNESS 


Figure 5-24A illustrates a perfectly flat panel supported on all four sides and 

loaded uniformly on two opposite sides. The panel will behave in a linear fashion 

until the stress, cr , reaches cr where 

c 


_ k T? 

c 12 


( 1 - 1 / 2 ) 


fe) ! : 


(5-43) 


The term cr is the Euler buckling stress, and k is a constant that depends on 
c 

the aspect ratio of the plate. If the load is increased past cr the center portion of 
the panel becomes ineffective (Fig. 5-24B), and results in a non-uniform distribution 
of stress. If, as suggested by VonKarman, this distribution is replaced by the dis- 
tribution shown in Fig. 5-24C where ( ei is the effective width adjacent to each side, 

\ 2 


then Equation 5-43 may be modified to give 

2 


a e = 


kir ' 


12 (1- v) 


.ft! 


(5-44) 


where cr is the edge stress. (The rationale behind this is illustrated in Fig. 5-25 
where the uniform stress, cr , acts on an equivalent plate of width b . The ineffec- 

c c 

tive portion of the real plate simply ties the two strips together where each strip has 
deformed into a one-quarter sine wave.) In this expression, k is usually taken equal 
to 4.0 (long plate). Solving Equation 5-44 for b g , 

v _ /E 


b = 
e 


\/ 3(1- v 2 )J a e 


t = 1. 92iy 

V a e 


(5-45) 


If the loading is increased until the edge stress approaches the yield value, then for 

0 

aluminum (with E = 10.5 x 10 and “ 42,000 psi) this expression equals 


- -I QOC /1Q« 5 
1,9Z V 42, 


5 x 10 
000 


t =30.4 t 


(5-46) 


which is the expression that is often used to "lump" sheet material into adjacent 
caps. 


Equations 5-43 and 5-44 may be written as 

,2 k ir ^E 

a c b = 2 

C 12 (l-i/ ) 


t 


2 
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t-m Fig. 5-24 Plate Behavior Under Compressive Load 



Fig. 5-25 Von Karman Model of the Effective Width of a Post-Buckled Plate 




n 


or 

. 2 k 7T 2 E .2 
O D 7T t 

e e 12 (1- i/ 2 ) 

from which 

,2 .2 
cr b = a b 
c e e 


or 



(5-47) 


This is the commonly used expression for relating the effective width, b to 

6 

the total width, b. The effective width, b , is load dependent, and must be obtained 

v 

by an iteration procedure where it enters into the calculation of the overall moment 
of inertia. 


If the total load on the sheet is P, then 


P = o’ t b 

e e . 


(5-48) 


Substituting Equation 5-47, 


P=bt (o- c <r e ) 


1/2 


Defining the average stress by. 


A bt 


leads to 


o= ( cr o) 

' c e' 


1/2 


(5-49) 


The stress, o , may be related to the overall strain, E, by 

o = E« . 

e 

Here, it is assumed that the plate is loaded by keeping the edges straight. The 
strain, « , is thus equal to the end shortening divided by the original length. Sub- 
stituting this into Equation 5-49, 

<T- ( O- E) 1/2 . 1/2 

' C ' € 


(5-50) 


n 


Differentiating, 


do - 


_i, „ r .l/2 -l/2_ 

d« -2 (,t c E) ‘ " E 


and dividing by E = -j- gives 


/ dcr \ 1/2 

*_W_ i fe.\ 
■gT"' 2 ' \ £ / 


- 1/2 


E' 

E 


At buckling, e = e and « = E € ; hence 

V V V 


E_'__l_ 

E “ 2 . (5-51) 

A plot of the stress strain curve is shown in Fig, 5-26. The curve is linear until 
the Euler stress is reached; at this point, the modulus drops immediately to one- 
half of the linear value. With increasing load the modulus continues to drop, but not 
as drastically as the sudden change in going past the Euler stress. 

Gerard and Becker (Reference 9-7) have given formulas for the postbuckling 
behavior of elastically compressed flat plates, which are 


For: 



TF = 0 +(1-/3) — 
e °e 


where /3 = 0. 50 (When the unloaded edges are 

held straight but free to move 
laterally) 

0 = 0 . 746 (When the unloaded edges are held 

straight and restrained from moving 
laterally). 


For 


w > 


» 
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1/2 

— = 0. 19 + 0. 81 (— ) (5-53) 

°e VV 

(Long, simply supported plate - 
edge held straight, free to move 
laterally.) 

Differentiating Equation 5-52 

dcr 

— = /5 E . (5-54) 

For the case of no lateral restraint 0 =0.5, which agrees with the previous 
expressions using the Von Karman assumption; however, for lateral restraint 
/3 = 0.746, indicating that the elastic postbuckling stiffness is greatly dependent on 
the conditions of lateral restraint. 

If the plate contains an initial imperfection, then the behavior illustrated in 
Fig. 5-26 is altered. In Appendix E of this volume, a large deflection analysis for 
axially loaded plates with an initial bow is formulated. The formulation follows 
closely the analysis given by Timoshenko (Reference 9-8), the addition being to in- 
clude into the lateral deflection term the displacement function 

w =W cos Fcos^ (5-55) 

o o 2a 2b ' ' 

where W q is the maximum initial imperfection at the center of the plate. A com- 
puter program was written and used to solve the resulting equations for various plate 
aspect ratios and values of W q . The results are plotted in Fig. 5-27 through 5-30. 

In Fig. 5-27 compressive load-strain curves have been plotted for a square plate 
having various initial imperfections. The curves have been normalized to the elastic 
buckling load 

2 

P = 10.07 and critical strain 

c a 

, with v = 0.316 and E - 10 7 . 

The end load, P, has been obtained by integrating the nonlinear end stress distribu- 
tion while the strain, e , is the end shortening divided by the origianl length. Note that 
the effect of the initial imperfection is to translate the conventional post buckling 
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Fig. 5-27 Load - Strain Curves for Square Plate with Various Initial Imperfections 





FUSELAGE SIDE WALL PANEL TYPICAL ASPECT RATIO 


n 



LINEARLY ELASTIC 






FUSELAGE BOTTOM DECK TYPICAL ASPECT RATIO 



Fig. 5-29 Load vs End Shortening for Typical Bottom Deck Panel 




WING TOP AND BOTTOM COVERS TYPICAL ASPECT RATIO 160 ' 
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1/4 LINEARLY ELASTIC (W 0 - 0) -200-1 

T-88 Fig- 5-30 Load vs End Shortening 




n 


behavior into the elastic range (compare Fig. 5-26 and 5-27). For (y 2 -) equal to 
0.04, the slope of the curve at the critical strain is approximately 0.75; while for 
Pr 2 -) equal to 2.0, the value is approximately 0.64. The curves tend to become 


parallel at larger strains where the slope is approximately 0.63. 


In Fig. 5-28 through 5-30, load strain curves are plotted for a typical fuse- 
lage side wall panel, bottom deck panel, and wing cover. These curves extend into 
both the tension and compression range. Considering Fig. 5-28, which represents a 
17x9. 5-in. fuselage sidewall panel, observe that the initial tension and compression 
slope is approximately one-half of the linear value for equal to 1.0. On the 

compressive side the curves are closely banded for a wide range of initial imper- 
fections. For large compressive strains the slopes of these curves approach a value 
slightly under one-half. On the tension side the slopes of the curves approach the 
linear value, where the panels with smaller initial imperfection approach this value 
faster. The plots in Fig. 5-29 and 5-30 exhibit the same behavior, the shift in 
values being due to different aspect ratios. In Fig. 5-31, the results are cross- 
plotted for various aspect ratios. Two curves are shown representing the initial 
effectiveness factor and the limiting values that occur with large strains. As would 
be expected, the smallest effectiveness factor is associated with the panel having the 
largest aspect ratio (short, wide panel loaded along long edges). 

Section 6 contains a discussion on the use of these results in combination with 
some practical considerations to remodel the fuselage and wing structures. Basi- 
cally, the following effectiveness factors, e, were selected: 

• Fuselage: 

- Axial and bending behavior in bottom deck and sidewall e = 0.5 

- Shear effectiveness factor e = 2/3 

• Wing: 

- Axial behavior in top and bottom covers e =0.85 

- Shear effectiveness in covers e =0.85 

* $ 

- Shear effectiveness in webs e=1.0. 

* 

In Subsection 6. 1. 3, the rationale f<>r selecting the shear and bending effectiveness 
factors (not covered in this section) is discussed. Appendix F contains detailed de- 
scriptions and calculations of these effectiveness factors. 
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5.11 EVALUATION OF TEST DATA 


5.11.1 PRELOAD 

As indicated in Subsection 4.3, the preload had a minor effect on the frequency. 
Considering a compact section, the natural frequency, f, varies with the axial load, 
P, according to the formula (Reference 9-9), 



(5-56) 


where f is the frequency with no preload and P is the Euler buckling load, 
o c 



(5-57) 


If L is taken as 97 in. (distance between interstage attachments, also the 

approximate distance between node lines for the first mode) and I is assumed to be 

85.2 in. 4 (the fully effective value at Station 116) then P = 940, 000 lb. Passing a 

c 

parabola through this point and the calculated value of 53. 20 Hz under no load, leads 
to curve "A" in Fig. 5-32. If an effectiveness factor of one-half is used in the wall 
and bottom deck panels, then the value of I at Station 116. 0 is approximately 69. 0 

4 

in. , and P = 810, 000 lb. Passing a parabola through this point and the zero pre- 
c 

load frequency (obtained from the horizontal test) leads to curve "B" in Fig. 5-32. 
The average axial dead load is approximately 150 lb for the model half structure 


FREQUENCY, Hz 
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(see Fig. 2-9), and has been added to the preload to obtain the total value. The test 
values are shown plotted on Fig. 5-32, where it can be seen that they lie between 
curves "A M and "B", being very close to curve "B. " Analytical values are shown 
for curves "A" and "B, " and a comparison of these values with the test data indicates 
that a curve passed through the test data is drifting slowly from curve "B" (partially 
effective section) to curve "A" (fully effective section). Although one can question 
the accuracy of the test data for this type of plotting, and the validity of applying 
Equation 5-56 to this analysis, there is at least a qualitative trend to the results. 
Each test point having a different preload represents a section with a different 
effectiveness factor. Thus there is a family of parabolas bounded by the curves 
"A" and "B, " where for any given parabola passing through a test point the effective- 
ness factor is constant. 


Considering two parabolas "x" and 

"y", then, 


f 2 -! 2 



X ox 

\ ^cx/ 

(5-58a) 

«■-«■ 
y oy 

( 1+ ^) 

\ cy/ . 

(5-58b) 


The frequency at zero preload is proportional to the square root of the moment 

of inertia 


f ox= C lA 

(5-59a) 

f oy= C l/V 

(5-59b) 

while the Euler load is directly proportional to I 


P =C„I 
cx 2 x 

(5-60a) 

P = C 0 I 
cy 2 y 

• 

Therefore 

P -p ft) 
cx Vy. 

(5-60b) 

(5-61) 
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From Equations 5-59a and b, 



Substituting Equation 5-62 into 5-61, 


P = P 
cx cy 



2 


Substituting Equation 5-63 into 5-58a, 



Solving for 




(5-62) 


(5-63) 


(5-64) 


(5-65) 


If parabola "y" corresponds to curve 
through the test value at 900 lb, then, 


"B" and "x" corresponds to the parabola 


4 = < 44 - 25 ) 2 - fl 57 ooo < 43 - 60 > 2 

from which 

f =44.226 Hz. 
ox 

Substituting into Equation 5-62 and 5-61, 


I 


x 



I = 1.0289 I 

y y 


and 

P = 1. 0289 P =833,425 lb. 
cx cy ’ 

The 3% increase in the moment of inertia requires a 10% increase in the 
effectiveness factor. Values of I for a fully effective section, half effective, and the 
condition just calculated, are shown in Table 5-11. For all values, the ratio of I 
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Table 5-11 Fuselage Moments of Inertia for Various Effectiveness Factors 


Effectiveness 

Factor 

Door On 

Door Off 

^door off) 

■ 

^ partial 

I 

^partial* 

(Ifull) 

<Ifull> 

IMI 

1.0 

85.2 

1.0 

71.2 

1.0 

0.836 

0.50 

69.0 

0.81 

57.8 

0.81 

0.836 

0.55 

71.0 

0.835 

59.4 

0.835 

0.836 
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(door off) to I (door on) has been kept constant. The preload thus caused a minor 
increase in the effectiveness factor for the section. A substantial increase would 
have required a prohibitive preload. 

5.11.2 STRAINS 

Measured strain data is shown in Fig. 4-9 through 4-13. The location of the- 
strain gages is shown in Fig. 4-8. The strains for the gages mounted in pairs on 
the inside and outside may be split into axial and bending strains, « A and « B , 



(5-66) 



(5-67) 


where « and c 0 represent the strains on the extreme fibers in the middle of a panel 
X 2 

(see Fig. 5-33). Assuming plane sections remain plane before and after bending, 

the strain, « may be related to the curvature (see Fig. 5-34) by 
13 

dx (1 + € B ) = (R + 7> ) . 


therefore 


1 

R 


2c 
B 

t 



(5-68) 


If it is assumed that the initial imperfection is in the form of a (1 - cos) wave 
and retains this shape under dead load and preload (see Fig. 5-35), then the change 
in deflection caused by the preload is 

y - (l - cos -?p) (5-69) 


5-55 































































n 


where a Q is the change in the maximum deflection at the center of the panel. Differ- 
entiating and substituting into Equation 5-68 



(5-70) 

(5-71) 


This expression is somewhat conservative. If a sine wave were chosen 

(y = a sin ^Qr-), where the bow is "out" in one panel and "in" in adjacent bays, 
ox a 0 

then the value of a Q is twice that given by Equation 5-71. The ratios of — for 

various panels are shown in Table 5-12. Values range from 0.07 to 0.54. Recog- 
nizing that these values are conservative; that they only represent the change due to 
preload; and that the preload has not straightened out the panel (since the plots of 
the inner and outer gages have not become parallel in Fig. 4-9 through 4-12); it is 


safe to conclude that the initial imperfections are equal to or greater than the panel 
thicknesses for a large number of panels. 

5.11.3 INFLUENCE COEFFICIENTS 

The procedure for obtaining influence coefficients was discussed in Subsection 
4. 4. Table 4-4 is a summary of the influence coefficients obtained from plots of the 
loading and unloading curves in both positive and negative directions. Typical load 
deflection curves are shown in Fig. 4-19 and 4-20. (Additional plots are given in 
Volume mA, Appendix A3.) 


For the mid-fuselage (Station 116), the ratio of analytical and test influence 
coefficients for the door-off condition is 


6 33 ( Model J ) _ 35.61 

$ 33 (Test) 45.87 


78 % . 


This ratio is in good agreement with the ratio of I for an effectiveness factor of one- 
half compared with I for a fully effective section (Table 5-10). 


I (1/2 effectiveness factor) _ 69.0 _ 
I (fully effective) 85.2 
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The deflection coefficient ratio for the nose is 


« n (Model I) 
« n (Test) 


55.4 

83.05 


66.6 %. 


The larger disagreement for this case has been attributed to the additional flexibility 
of the forward mid-fuselage splice of the test model. 

The deflection coefficient ratio for the wing tip is 

6 (Wing Tip - Model I, Symm) _ 290. „ n n 

{(Test) “ 362 8U 

With the structure supported on the interstages, loads on the wing tip do not 
produce appreciable bending in the fuselage. Hence the 20% discrepancy between 
analysis and test has been attributed to the initial imperfections in the test model wing 
skins and the additional flexibility associated with cutouts in the wing carry-through 
structure. 


Deflection data for the fin is given in Appendix G of this volume. For a hori- 
zontal load at the fin ballast the ratio of fin rotations with respect to the fuselage 
interface is 

0 (Model I) 534. 0 x 10 -6 _ _ n . „ 

a = r - 30. 4 70 . 

(Test) 1751.0 x 10“ 6 


This large discrepancy has been attributed to the flexibility of the forward fin- 
fuselage attachment of the test model. 

The ratio of deflection coefficients for the payload with respect to the inter- 
stage attachments is 

{ (Model I) = 27.1 = 64 « 

{ (Test) 42.4 /0m 

This discrepancy has been attributed in part to fuselage effectiveness, but mostly 
to local flexibility associated with the payload support structure. An intermediate 

analysis correcting all items except the payload yielded a payload deflection coeffi- 

_6 

cient of 30.8 x 10 . This represents an increase of 13. 5 %. It is not as large as 

the mid-fuselage discrepancy because the payload supports are rather close to the 
interstage attachments; hence, less fuselage bending is involved. 
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5.12 EFFECTIVENESS OF CARGO DOOR ATTACHMENTS 
5.12.1 SYMMETRIC BENDING CASE 


The door shell is primarily ineffective for symmetrical loading because of the 
breather joints. However, the edge longeron, which is connected to the fuselage at 
four stations, is partially effective. The effectiveness of the door longeron can be 
determined through the use of some simple analytical models. The NASTRAN Model 
I static load analysis with a unit load at Station 116. 0 can be approximated with a 
simple beam having a unit load at its mid-point. The length of the beam is 97 in. , 
which corresponds to the distance between the interstage attachments (points of sup- 
port for the Model I static analysis and for the static test). Assuming the beam 
properties to be constant and equal to the values at Station 116.0, the deflection at the 
mid-point is 


6 = 


PL 
48 El 


+ 


PL 

4 A G 
s • 


(5-72) 


Substituting L = 97, P = 1, G = and E = 10.5 x 10^ yields 

«- (“P + -t) IQ' 6 (5-73) 

s • 

With the cargo doors off, the fully effective section properties are 

I = 71.4 in. 4 
A g = 0.66 in. 2 

where the shear area, A , is taken as the cross-sectional area of the side walls. 

s 

Using these values, 6 = 34.3 x 10"® which agrees rather well with the Model I 
analysis where 6 = 35.61 x 10"®. 

In Subsection 6. 1. 3 following it is reasoned that the fuselage skins are approx- 
imately one-half effective for direct stress, and two-thirds effective for resisting 
shear. For this case, the reduced section properties are 

4 

I = 57.8 in. (refer to Appendix F of this volume) 

A s = | (.66) = 0.44 in. 2 
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Substituting these values into Equation 5-73, 6 = 44.9 x 10~®, which is in 
rather good agreement with the static test result, 6 33 = 45.87 x 10 -6 (refer to 
Subsection 4. 4). 

Figure 5-36 shows the effect on the moment of inertia when the door longeron 
is added. Since the door longeron shifts the neutral axis by only a small amount, 
compared to the total height of the section (see Table 5-13), the following approxi- 
mate expression can be used for the moment of inertia: 


I = I F + 2 A D (d) 2 . 


Substituting this expression into Equation 5-74 and solving for A D 



(5-74) 


(5-75) 


A 0 (EFF. DOOR LONGITUDINAL AREA) 



LA = A + 2A d 

EA*Z = Ax c + 2 A d (d + c) 

N A. (FULL SECT) = Z = £A * 2 = A x c + 2 A p (d + c) 

LA A + 2 A d 

I (FULL SECT) = I = l F + A (Z - c) 2 + 2A D Id + c-Z) 2 

IF Z - c SMALL, Z = c 

APPROX. I (FULL SECT) I - l F + 2 A Q (d) 2 
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Fig. 5-36 Moment of Inertia for Fuselage and Effective Door Longeron at Sta 116.0 


Table 5-13 Moment of Inertia Using Exact and Approximate Formulas for Sta 116.0 



•f 

A 

a d 

C 

Z 

1 

1 Approx 

Exact 

Approx 

1 Exact 

100% Eff. Fus 
0% Eff. ad 

71.4 

1.766 

0 

4.95 

- 

71.4 

- 

- 

100% Eff. Fus 
100% Eff. ad 

71.4 

1.766 

.06 

4.95 

5.73 

87.1 

88.2 

1.02 

50% Eff. Fus 
0% Eff. a d 

57.8 

1.124 

0 

5.35 

- 

57.8 

- 

- 

50% Eff. Fus 
100% Eff. ad 

57.8 

1.124 

.06 

5.35 

6.47 

72.6 

71.9 

0.99 
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Equation 5-75 is approximate and yields results which are about 5 % lower than the 
correct values when the shift in the neutral axis is appropriately accounted for. 
Substituting into Equation 5-74 the following values: 

S = 31.0 x 10"® (Model I deflection with cargo door attached) 

Ip = 71. 4 in. 4 (Estimated fully effective Model I fuselage) 

2 

A = .66 in. (Estimated fully effective shear area) 

s 

d = 11.55 in. (Refer to Appendix F of this volume) 

yields A D = 0. 042 in. 2 . The actual door longeron area used in Model I is 
0.056 in. 2, hence, the effectiveness factor is 


0.042 

0.056 


75 %. 


The door longeron is not fully effective, even in Model I, because it is not attached 
continuously. 

The effectiveness of the door longeron in the static test can be estimated by 
substituting the following into Equation 5-75: 

_6 

6 = 40. 7 x 10 ($ 33 from static test with the cargo door on) 

I_ = 57.8 in. 4 (I for 1/2 effective skins) 

r 

A g = 0.44 in. 2 (2/3 effectiveness factor on shear area) 


d = 16.5 - 5.35 + 0.31 = 11.15 (refer to Appendix F; the addition of 0.31 

is to account for the offset between the 
fuselage and door longerons) 

2 

from which A Q = 0.035 in. . Hence the effectiveness factor for the door longeron 
in the static test is 


0.035 

0.063 


55 %. 


The area of 0.063 is a revised value for the door longeron (see Fig. 5-37). 


To verify these results a separate, small, simplified model (Fig. 5-38) was 
analyzed using NASTRAN. The model consists of an assemblage of CBAR elements 
which incorporate the changes that were eventually made in Model n. Figures 5-37, 
5-39, and 5-40 show the details of the revised door longeron, door strap, and 
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Fig. 5-37 Door Longeron Properties 
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Fig. 5-39 Door - Fuselage Strap 


Fig. 5-40 Door - Fuselage Connection 


door-fuselage connections, respectively. The results of the NASTRAN analysis are 
shown in Fig. 5-41. The center deflection 6 = 41. 24 x 10 -6 agrees closely wijjji the 
test value of 40. 7 x 10 -6 . Substituting the calculated deflection and the properties of 
the simplified model (6 - 41.24 x 10“®, Ip = 57.8 in. A g = 0.44 in. 2 and 
d = 11.61) into Equation 5-75 gives A D = 0.029 in. 2 and, therefore, the effective- 
ness factor for the door longeron is 


0.029 

0.063 


46 % 
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0.0039 IN. -LB 
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Fig. 5-41 Simplified NASTRAN Model Results 
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The reduction in the effectiveness factor is due to inclusion of the cargo door con- 
nection details. 

5.12.2 ANTISYMMETRIC BEHAVIOR CASE 

The torsional stiffness of the fuselage is greatly dependent on the local flexi- 
bility of the door attachments. If these attachments are rigid (for example, the door 
longeron is continuously connected to the fuselage upper longeron by a rigid shear 
panel), then the combined "U" shaped fuselage and circular cargo door form a closed 
cell. If, however, the connections are very flexible, then the behavior of the com- 
bined section will approach the behavior of an open channel and the torsional stiffness 
will drop significantly. To study this effect, a typical fuselage section was considered 
(Fig. 5-42). The factors k^, kg, kg, and k 4 represent thickness effectiveness fac- 
tors. In Model I, points "a" and "b" (Fig. 5-42) are common, and the door extends 
from point "b" on one side to point "b" on the opposite side as a segmented arc. In 
Model II, points "a" and "b" are separated by 0.31 in. on both sides, and straps were 
inserted to connect the cargo door to the fuselage at four discrete points along the 
length of the door. For Model I, = kg = k 4 = 1.0, and kg represents the effective- 
ness due to discrete attachments although the connection was made at the same water 
line. In Model n, k^ = 1.0, kg = k 4 = 2/3, and kg represents the additional flexibility 
due to separating the longerons by a distance of 0.31 in. and including bending and 
shear in the straps and longerons. 



Fig. 5-42. Torque Cell Representation of Typical Fuselage Section 
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For a closed cell, the constant shear flow, q, is given by 


_ = _T T 

q 2A Z 2A 
i 

and the rotation, 9 , due to q is given by 



where 

A = enclosed area 

tj = effective thickness of segment i 

S. = length of segment i 

G. = effective shear modulus 
L = length of torque cell. 


(5-76) 


(5-77) 


With the effectiveness factors, k, associated with the thickness and not the shear 
modulus. Equations 5-76 and 5-77 may be combined to give the rotation, 0, for a 
unit torque 


9 = 


L 

4 A 2 G 



(5-78) 


Substituting the values given in Table 5-14 gives 

,-12 I 2460 . 40 


9= 14.95 x 10 


1600 


1000 




(5-79) 


In this expression, L =103 in. (the length of the cargo door), and G = 4. 04 x 10®. 
As mentioned previously for Model I, k^ = kg = = 1. 0. The value of kg was ob- 


Table 5-14 Parameters for Simple Torque Cell Calculations 


Segment 

2 A; 

Si 

ti 

Si/tj 

1. Door Shell 

491 

39.3 

.016 

2460 

2. Interface Region 

7.75 

.62 

.016 

40 

3. Side Wall 

400 

32.0 

.020 

1600 

4. Bottom Deck 

407.75 

25.0 

.025 

1000 

Z 

1306.5 

- 

- 

- 
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tained by equating the strain energy of an equivalent shear panel to the strain energy 
of the axial loads in the door and fuselage longerons plus the energy in an additional 
shear panel required to extend the door region (see Fig. 5-43). 

Thus, 

u (Equivalent Panel) ~ U (Shear Panel) + U (Longerons). (5-80) 

The shear flow of q lb/in. is transferred across the door connection by four fittings 
(Fig. 5-44), each assumed to carry equal shear forces 


H 



25. 75 q . 


(5-81) 


The axial force in the longerons thus builds and drops in a "saw tooth" pattern, Fig. 
5-45. The axial energy, U, is, 



(5-82) 


Evaluating this integral for the total length, L. = 103 in. , for both the door and 
fuselage longerons gives 



(5-83) 


where 


Ap = door longeron area 
A F = fuselage longeron area. 

The strain energy for the shear panel is 



(5-84) 


Adding Equations 5-83 and 5-84 and equating to the strain energy of the equivalent 
panel 


1 , 

2 




(5-85) 


From which 



(5-86) 
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Fig. 5-43 Equivalent Shear Panel to Represent Energy in Door and Fuselage Longerons 
Plus Physical Door Shell Extension (Model I Simulation) 
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Fig. 5-44 Assumed Shear Transfer Between Door and Fuselage for Model I 
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Fig. 5-45 Axial Load Distribution in Door and Fuselage Longeron 
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With A =(0.31) (103) =32.6, t =0.016, A D = 0.056, A p =0.1, and-| =2.6, 
Equation 5-86 gives kg = 0.032. 


Returning to Equation 5-79 


9 = 14.95 x 10 


-12 


2460 40 

1.0 0.032 


1600 1000 
1.0 1.0 


(8-87) 


—6 

= 0.095 x 10 radians. 

For Model I, the relative twist between Stations 64. 0 and 166.5 is 0.10 x 10 - 
radians. Comparing these two results 


0 (simple torque cell) _ 0.095 _ 95 ^ (5-88) 

9 (Model I) °* 100 ’ ' " 

From test data, the relative twist between Stations 64. 0 and 166.5 is 
9 = 0.17 x 10 -6 radians. Substituting this value, and = 1.0, k g = = 2/3 into 

Equation 5-79 and solving for kg, yields kg = 0.008. This value is one-quarter of 
the value of kg predicted by Model I analysis. 

The Model n analysis described in Section 6, predicted a twist 9 equal to 
0.32 x 10 -6 radians (approximately twice the test value) for a unit torque applied at 
the nose fittings. The following analyses were performed to understand the deficien- 
cies that still lie in Model n. A simplified model (Fig. 5-46) was devised and analy- 
zed making the following assumptions. 


The shear flow, q (assumed to be constant), is transferred from the door to 
the fuselage through the four straps. Taking advantage of antisymmetry only half of 
the door region was analyzed. The straps are assumed to carry equal shear and are 
pinned on the fuselage side, although offset from the fuselage longeron, and moment- 
connected to the door longeron. The vertical interface forces on the door longeron 
caused by the "strap tie moments" were obtained by assuming the longeron to be a 
continuous beam on rigid supports. The results of the hand analysis are shown in 
Fig. 5-46. Note that the vertical interface forces on the side wall produced a minor 
change in the shear that was assumed in Equation 5-79, and hence was neglected. 

The kg value for Model II was obtained by equating the strain energy of an equivalent 
panel to the combined strain energy (see Fig. 5-47) of the following: 


5-68 



0.94 q 


BEGINNING OF DOOR REGION 
(STATION 64) 


c 


- 14 "- 

-q— 






20.8 q 


-L/2 = 515" 
24" 

q 


- 135 " *- 


25.75 q | 

f *0.66 q I 0.30 q 




208 q 


25.75 qj 


DOOR LONGERON 


28 

YO.29 




DOOR LONGERON LEVEL 

f 

FUSELAGE LONG. LEVEL 
1 


0.5" 

T~ 


0.30 q 

25.75 q 


0.26 q 


0.26 q 


28.3 q I 

—XX- 25.75 q 


\ 

0.30 q 


25.75" 


STRAPS 


026 q 


25.75 q 


1288 q 


> 0.30 q 
>TN 25.75 q 


1288 q l°' 26q 


£ 


25.75 q 


*0.57 q 1 .44 q 


0.28 q 1.51 q 2.57 q 

FUSELAGE LONGERON 


0.70 q 


■AXIS OF 
ANTISYM 


'■ f '■ 


1 87 q 027 q 



T-104 Fig. 5-46 Assumed Shear Transfer Between Door and Fuselage for Model II 
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L = 103 IN. 



T 
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T-105 Fig. 5-47 Equivalent Shear Panel to Represent Strain Energy in Straps and Longerons 


• Axial energy in the door and fuselage longerons, same expression as 
Equation 5-83 

• Bending energy in the door and fuselage longerons 

• Bending and shear energy in the straps (the axial energy is negligible). 
The bending energy in the longerons and straps is 



(5-89) 


(5-90) 
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Adding Equations 5-83, 5-89, and 5-90 yields 

2 


17-l . 

i- -i- 1 cr7»7c; a. 

( 3700 ^ 734 J H 73 V 

L\a d 

• 1 lottO T 

A P / 

vv t; vAI 


2 

25" 


2920 ~ 
A 

L s J 


(5-91) 


Equating this expression to the equivalent shear panel energy given by the left 
side of Equation 5-85 

(A/t) 

*2 = 


G_ / 2920\ G / 5775 

L G s \ A s) E \ A D 


. 5775 , 3700 . 734 A 1173 

a f *f ^ 


). 


(5-92) 


Substituting the following into Equation 5-92: 
G = G. 


E/G 

A 


D 

X F 

A 


s 

= 2.6 

2 

= shear area of straps = 0. 035 in. 

2 

= area of door longeron = 0. 063 in. 

2 

= area of fuselage longeron = 0. 16 in. (includes equivalent skin area) 

4 

= inertia of door longeron = 0. 004 in. 

4 

= inertia of fuselage longeron = 0,002 in. 

4 

= inertia of strap longeron = 0. 002 in. 

2 

= area of equivalent shear panel = (0.31) (103) = 32.6 in. 

= thickness of equivalent shear panel = 0. 016 in. , 


yields kg = 0.0024. Substituting this value of kg into Equation 5-79 with k^ = 1.0 
and kg = k 4 = 2/3 yields 

6 = 14. 95 x 10" 12 [2, 460 + 16, 660 + 2, 400 + 1, 500] 


— 

= 0.344 x 10 radians. 


(5-93) 


Comparing this value with Model II results 


6 (Simple Model) _ 0, 344 
9 (Model II) 0. 320 


1.07. 


(5-94) 


The rather good agreement lends credibility to the assumptions employed in the sim- 
plified model as shown by Table 5-15. Note that the kg term in Equation 5-93 con- 
tributes 72% to the twist, 6 . If Equation 5-92 is normalized, then 
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_ 2040 
*2 854980 


strap shear. 


[o.i 


098 + 0. 041 + 0. 016 + 0. 416 + 0. 165 + 0. 264 


axial, door longeron 

axial, fuselage longeron, 
bending, door longeron— 


bending, fuselage longeron, 
bending, strap 


(5-95) 


Here it can be seen that the bending terms dominate the denominator. If the 
inertia of the door longeron is increased because of including more effective shell 
material, then k 2 will be smaller and the calculated twist will be more in line with 
test data. 


Table 5-15 Torque Cell Effectiveness Factors and Rotations for a Unit Nose Torsion 



Modal 1 

Model II 


Test 

ki 

1.0 

1.0 


1.0 (assumed) 

k 2 (calc) 

0.032 

0.0024 


0.008 

^3 

1.0 

2/3 


2/3 (assumed) 


1.0 

2/3 


2/3 (assumed) 

0 (Analysis-Test) 

0.100 x 1 0' 6 

0.32 x 10* 


0.170 x 10* 

9 (Calculated) 

0.095 x 1 0' 6 

0.344 x 10* 


- 
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FINITE ELEMENT ANALYSIS FOR MODELH 


SECTION 6 
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6 - FINITE ELEMENT ANALYSIS FOR MODEL II 

6. 1 FINITE ELEMENT MODEL H 

6.1.1 GENERAL 

The second finite element model was constructed using information obtained 

from 

• Behavior of Model I (Section 3) 

• Experimental results (Section 4) 

• Separate analytical investigations (Section 5). 

Various methods are available for tuning the parameters of a finite element 
model so that the analytical results are forced to agree with test data. If the number 
of individual parameters (areas and gages) is too large, then it is possible to "slave" 
groups of parameters into a fewer number of unknowns. This was not done because 
on a typical large aerospace project, the finite element model is the chief source of 
a large variety of information - including data such as: static deflections due to 
applied loads; thermal distortions; internal loads due to flight maneuvers, landings, 
taxing and take offs; thermal stresses; and influence coefficients for flutter analysis, 
transient dynamic analyses, and aeroelastic corrections. An attempt to "tune" a 
finite element model to accurately predict one type of behavior might jeopardize its 
ability to predict other types of behavior, thus degrading the overall usefulness of 
the model. The alternative is to create several "tuned" models; clearly this is not 
desirable, at least in the design phase of a project. 

Note that at this stage, a good deal of thought was given to modifying the 1/8- 
scale experimental model. No amount of analytical effort can change basic undesir- 
able behavior; only by changing the actual design can the situation be rectified. After 
reviewing a number of proposed changes it was decided not to incorporate any 
physical model changes because modification would have produced too drastic a 
schedule slippage for the project to tolerate. 

With the foregoing considerations as prime factors, the following course was 
decided: 

• Not to change the l/8-scale model 

• Not to "tune" the analytical model to specific test data 
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• To establish a second finite element model that would recognize test 
results and incorporate the findings obtained in the separate analytical 
studies. 

In the next subsection overall remodeling guide lines are discussed. These 
are elaborated on in the ensuing detailed subsections covering modifications to the 
individual components. 

6. 1. 2 REMODELING GUIDE LINES 

The guiding rules described in the following paragraphs for modifying Model I 
were established and observed in developing the second finite element model (Model 

n). 

FUSELAGE - A bar and shear panel idealization was used in place of the grid of 
direct stress elements. The areas of the bars include actual areas plus lumped 
effective skin material. An effective shear modulus was used in place of the actual 
material value. This accounts for additional shear flexibility due to the presence of 
initial imperfections. A typical effective width and effective shear modulus was used 
throughout the main portion of the fuselage, instead of attempting to use separate 
values for each panel based on test or measurements. It was felt that this repre- 
sented a more realistic design approach. 

The flexibility of the forward fuselage strap (refer to Volume HI, Appendix A, 
pg Al-1; see part AD 383-526-41 on view C-C) was included by double noding 
Station 64. 0 and inserting an axial spring between the nodes. The spring constant 
was determined using test data. 

WING - The wing carry through was modified to more accurately account for cut- 
outs. A bar and panel idealization was used on the wing covers, where the cover 
panel effective widths were lumped into equivalent bar areas. An effective shear 
modulus was used in modeling the rib and spar webs. 

FIN - Springs were added between the fin and the fuselage to account for the 
flexibility of the connections in the l/8-scale model. The values for these springs 
was based on static test data. No attempt was made to refine the model in this 
region. Early in the project, but after the model was constructed, this connection 
detail was recognized to be rather poor from a flexibility point of view. An alter- 
nate design was developed but not implemented on the constructed model. 
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CARGO DOOR - The door longeron was shifted to the actual 1/8 scale model 
location. In Model I, the door longeron had been set at the same water line as the 
fuselage upper longeron. Additional modeling detail was also added to more accu- 
rately account for the flexibility of the door-fuselage connections. 

PAYLOAD - Springs were added at the payload supports to account for the local 
flexibility of the connections. These values were determined by test data. 

6.1.3 FUSELAGE 

In Subsection 5. 10 the effective width for uniformly strained panels containing 
initial imperfections was discussed. Time prevented the extension of this study to 
include bending and shear. For a fuselage section containing many longerons the 
uniaxial straining case is sufficient. The conventional procedure is to lump adjacent 
skin effective widths with the longeron areas using an effective width formula 



where tr,. is the Euler stress for the ith panel and o\ is the bending stress obtain- 

Myj 

ed from o\ = -y— . The moment of inertia, I, is dependent on the lumped areas, 

hence the procedure is iterative. Since the stress depends on the sign and value of 
the bending moment, the moment of inertia is load -condition dependent. When the 
constant n equals one-half, the formula is of the form suggested by Von Karman. 
Marguerre suggested using two-thirds for the constant. Many other formulas have 
been suggested for calculating effective widths; however, again they are based on the 
uniaxial case. Where there are relatively few longerons the conventional effective 
width concept must be modified because a large stress gradient will exist across a 
given panel. The mid-section of the l/8-scale model fuselage contains five longe- 
rons (two upper, two lower, and one on the center line). Thus, for the symmetric 
case a bending stress distribution exists in the sidewall skins; for the antisymmetric 
case, a bending stress distribution exists in the bottom deck. 

The following rationale has been used to develop an effective amount of 
material to be considered as working in bending. In Fig. 6-1, the effective width 
is illustrated for a plate subjected to uniform end shortening. The stress drops off 
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T-106 


Fig. 6-1 Stress Distribution Due to Uniform End Shortening 


rapidly a short distance from the edge. Figure 6-2 illustrates an assumed stress 
distribution due to bending strains (the ends remain straight after bending). Again 
the stress is assumed to drop off from the linear (-p) distribution such that the 
inner portion becomes ineffective in resisting bending. The non-linear stress 
distribution is replaced with two trapezoidal blocks each having a height such 

that the moment of the trapezoidal blocks about the neutral axis is equal to the 
moment of the nonlinear distribution. The moment of the inertia of the partially 
effective rectangular section is 



( 6 - 2 ) 


where d Q = ^b - b g j is the height of the ineffective material. Since the effectiveness 
factor is 


e = 



(6-3) 


then the ratio of the effective I to the fully effective value I is, 

e o 



= 1 - (1 - e) 3 . 
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T-107 Fig. 6-2 Assumed Stress Distribution Due to Linear End Bending 

This ratio is relatively insensitive for values of e above one-half as can be seen 
here: 


e 

W 

0.4 

0.784 

0.5 

0.875 

0.6 

0.936 

0.7 

0.973 


It was decided to choose a common effectiveness factor for both the sidewall 
and bottom deck. This was done to simplify the situation and also recognize the 
fact that the symmetric and antisymmetric cases will produce reversed behavior in 
the side wall and bottom deck. (For the symmetric case, the bottom deck is loaded 
axially and the side wall is in bending; for the antisymmetric case, the side walls are 
loaded axially and the bottom deck is in. bending). The common effectiveness factor 
leads to one overall model for representing both types of boundary conditions. 

From the analytical work discussed in Section 5 the common effectivness factor was 
chosen to be one half. Applying this factor to the center of the fuselage, and pro- 
jecting it along the length, resulted in the distribution of effective material shown 
crosshatched in Fig. 6-3 (Note that the wide crosshatch in the forward and aft regions 
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FUSELAGE BOTTOM DECK (t - 0.025") 




II 


indicates that no change was made between*Model I and Model n. These areas were 
left as an assemblage of direct stress elements). Since there is no center longeron 
between Stations 64 and 78, and 162 to 186.25, no effective material has been in- 
cluded in the bottom deck central strip. A transition in the sidewall effective 
material strips was provided between Stations 125. 5 and 141.75 to recognize the 
sidewall cutout where the wing is attached, and the start of the angle attaching the 
top of the wing to the fuselage at WL 51.5. 

Figure 6-4 shows a typical fuselage mid-section and illustrates how the effec- 
tive widths were lumped into equivalent caps. In the NASTRAN model, the actual 
longerons were kept separate from the equivalent bar areas designated a. , a 2 and a g 
to facilitate possible future adjustment. The values of these areas were determined 
by solving three equations: the first equated the I of the total equivalent section 
(Fig. 6-4A) to I of the partially effective section (Figure 6-4B) the second equation 
preserved the location of the neutral axis; and the third employed an assumption, 

a 2 = * a l + a 3*‘ 

This assumption is based upon an even proportioning of the effective width material. 
The lateral moment of inertia, I , of the equivalent section was then checked against 
the I zz of the partially effective section. For an effectiveness factor e = 0. 5 this 
disagreement is approximately 10%. 



A. EQUIVALENT FUSELAGE SECTION B. PARTIALLY EFFECTIVE FUSELAGE SECTION 


T ' 1 09 Fig. 6-4 Method for Lumping Effective Panel Widths into Equivalent Caps 
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If an attempt is made to preserve the I , in lieu of the assumption that was 

employed, then a 2 becomes much larger and a g drops significantly. The results of 

three methods for calculating a^, a 2 , and a^ are shown in Table 6-1. It was felt 

that to preserve both I and I required too drastic a shift in material; since this 

yy zz 

approach would introduce extraneous local shear lag effects, it was not used. Curves 
and tables in Appendix F of this volume give section properties at various fuselage 
stations for ranges in the effectiveness factor, e. 


The value of 1^ for e = 0. 5 is 57. 8 (not including the doors); the fully effective 
value is 71.2. The ratio is 


I (e = 0. 5) q 

yy v _ 57. 8 

Iyy(Fully Effective) 71.2 


= 0. 81 


From Subsection 4.4, the ratio of mid-fuselage influence coefficients was 


a 3_3 (ModelI) _ 35.61 
5 33 (Test) 45.87 ' 

The rather good agreement tends to lend credibility to the analytical work and 
reasoning used in determining the effectiveness factor. 

A shear effectiveness factor of two-thirds was used on all shear panels in the 
mid-fuselage region. This value is recommended in References 8, 9-10 and 9-11. 
Test results (Reference 9-12) on very thin diagonal tension beams gave ratios of 
G e ff/G ranging from 0. 45 to 0. 7, with an average of 0. 6 for 12 tested panels. This 

substantiated to some extent the analytical result of two-thirds recommended by Stein 
(Reference 9-11). Although the aspect ratios of the tested panels differ from the 1/8- 
scale model, the value of two-thirds is believed to be representative of the actual 
behavior. 


Model II contains the same fuselage grid refinement as Model I. Thus, a single 

shear panel was not used to cover the entire side wall or half of the lower deck. 

2 

Phantom members having an area of 0. 001 in. were used as intermediate rod elements 
between the panels since the use of multiple panels permits variation in load at each 
corner attachment. Subdividing the shear panels allows for a more correct interaction 
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Table 6-1 Typical Fuselage Section Properties at Sta 116.0 
(No Doors, Effectiveness Factor e = 0.5) 



Method 1 

1. lyy Preserved 

2. N.A. Preserved 

3. Assume A 2 = (Aj + A 3 ) 

Method 2 

1. lyy Preserved 

2. N.A. Preserved 

3. I zz Preserved 
Not Used 

Method 3 
Summation of 
Effective Width 
X Thickness 
Not Used 

mm 

0.057 

0.061 

0.083 

H 

0.121 

0.166 

0.161 

SgS.'l 

0.064 

0.027 

0.078 


57.8 

57.8 

- 

an 

109.6 

124.6 

- 


TT-26 


with the frames. However, when this is done, rods must be placed around the panels 
for stability. Of course, MFC's could have been used in lieu of the phantom rods, but 
experience shows that using rods is more convenient. 

The fuselage was double-noded at Station 64. 0, and a spring inserted in the axial 
direction to account for the flexibility of the joint (refer to Volume HIA, Appendix A2). 
The value for the spring constant of 148, 000 lb/in. was back-calculated from static 
test data. This value is approximately equal to a joint effectiveness factor of 70 %. 

6.1.4 WING 

Remodeling of the wing followed, in principle, the procedures used on the fuse- 
lage. A sketch of the wing effective material strips is shown in Fig. 6-5. A direct 
stress effectiveness factor of e = 0. 85 was used in both the rib and spar directions for 
the upper and lower covers. This effectiveness factor is in line with the value obtained 
in Subsection 5. 10 for a panel having an initial imperfection of (¥) = 1. 0. An 
effective shear modulus of = 0. 85 G was also used in the covers. The top and 
bottom covers were thus replaced with a rod and shear panel idealization. The 
equivalent rod areas were obtained by simply adding the effective cover material to the 
actual spar and rib caps. 

A fully effective value of G was used for the spar and rib shear panels. Here, it 
was felt that the wing ballast weights (which are attached to the webs, at the midpoints, ) 
would tend to minimize any initial imperfection, and also would act to prevent any out 
of plane deformation (panel bowing) when the structure is loaded. 
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The grid system for Model n Is much coarser than for Model I, with the mesh 
being set by the location of ribs and spars. 

The wing effective width strips were extended into the carry-through structure 
as shown in Fig. 6-5. These widths were reduced where appropriate to account for 
cutouts. Note that the top cover of the carry-through does not contain an effective width 
strip along the fuselage side wall in the longitudinal direction. Here, the cutouts for the 
fuselage frames completely destroy the longitudinal continuity of the deck. 

Static test results for loading at the wing tip were discussed in Subsection 5. 11. 
There the ratio of Model I to test data was calculated to be 0. 80. This is in good agree- 
ment with the theoretical initial imperfection calculations which indicate that an effective- 
ness factor of 0. 85 whould be used for Model n. The difference between 0. 80 and 0. 85 
can be attributed to additional flexibility in the fuselage and wing carry-through not 
present in Model I, but being accounted for in Model n. 

6.1.5 FIN 

Springs were inserted between the fin and the fuselage to account for the local 
flexibility associated with the fin/fuselage connection. The values for the spring con- 
stants were based on the difference between static test and analytical values calculated 
using Model I. In Fig. 6-6, the local springs are identified as k^, kg, and k g . The 
equivalent fuselage spring constants K^p, Kgp and K^p were obtained by evaluating the 
fin fuselage interaction forces for a unit load applied at the fin ballast. This loading 
produces a local angular distortion, 0p A (see Fig. 6-7), of the fuselage interface equal 
to 5.34 X 10 6 radians. Using this value the fuselage spring constants were found to 
be K 1f = 18, 150 lb/in. , Kgp = 18, 600 lb/in. , and K^ = 3, 770 lb/in. The difference 

in magnitude between points 1 or 2 and 3 indicates the relative softness of the aft 

_6 

frame. In the test, the relative angle, Op T » was found to be 17.51 A 10 radians. 

The difference between these two angles is a measure of the local springs that were 
added. To evaluate this redundant system the springs k^, kg and kg were placed in 
series with K 1f Kgp and Kgp and the total angle, 0p T used to relate the applied 
force and the rotation and the combined spring system. This yielded one equation with 
three unknowns. The second equation was obtained by performing the same analysis, 
but for antisymmetric loading and assuming that all load was taken by the middle support. 

For the antis ymmetric case, the local lateral analytical rotation was calculated to be 

-6 -6 

17. 5 X 10 radians, while the static test gave a value of 43. 3 X 10 radians. The 
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EQUIV. FUSELAGE STIFFNESS 


A. ANTISYMMETRIC LOADING B . SYMMETRIC LOADING 

T-i 12 Fig. 6-6 Local Springs to Account for Fin Connection Flexibility 

final equation assumed the local springs and k g to be equal. Values of the various 
spring constants are shown in Table 6-2 where the third column gives values for the 
revised equivalent local fuselage springs, including the flexibility of the connection. 
Observe that a rather large change has been made to the forward connection (approxi- 
mately 80 %). 

Recommendations for changing the design of the forward attachment were made 
early in the study. However, as mentioned previously, these suggestions were not 


Table 6-2 Local Spring Constants to Account for 
Fin Connection Flexibility 
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implemented because the workload in Langley shops indicated a large slippage in 
testing schedule would result. On balance, it was decided that more productive infor- 
mation could be obtained from continued tests of the existing experimental model to 
uncover other possible sources of disagreement. 

A side investigation (not reported in Section 5) entailed analytically "cutting” 
the forward fin connections. The effect of doing this was to drop the frequency of the 
first fuselage bending mode from 53. 2 Hz to 48. 05 Hz. Analytically "cutting" the aft 
connection produced a minor change in the frequency. 

Thought was also given to modeling the connection in detail (again as a side study) 
to obtain equivalent spring constants to incorporate in the overall analysis. This idea 
was also abandoned because of joint nonlinearities. The simplest and most straight- 
forward approach was to use the test data. 

6.1.6 CARGO DOORS 

In remodeling the cargo door region, particular attention was paid to the actual 
location of the fuselage and door longerons and the details of the straps that connect 
the door to the fuselage. These details are shown in Fig. 5-37, 5-39 and 5-40. These 
changes added greatly to the flexibility of the cargo-door/fuselage interface as dis- 
cussed previously in Subsection 5. 12. 
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6.1.7 PAYLOAD 

Static test No. 26 (a Z load on the payload at Station 117. 5) revealed relative 
motion between points A and B (Fig. 6-8) on the payload and points A' and B* on the 
fuselage. The absolute and relative deflections are shown in Table 6-3 for a unit load 
at point C on the payload. Note that Model I did not account for this relative motion. 
To account for the discrepancy, it was decided to add vertical springs equal to 21^ 
at the forward support (k^ for the half structure) and kg at each of the aft supports. 
Since the payload is determinately supported, a load, P, at Station 117. 5 produces 
reactions of 0.46 P and 0.54 P at stations 78.0 and 151. 88, respectively. 

The values of the spring constants are thus: 

k i = /own qq v = 23, 150 lb/in. and k 0 = &H°- 54) - = 32, 250 lb/in. 

1 (2)(U.99 x 10 -5) 2 (0 , 84 x 10 - ) 
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Table 6-3 Comparison of Analytical and Experimental Deflections 



P = 1 lb 

Vertical Deflection in. 


Payload Deflection Relative 
To Fuselage, in. 

Point 

Test 

Model 1 


Test 

Model 1 

A 

3.03 x 10 s 

1.39 x 10 s 


0.99 x 10 s 

0 

A' 

2.04 x 10 s 

1.39 x 10 s 


0 

0 

B 

1.87 x 10 s 

0.69 x 10 s 


0.84 x 10 s 

0 

B' 

1.03 x 10* 

0.69 x 10 s 


0 

0 

C 

4.24 x 10 s 

2.71 x 10 s 


2.74 x 10* 

1.69 x 10* 
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6.2 NASTRAN ANALYSIS FOR MODEL II 
6. 2. 1 COMPONENT RESULTS 

The NASTRAN analysis of Model n followed the procedures outlined in Subsection 
3.3, except that the springs between the payload and fuselage, and the fin and fuselage, 
were added in Phase II of the analysis. The frequencies of the first three modes for 
the various components are given in Table 6-4 together with the results obtained using 
Model I. 


Table 6-4 Frequencies for Orbiter Structures — Comparison of Models I and II 



Frequency Hz 

Substructure 

1 st 

Mode 

2nd 

Mode 

3 rd 

Mode 

1 

II 

1 

II 

1 

II 

Fuselage 

• Free, Symmetric 

• Free, Antisymmetric 

62.2 

28.8 

55.2 

25.1 

129.9 

89.1 

109.9 

73.8 

191.3 

128.3 

154.3 

110.0 

Payload 

• Restrained, Symmetric 

• Restrained, Antisymmetric 

81.2 

68.6 

8i 2 < a j 
68.6 (a > 

268.5 

175.4 

268.5* a ! 

175.4< a > 

627.7 

462.8 

627.7 (a) 

462.8 (a > 

Cargo Doors 

• Free, Symmetric 

• Restrained, Antisymmetric 

4.6 

156.4 

4.3 

151.7 

10.7 

622.2 

10.9 

319.8 

17.6 

1054.6 

19.9 

384.3 

Wing 

• Restrained 

77.6 

74.5 

158.3 

148.7 

259.9 

254.5 

Fin 

• Restrained, Symmetric 

• Restrained, Antisymmetric 

264.2 

107.8 

281. 2< b | 
137.5 (b) 

841.3 

407.2 

852.2| b | 
462. 7' b ' 

1263.3 

1018.7 

1269.4 ( (b | 

1307.1< b > 


a Springs Added in Phase 2 

^Supports Changed, Springs Added in Phase 2 
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6. 2. 2 ORBITER RESULTS 

The frequencies for the total symmetric and antisymmetric Orbiter are given in 
Table 6-5 together with the results obtained using Model I. These results and a com- 
parison with test data will be discussed in Section 7. 


Table 6-5 Frequencies for Symmetric and Antisymmetric 
Modes for Models I and II 


Mode 

Frequency, Hz 
Symmetric Case 


Frequency, Hz 
Antisymmetric Case 

Model 1 

Model II 


Model 1 

Model II 

1 

53.2 

44.2 


52.9 

42.2 

2 

62.6 

54.4 


72.6 

57.0 

3 

75.2 

63.0 


85.1 

58.6 

4 

108.5 

80.2 


92.0 

71.6 

5 

133.7 

103.5 


101.5 

78.9 

6 

156.3 

115.9 


135.3 

103.5 

7 

162.2 

121.6 


- 

120.3 

8 

175.2 

139.7 


- 

146.5 

9 

216.4 

170.9 


- 

159.2 

10 

- 

185.0 


- 

196.1 
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7 - CORRELATION OF ANALYTICAL AND EXPERIMENTAL RESULTS 


7.1 STATIC TEST CORRELATION 

The static testing procedure was described in Section 4, together with some 
Limited correlation with the Model I results. A list of the various static tests is 
given in Table 4-3 and the location of the dial gages that were used to measure de- 
flections is shown in Fig. 4-7. Table 7-1 indicates the static tests used as the basis 
of correlation between test data and both Model I and Model II analytical results. De- 
flection plots of the structure, obtained from test and analysis, are shown in Fig. 

7-1 through 7-10. Individual load deflection plots are contained in Volume niA, 

Appendix A3, and Appendix G of this volume. 

Stresses in the fuselage panels at the upper longeron, lower longeron, and 

bottom deck are also plotted in Fig. 7-1 through 7-7. The stresses have been cal- 
/Mc\ 

culated using l — j — ) where a reduced value of the moment of inertia has been used. 

The values of the moment of inertia, extreme fiber distances c^ and c 2 , and the 
reciprocal of the section modulus (rpj , are plotted in Fig. 7-11. The following 
nomenclature has been used in Fig. 7-1 through 7-8: 

/Mc\ 

0 - T (D.L.) = bending stress at top longeron level ( — 1 » due to 1-gLoad 
cr T (L-L- ) = bending stress at top Longeron level ^V|c^ ^ due tQ Live Load 

/Mc\ 

a-rj, (TOT) = bending stress at top longeron level l-y J , due to 1-gand Live Load 

/Mc\ 

o-g (D.L. ) = bending stress at bottom deck level due to 1-gLoad 

/Mc\ 

o- B (L.L.) = bending stress at bottom deck level due to Live Load 

/Mc\ 

o-g (TOT) = bending stress at bottom deck level , 1-g and Live Load 

cx c (Bot. Deck) = elastic buckling stress of bottom deck 

c j c (Sidewall-Top) = elastic buckling stress of side wall for positive bending 

(compression top) 

& (Sidewall-Bot. ) - elastic buckling stress of side wall for negative bending 

L 

(compression bottom) 
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Fig. 7-1 Static Test No. 3 - Stress and Deflection Correlation 
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Fig. 7-7 Static Test No. 21 — Stress and Deflection Correlation 
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The critical stresses, <r c » are based on linear elastic values given in Tables 
5-9 and 5-10. The cross hatched regions in Fig. 7-1 through 7-8 indicate where the 

values of a (Total) =<J'(1 ) + <j (Live Load) exceed the linear elastic critical stresses. 

( Mc\ ® 

— — Istresses are perhaps 20 % larger than the actual values as indicated in 
Fig. 7-12. Here, it is assumed that the actual stress varies as 

<t = <r o ^1 - cos^ because of the presence of an initial bow. If the nonlinear dis- 

tribution is replaced with trapezoidal blocks having an equal bending moment, then 
/Mc\ 

the linear J - ) v& l ue » represented by the trapezodial distribution, will exceed the 

maximum nonlinear value by 20 %. The critical stresses are based on the assump- 
tion that the panels are simply supported; therefore, any restraint provided by the 
frames will incr ease these values. The orosshatched regions thus tend to over- 
estimate the behavior, although they do indicate a basic softening of the structure 
that is reflected in the values of the influence coefficients obtained by the method 
illustrated in Fig. 7-13. The lower portion of the load deflection curves (Fig. 7-13) 
between zero and 50 % of the maximum was used to determine the loading influence 
coefficient, iTj» anc * unloading influence coefficient, & Note that the 
load deflection curves exhibit hysteresis effects (refer to Volume IHA, Appendix A3), 



T-208 Fig. 7-12 Equivalent Trapezoidal Stress Blocks for an Assumed Nonlinear Stress Distribution 
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APPLIED LOAD, LB 


n 


and that the plots for the gages at the point of load application are more nearly 
linear. The deflection plots for Model II (cargo doors on) agree rather well with 
test data (Table 7-2), as can be seen from Fig. 7-1, 7-2, 7-7, 7-9, and 7-10. 

Figure 7-8 illustrates the wing bottom cover bending stress distribution. 

✓ 

These values were obtained by fairing curves through panel average stresses. The 
drop in the aft portion of the cover stress is due to the redistribution caused by the 
cutout provided for the interstage fitting. The wing deflections for Model II agree 
closely with test measurements as shown in Figure 7-9. 



Fig. 7-13. Method Used to Obtain Experimental Influence Coefficients 
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7.2 DYNAMIC TEST CORRELATION 

Analytical and experimental results for the first five symmetric modes are 
shown plotted in Fig. 7-14 through 7-18. Table 7-3 contains a comparison of the 
analytical and experimental frequencies for these same modes. Note that Table 7-3 
is aligned in ascending order of the Model II analytical results. Descriptions of the 
basic behavior of these analytical results were developed, then correlated with a 
test mode having the same behavior, and its frequency was entered in the table. The 
mode corresponding to the test frequency of 104.1 Hz is somewhat questionable, 
being very similar to the 106.8-Hz mode. The 104.1-Hz mode was not measured in 
the second (vertical suspension) mode survey. Note, also, that two modes were 
measured for the wing first bending mode. The 51.2-Hz frequency corresponds to 
resonance of the left wing, while the 54.2-Hz frequency corresponds to that of the 
right wing. Mode survey data is contained in Appendix C of this volume; additional 
experimental data is contained in Reference 9-5. 

Referring to Table 7-3, the analytical results for the first five symmetric 
free-free frequencies of Model n are within an error of 8.3 % of the experimental 
results, and four of these frequencies agree within an error of 3.2 % or less. Table 
7-3 also indicates the large improvement between Model I and Model n results when 
compared with test data. 

The mode shapes are presented in Fig. 7-14 through 7-18. In reviewing these 
figures it may be noted that there is fairly good agreement between analysis and test. 


Table 7-3. Comparison of Analytical and Experiment Frequencies for the Symmetrical Free-Free Normal Modes 


Mode 

Model 1 
Freq, Hz 

Model II 
Freq, Hz 

Test 

Freq, Hz 

Model 1 
to Test 
% Error 

Model II 
to Test 
% Error 

Description 

1 

53.2 

44.2 

43.6 

+22.0 

+ 1.4 

Fuselage 1st bending 

2 

62.6 

54.4 

51.2 1 

54.2 I 

+ 18.8 

+3.2 

Wing 1st bending (vs payload vertical) 

3 

75.2 

63.0 

58.2 

+39.2 

+ 8.3 

Wing 1st bending (vs aft fuselage vert) 

4 

108.5 

80.2 

80.1 

+35.4 

+ 0.1 

Fin fore-aft 

5 

— 

103.5 

104.1 <a) J 
106.8 ) 

- 

-3.1 

Payload - aft vertical 

6 

133.8 

115.9 

- 

- 

- 

Payload fore-aft (vs fwd fuselage fore-aft) 

7 

162.3 

121 5 

~ 

- 

- 

Aft fuselage pitch 

8 

133.8 

139.7 

- 

- 

- 

Wing 1st torsion (vs fwd fuselage fore-aft) 

9 

175.3 

170.9 

- 

- 

- 

Fuselage 2nd bending & wing fore-aft 

10 

216.5 

185.0 

- 

- 

- 

Wing fore-aft bending 


Questionable mode; very similar to 106.8 Hz mode, also not measured during the second survey on 1-3-74. 
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T -210 Fig. 7-14 Normalized Mode Shape - First Symmetric Mode 





Normalized Mode Shape - Second Symmetric Mode 
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Fig. 7-17. Normalized Mode Shape - Fourth Symmetric Mode 
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Table 7-4. Substructure Contributions to the Z Component of Momentum 
for the First Free- Free Mode for Model I and Model II 


Substructure 

Weight 

% 

Model 1 

Model II 

m) z 

% + 

(M 0) 2 

% + 

Fuselage 

100.91 

48.34 

-0.1219 

— 

-0.1180 

— 

Payload 

64.62 

30.96 

+0.0661 

46.6 

+0.0970 

69.5 

Wing 

32.86 

15.74 

+0.0675 

47.5 

+0.0367 

26.3 

Cargo Doors 

6.52 

3.12 

+0.0084 

5.9 

+0.0061 

4.2 

Fin 

3.85 

1.84 

-0.0201 

- 

-0.0218 

- 

2 


100.00 

0.0 

100.0 

0.0 

100.0 
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These mode shapes have been normalized to the most representative displacement 
component. Thus, the first mode has been normalized to a unit motion of the mid- 
fuselage (Station 125.5), and not to the maximum displacement which occurs at the 
fin ballast. From Fig. 7-14 observe that the first mode for Model I agreed closely 
with test data, the exception being the fin ballast motion, although the frequency was 
in error by 22.0%. The mode shape for Model II agrees rather well with test data 
for the fuselage, has better agreement for the fin, but shows less agreement for the 
wing tip. The Z components of the momentum for the five substructures are tabu- 
lated in Table 7-4 for Model I and II. Although the fuselage displaces up through the 
center portion, the total fuselage momentum is negative because almost all of the 
mass is concentrated at the ends (see Fig. 5-4), which displace down. The fin adds 
slightly to the negative momentum of the fuselage, and this is balanced by the posi- 
tive momentum of the payload, wing and cargo doors which in Model I contribute 
46.6, 47.5, and 5.9%, respectively, to the total positive momentum. In Model n 
the negative momentum of the fin and fuselage, which is approximately the same as 
Model I because of the similarity of the mode shape, is balanced by different per- 
centages for the other three substructures, predominantly the payload and wing 
(refer to Table 7-4). The moment of inertia of the fuselage for Model II is approxi- 
mately 80 % of the Inertia for Model I, while the Model II wing moment of inertia is 
approximately 85 % of the Model I value (refer to Subsections 6.1.3 and 6.1.4). 
Springs calculated from static deflection data were added to the Model II payload 
support structure (Subsection 6.1.7) which were not present in Model I. Thus, on 
a relative basis, the Model n wing is stiffer than the Model I wing, while the Model n 
payload is softer compared to the respective fuselage stiffnesses. Because the pay- 
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Table 7-5. Estimated Stresses for Static and Dynamic Tests 


L< 


Horizontal Suspension 
Dynamic Test* 

Static Test No. 3 

Static Test No. 4 

Static Test No. 13 


Ifl 

Dyn 

Max + 

Max - 

Ifl 

P z = -80 

Total 

Ifl 

p z “ +120 

Total 

Ifl 


Total 

1 

9i 

3 

LL 

Upper 
Longeron 
Sta 116 

a c = -270 

+200 

±321 

+521 

-121 

-200 

-320 

-520 

-200 

+490 

+290 

B 


■ 

Lower 
Longeron 
Sta 116 

o c --533 

o c (Bot 
Deck) 

- -167 

-100 

±193 

+ 93 

-293 

+ 100 

+150 

+250 

+ 100 

-230 

-130 


8 

1 

Wing 

Upper 
Cover 
Sta. 153, 
B.L. 21 

o c -- 192 

- 80 

± 71 

- 9 

-151 

8 


/ 

\ 


8 

+80 

+240 

+320 

Lower 
Cover 
Sta. 153, 
B.L. 21 

o c = -192 

+ 80 

± 71 

+ 151 

1 

8 

| 

8 

8 




-200 

-280 


*3 lb Harmonic Force at Nose Fitting - 1st Symmetric Mode 
TT 195 


load weight is twice the wing weight, a small change in the payload motion must be 
offset by a larger change in the wing motion. There is evidence that the payload 
springs and the wing stiffness have not been correctly "tuned" to produce the correct 
behavior. 

Table 7-5 lists estimated stresses for the static and dynamic tests for the fuselage 
longerons and wing cover. In the inverted dynamic test, the dynamic and maximum 
stresses in the bottom deck exceeded critical stresses. In static test No. 3, the 
stress in the upper portion of the side wall panel due to the applied load exceeds the 
critical stress of -270 psi. The applied load in static test No. 4 was not large enough 
to produce a net compressive stress in the bottom deck in excess of the linear elastic 
critical values. Hence, there is some question on the suitability of static test data 
for determining actual dynamic elastic characteristics and, in the case of the fuse- 
lage, some questions as to whether static tests fully covered the range of dynamic 
behavior. For the wing, results of static test No. 13 appear to have adequately 

covered the dynamic test range. All of the static applied load and 1 g stresses in 

/ Mc\ 

Table 7-5 have been estimated using ( — j — ) with a reduced value for the moment 
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of inertia. The dynamic stresses were obtained by multiplying NASTRAN strains 
by the test amplitude. The static load deflection plots do not show an abrupt change 
in stiffness that would be expected for linear buckling; they do show a slight "bending 
over" as might be expected for a bowed panel. In addition hysteresis is present in 
the loading and unloading cycle. 

Table 7-6 lists numerical values for strain and deflections for both analysis 
and test results for the first free-free symmetric mode. These values verify the 
good agreement noted previously. 


Table 7-6 Selected Analysis and Test Results 
for First Free-Free Symmetric Mode 



Model 1 

Model II 

Test 

Frequency 

53.2 

44.2 

43.6 

e x (Top Longeron) 

0.0102 

0.0112 

- 

e v (Wing Top Cover) 
Sta 153.4 (BL 21) 

-0.00302 

-0.00098 


5 z (Sta 125.5) 
Node 1105 

1.0 

1.0 

1.0 

6 X (Fin Ballast) 
Node 4400 

3.36 

3.73 

4.0 

6 z (Sta 116) 
Node 905 

1.04 

1.05 

1.0 

6 Z Wing Tip 
Sta 162 
Node 3017 

1.53 

0.40 

1.1 

5 2 Payload 
Mid-Point 
Node 4886 

0.44 

0.63 

0.5 

6 Z (Nose) 
Node 110 

-1.43 

-1.60 

-1.5 

6 Z (OMS Ballast) 
Node 2200 

-0.76 

-0.64 

-0.85 ± 

6 X (OMS Ballast) 
Node 2200 

1.18 

0.97 

1.1 ± 

6 Z (Approx. Wing CG) 
Node 3413 

0.69 

0.28 

0.55 ± 


TT-196 
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8 - CONCLUSIONS AND RECOMMENDATIONS 

Comments and recommendations concerning the analysis and correlation fall 
within three basic categories: modeling procedures, correlation, and experiences 
with NAS TRAN. 

8.1 MODELING PROCEDURES 

In NASTRAN Model I, the CQDMEM2 element was used to model a fully ef- 
fective structure. Lack of correlation should not be attributed to the behavior of this 
element or to NASTRAN. The major cause for lack of correlation was traced to the 
ineffectiveness of the skin panels due to the presence of initial bows. Model II used 
bars and shear panels to represent the behavior of the physical model. Here, ef- 
. fective bar areas and an effective shear modulus were used. It is felt that this type 
of modeling is more descriptive of the actual characteristics of the structure. The 
effectiveness factors were obtained through the use of a computer program that 
solved the large deflection problem of in-plane loading of plates with initial imperfec- 
tions. Little information exists in published form that contain the parameter range 
and type of loading of interest of the present study. Work was done by NASA some 
20 years ago. It is recommended that this work be updated to cover a broad range 
of aspect ratios, initial imperfections and loadings. Publication of this data in chart 
form would be beneficial to modeling efforts on many structural projects. 

Some remodeling of major joints could have been made in the analysis. How- 
ever, for some joints such as the fin/fuselage connection, only static tests could give 
the exact behavior - no matter how careful an analysis was made. 

Consistent mass concepts tend to lose their significance for structures that do 
not behave in a linear, fully effective fashion. In fact, for structures of this type, 
it becomes desirable to control the mass and stiffness properties independently. 

8.2 CORRELATION 

• The results of the vibration tests reported here are limited in scope, and 
have not been fully completed 

• The static test data have a ±10% error in reproducibility 

• It is felt that the agreement of Model n results with test data is sufficiently 
close to pursue a coupling analysis. 
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8.3 COMMENTS ON AND EXPERIENCES WITH NAS TRAN 

The analysis of the 1/8-scale model was viewed as a pilot study of the use of 
NASTRAN on realistic aerospace projects. Therefore, some of the comments extend 
past specific application to the 1/8-scale model. Many areas must be considered in 
attempting to judge the NASTRAN system, some of which are not necessarily concerned 
with NASTRAN itself. 

One of the prime reasons for Grumman's interest in NASTRAN is the work that 
NASA is doing in attempting to establish NASTRAN as the industry standard. For 
multicorporation aerospace projects, standardization is a necessity. Standardization 
of computer programs, type of data, etc, is important from a contractual point of view, 
where subcontractors interface with prime contractors and the prime contractor inter- 
faces with the principal agency. However, with standardization there must be sufficient 
technical flexibility - in this regard NASTRAN has room for improvement. 

8. 3. 1 FLEXIBILITY OF THE NASTRAN SYSTEM 

NASTRAN contains a number of rigid formats that allow for specific types of 
problems to be solved. Where these formats do not fit the specific application, the 
rigid format may be changed with ALTERS. In analyzing the 1/8-Scale Orbiter Model, 
numerous ALTERS to Rigid Format 3 were made. The system was modified to handle 
substructuring using level 15.5. It was found that, although the basic system can be 
modified, it requires an extensive learning period to become sufficiently proficient 
with the system to make use of ALTERS. After making a series of extensive ALTERS, 
the modified program must then be checked; if altered again, it must be rechecked. 
Because this procedure is cumbersome, it tends to produce a small group of NASTRAN 
experts. Our own in-house system (COMAP ASTRAL; Reference 9-13) contains no 
rigid formats; only one structural command generates all the required matrices (or 
those that are requested) to solve the problem. All matrix operations are coded in a 
simple, interpretive language that may be learned in a few hours. This system re- 
quires that the engineer understand the physics of the problem. Altering NASTRAN, 
on the other hand, requires a good understanding of the systems aspects of NASTRAN. 
Thus, there is diversion of engineering talent from engineering problem-solving to 
overcoming NASTRAN system complexities. 
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8.3.2 CAPABILITY OF NASTRAN 

The original analysis plan for the 1/8-scale model required that all components 
(Orbiter, ET and two SRB's) be coupled to determine mated vehicle modes. This 
combined hydro-visco-elastic analysis is theoretically possible within the present 
NASTRAN system, but from a practical stand point it is not. Although modes were 
successfully calculated for the individual components, it was not possible to couple 
these components and analyze the total (mated) vehicle due to extremely large computer 
time requirements. Note that the lack of correlation of the Orbiter analytical and test 
modes for Model I is not a fault of NASTRAN. This lack of agreement did, however, 
cause a shift of emphasis from studying the coupling problem to examining the Orbiter 
modeling, design, and fabrication in more detail. Although the Orbiter correlation 
has been cleared up by the use of Model n, coupling is still a major problem. Two 
approaches seem feasible: 

• Model coupling in lieu of static coupling to reduce the size of the final problem 

• Incorporation into NASTRAN of approximate reduction schemes that employ 
the automatic tri-diagonal reduction algorithm. (For example, FEER (Ref- 
erence 9-3; or ALARM, being developed by Grumman under Master Agreement 
NAS 1-10635, Task 17. FEER and ALARM would have to be extended to 
handle complex modes). 

Master Agreement NAS 1-10635, Task 21, will pursue the model coupling approach 
by making the appropriate ALTERS to the NASTRAN rigid formats. 

Many complaints of NASTRAN are associated with the form of the stress output. 
Average stresses in elements are of little use to the engineer in the design mode. 
Traditional cap loads and shear flows are preferable. At this stage the designer is 
looking for the best load paths, and will rearrange the framing to the best of his ability 
to obtain it. The elements of the structure are proportioned to withstand the load 
imposed on them using allowable stresses, therefore the average stress in an element 
is not a convenient quantity. NASTRAN level 15. 5 produces element corner forces 
which may be used as input to a post processor to produce equivalent cap loads and 
shear flows. The Grumman post processor produces listings of member loads not only 
by condition, but by critical condition ranking as well. 
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Tabl« 8-1 1/8-ScaU Model Obiter NASTRAN Analysis Computing Statistics 


Anal. 

Description 


1 . Wing - Phase 1 


2. Cargo Doors - Phase 1 (Symmetric Case) 


( Antisymmetric Case) 


3. Payload Beam - Phase 1 (Symmetric Case) 

_ 

( Antisymmetric Case) 

*0 

4. Fin — Phase 1 (Symmetric Case) 

o 

( Antisymmetric Case) 


5. Fuselage — Phase 1 (Symmetric Case) 

c 

( Antisymmetric Casel 


6. Orbiter - Phase 2 (Symmetric Case 

w 

( Antisymmetric Case) 

o 



(Took 2 % mo. calendar time starting from uncor- 


ractad data) 


Good Runs 



1 . Wing - Phase I (Not Run, Same as 

2. Cargo Doors — Phase I Original Model) 

3. Payload Beam — Phase I (Symmetric Case) 

( Antisymmetric Case) 


4. Fin - Phase I 

5. Fuselage - Phase I 

6. Orbiter - Phase 2 


(Symmetric Case) 

( Antisymmetric Case) 
(Symmetric Case) 

( Antisymmetric Case) 
(Symmetric Case) 

( Antisymmetric Case) 



Total Time Used (1 week) 


6. Symmetric Orbiter — Phase 2 (Fwd Fin Sup- 
port Cut) 

(Aft Fin Sup- 
port Cut) 

Total Time Used (1 week) 


1. Wing — Phase I (Skin 50% Eff) 

(Skin 85% Eff) 

2. Cargo Doors - Phase I (Symmetric Case) 

( Antisymmetric Case) 

3. Payload Beam — Phase I (Symmetric Case) 

( Antisymmetric Case) 

4. Fin - Phase I (Symmetric Case) 

( Antisymmetric Case) 

5. Fuselage Phase I — G = 50% Eff (Symmetric 

Case) 

( Antisymmetric Case) 
Fuselage Phase l-G = 2/3 Eff (Symmetric 
Case) 

( Antisymmetric Case) 

6. Orbiter — Phase 2 (Symmetric Case) 

( Antisymmetric Case) 


Total Time Used 

(Took 3 mo. calendar time which includes 
extensive data revision) 


Total Time Used For Orbiter Mode Analyses - 


CPU | Sys No. of I CPU Sys No. of 

Min Run* I Min Min Runs 
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8.4.3 COMPARISON OF COMPUTER RUNNING TIMES 

A complete record of the amount of computer time needed to debug and analyze 
the various Orbiter models using NASTRAN is presented in Table 8-1. Table 8-1 
shows the NASTRAN computing statistics for obtaining the Orbiter symmetric and 
antis ymme trie dynamic modes by the substructuring method of analysis. An orginal 
NASTRAN analysis of the l/8-scale Orbiter (prior to minor structural changes), and 
two revised analyses known as Model I and Model n, are shown. In addition, two 
Phase 2 analyses of Model I were made. The NASTRAN Rigid Format 3 which included 
extensive DMAP ALTERS for checking purposes was used throughout. For each item 
in Table 8-1, undeformed and deformed plots were also obtained. The original analysis, 
as expected, used up most of the debugging time. The model I analysis profited con- 
siderably from the availability of the corrected input data generated in the original 
analysis, as can be seen in Table 8-1, since only minimal changes were made. Model 
II had extensive data revisions and thus a somewhat larger debugging time than Model I. 
In all of the foregoing three analyses however, production run times were about the 
same. This is because there was no major change in idealization size or number of 
dynamic DOF 's in the three analyses . 

A summary of the total amount of computing time used on the NASTRAN analysis 
is, therefore, 411 CPU (1522 System) minutes for 88 debugging runs, and 524 CPU (2736 
System) minutes for 40-good runs. In addition time to develop ALTERS for Rigid Format 
7 (for the ET and SRB analyses took 40CPU (310 System) minutes. 
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The flange mass has been included as non-structural mass of .0026 #/in. 

3) Rods 351 9 -* 3580, 3583 -* 3587 are vertical web bars. They represent 
actual web connection angles, or if none, represent fictitious members 
to transfer load into shear panels. No mass assign to these members, 
since ballast was assumed to cover whole webs. 

4) Rods 3581 and 3582 represent Interstage connection Lug. 
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1/2 BALLAST WEIGHT CONCENTRATED AT 4400 
W = 3.0 LBS 

l x (ROLL) “ 1.563 LB IN. J 
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END RIB (ALL ALUM, p = .1 #/IN 3 ) 


SHEAR WEBS t = .04 
RODS A - .034 IN 1 
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NOTES ot m w 

1) ALL SPAR WEBS ARE SHEAR ONLY ELEMENTS (t = .032 A8, p - .1 #/IN 3 

2) OUTER SPAR CAP AREAS ARE EQUIVALENT AREAS TO GIVE CORRECT 
MOMENT OF INERTIA OF ACTUAL CHANNEL. NO DENSITY p ASSIGNED 
TO THESE RODS ONLY NON-STRUCTUR AL WEIGHT TO INCLUDE FLANGES 
AND OTHER NEGLECTED MATERIAL. 

3) RODS 4472 -* 4486 REPRESENT FICTITIOUS RODS TO TRANSFER LOAD 
INTO SHEAR PANELS (A = .034 AL, NO p) 


4) 


5) 

6 ) 


< 

H 

</> 

CENTERLINE SPAR CAPS ARE EQUIVALENT AREAS SO THAT WHEN 
ADDED TO OUTER CAPS, GIVE THE CORRECT CROSS SECTION AREA OF 
CHANNEL. NO WEIGHT ASSIGNED TO THESE RODS. 

COVER MPC'S RELATE DEFLECTIONS NORMAL TO COVER TO ADJACENT 
SPAR DEFLECTIONS. 

FOR SYM. CASE ONLY, THE LOCAL D.O.F. NORMAL TO PLANE OF SPAR 
WILL BE MPC'ed TO OUTER SPAR NODE D.O.F. 
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NOTES 

1) ORBITER PAYLOAD BEAM REPRESENTED BY BEAMS 4882 -» 4889 (ALUMINUM) 

SECTION PROPERTIES ARE FOR X STRUCTURE (A = 5.775 IN’, l yy - 23.23 IN 4 , l„ -= 21.87 IN 4 . J « 31.0 IN 4 ) 

2) STRUCTURAL DENSITY OF P - .1 #/IN* WAS USED. THE 54 BALLAST WEIGHT OF 18.1 LBS WAS 
DISTRIBUTED AS NON STRUCTURAL WT OF .245 #/IN. 


3) 

4) 


ALL SUPPORT BRACKETS, ASSUMED RIGID. THEIR WEIGHTS WERE DISTRIBUTED AS CONCENTRATED 
WEIGHTS AS SHOWN IN SKETCH 

THE TORSIONAL WEIGHT INERTIA OF PAYLOAD BEAM & BALLAST WAS INCLUDED AS CONCENTRATED MASS. 
THIS MASS. AT EACH NODE BASED ON 4.3 #/IN* 

iFT 


to 


CONCENTRATED WT 
BEAM OVERHANG - .55# 
BRACKET - .19 
.74# 

z 

Y^U 



10 

8 

< 

h 

CO 

<t BEAM (WL 62.5) 





CONCENTRATED WT 
BRACKET - .7# 



CONCENTRATED WT 
BEAM OVERHANG = 55# 

BRACKET = .92 


1.47# 


RIGID BRACKET 


C.G. ON CAT WL 51.24 
WT OF 1.31# BEAMED TO 4890 & 4892 


NODE 

TORSIONAL MASS 
(# IN*) 

4882 

25.13 

4883 

40.85 

4884 

40.85 

4885 

44.08 

4886 

40.85 

4887 

37.63 

4888 

38.16 

4889 

36.28 

4890 

23.24 


INTERFACE FORCES 

RIGID BRACKET 

WT OF 1.4# ASSUMED ON 8 PLANE 
C.G. LIES APPROX. MIDWAY BETWEEN 
NODES 4889 & 4892. THEREFORE BEAM 
WEIGHT TO THESE NODES 


STA 151.875 
WL 56.7 
SsTBL 10.125 

PAD WEIGHT OF .17# 

^ CONCENTRATED AT 4891 

INTERFACE FORCE 
PAD WEIGHT OF .15# 

CONCENTRATED AT 4892 

CONCENTRATED WT AT 4.12 
PAD -.16 

BRACKETS - .39 & .7 - 1.04 
1.24 


AFT MPC'S (RELATE RIGID BODY MOTIONS) 


, — BRACKET C.G. AT WL 66.7 
' j BRACKET WEIGHT OF .43# BEAMED TO 4881 & 4882 
/ RIGID CONNECTION 


ASSUMED BETWEEN 
GRID PTS. 4881 & 4882 


MPC'S AT STA 78 (RELATE RIGID BODY MOTION) 


WL 51 .933 

CONCENTRATED WT ' 1,4 

FOR BRACKET - .24# 


INTERFACE 

FORCES 


CASE 


SYM 


SYM 


ANTI 


DEPENDENT 

D.O.F. 


A 4881 X 


A 4882 Z 


B 4882 X 


INDEPENDENT D.O.F. 


A 4881 Y 


.09463 


A 4881 Z 


1.0 


A 4882 X 


1.0 


A 4882 Y 


-.09463 


0 4882 Y 


-10.667 


CASE 


SYM 


SYM 


ANTI 


SYM 


ANTI 


SYM 


ANTI 


SYM 


SYM 


ANTI 


SYM 


ANTI 


SYM 


DEPENDENT 

D.O.F. 


A 4889 X 


A 4889 Z 


0 4889 X 


0 4889 Y 


A 4890 Y 


A 4890 Z 


0 4890 X 


0 4890 Y 


A 4891 X 


A 4891 X 


A 4891 Y 


A 4891 Y 


A 4892 Z 


INDEPENDENT D.O.F. 


A 4890 X 


1.0 


■ 78409 


■09091 


.09091 


.47273 


1.0 


04890 Z 


10.126 


A 4891 Z 


1.0 


-.09877 


1.08642 


1.0 


.09877 


.51358 


1.0 


A 4892 X 


-.78409 


.09091 


-.09091 


.52727 


A 4892 Y 


1.0 


.52727 


1.0 


ORBITER LONG. PAYLOAD 


A 4890 Y 
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Appendix B 

ORBITER STATIC TEST CASES - ANALYTICAL DEFLECTIONS FOR 
MODEL I SUPPORTED ON INTERSTAGES 


TABLE OF CONTENTS 


NOTE 

Appendix B consists of the 
following illustrations: 

Figure Page 

B-l Fuselage Load - Sta 46.75, Grid 120 B-2 

B-2 Fuselage Load - Sta 116.0, Grid 906 B-3 

B-3 Wing Tip Load - Grid 3017 B-4 

B-4 Fin Ballast Load - Grid 4400 B-5 

B-5 OMS Ballast Load - Grid 2200 B-6 

B-6 Payload Load - Sta 117.5, Grid 4885 B-7 

B-7 Payload Moment - Sta 117.5, Grid 4886 B-8 

B-8 Orbiter Fuselage - Dead Load Down B-9 


B-l 



B-2 


T-120 Fig. B-1 Fuselage Load • Sta 46.75, Grid 120 
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T-124 Fig. B-5 OMS Ballast Load - Grid 2200 
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T-127 Fig. B-8 Orbiter Fuselage ■ Dead Load Down 


Appendix C 

MODE SURVEY DATA 



TABLE OF CONTENTS 


NOTE 

Appendix C consists of 
the following illustrations: 

Figure Page 

C-l Preliminary Mode Survey Data (7-12-73), f = 43.6 Hz 

Symmetric Z- Excitation C-2 

C-2 Preliminary Mode Survey Data (7-13-73), f = 51.2 Hz 

(LH Wing Resonance) Symmetric Z -Excitation C-3 

C-3 Preliminary Mode Survey Data (7-13-73), f = 54.2 Hz 

(RH Wing Resonance) Symmetric Z -Excitation C-4 

C-4 Preliminary Mode Survey Data (7-12-73), f = 58.2 Hz 

Symmetric Z- Excitation C-5 

C-5 Preliminary Mode Survey Data (7-12-73), f = 80. 1 Hz 

Symmetric Z-Excitation C-6 

C-6 Preliminary Mode Survey Data (7-12-73), f = 95.7 Hz 

Symmetric Z-Excitation C-7 

C-7 Preliminary Mode Survey Data (7-12-73), f = 104. 1 Hz 

Symmetric Z- Excitation C-8 

C-8 Preliminary Mode Survey Data (7-12-73), f = 107.0 Hz 

Symmetric Z-Excitation C-9 

C-9 Preliminary Vertical Suspension Frequency Sweep, Effect of 

Preload C-ll 

C-10 Preliminary Mode Survey Data (1-3-74) 0-lb Static Load, f= 

43.67 Hz Symmetric Z-Excitation C-13 

C-ll Preliminary Mode Survey Data (1-3-74), 0-lb Static Load, f= 

80.06 Hz Symmetric Z-Excitation C-14 

C-12 Preliminary Mode Survey Data (1-3-74), 600-lb Static Load, f= 

44. 23 Hz Symmetric Z-Excitation C-15 

C-13 Preliminary Mode Survey Data (1-3-74), 600 lb Static Load, 

f = 80. 0 Hz Symmetric Z-Excitation C-16 

C-14 Preliminary Mode Survey Data (1-3-74), 600-lb Static Load, 

f = 96.85 Hz Z-Excitation C-17 

C-15 Preliminary Mode Survey Data (1-3-74), 600-lb Static Load, 

f = 106. 84 Hz Symmetric Z-Excitation C-18 

C-16 Preliminary Mode Survey Data (1-8-74), 0-lb Static Load, f= 

43.6 Hz Antisymmetric Y- Excitation C-19 

C-l 



C-2 












C-3 





C-4 









C-5 








C-6 







C-7 


Fig. C-6 Preliminary Mode Survey Data (7-12-73), f * 96.7 Hz 
Symmetric Z-Excitation 
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Fig. C-9 Preliminary Vertical Suspension Frequency 
Sweep — Effect of Preload 
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Fig. C-12 Preliminary Mode Survey Data (1-3-74), 600-lb Static Load, f-44.23 Hz 
T -1 39 Symmetric Z-Excitation 
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ERROR BOUNDS FOR APPROXIMATE FREQUENCIES OF 
FREE-FREE UNDAMPED VIBRATION SYSTEMS 
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Appendix D 


ERROR BOUNDS FOR APPROXIMATE FREQUENCIES OF 
FREE-FREE UNDAMPED VIBRATION SYSTEMS 

In very large dynamic systems, the dynamic degrees of freedom (DOF's) have 
to be reduced through a reduction scheme. This is necessary if the problem is to be 
solved economically and within a reasonable time. Therefore, It is important to have 
a check built into the program that would indicate the deviation of the approximate 
solution of the reduced system from the exact solution of the unreduced system. 

A method for error bounds was suggested by I. U. Ojalvo*; Isaacson and 
Keller use a similar method**. This method has been extended below so that error 
bounds for free-free systems can be obtained. 

The exact eigenvalue problem is given by 

[K] = X. [M] {*.} (D-l) 

where X. is the exact eigenvalue for the i th exact mode shape {*.} . 

For free-free systems [k] will be singular. To avoid inversion problems 
later on, add to each side of Equation D-l «[mJ {*.}, where € is a non zero in- 
crement added to the eigenvalue X. 


( [K] + € [M] ) =(X t + €) [M] 1^}. 

(D-2a) 

or 


[K m ] =<*,+«> [M] t*-,l 

(D-2b) 

where 


[k m ] = [k] + «[M], 

(D*3) 

Solving Equation D-2b, 


I*,! -<*, + «> [k m ] - 1 [M] )*,| 

(D-4a) 


♦Ojalvo, I. U. , "An Error Bound For Approximate Modes and Frequencies", 
Grumman Memorandum No. STMECH 73. 66 (November 1973) 

♦♦Isaacson, E. and Keller, H. B., Analysis of Numerical Methods , John Wiley & 
Sons, New York (1966, page 141) 


D-l 


or 

(rr^M - w 1 M I'll D - 4b 

for more than one mode Equation D-4b may be written as : 

[*] [a] -1 = [k m ] _1 [m] O] (°- 5 ) 

where [A] is the diagonal matrix given by 


M- 


v « 


and W is the modal matrix whose columns are the vectors 

Since it is always possible to expand an approximate mode, {♦j} . into a linear 
combination of arbitrarily normalized exact mode shapes, let 


I*, I - M )b,( 


(D-6) 


The element, 4>^., for the i^ approximate mode would be stated as 

♦, . = 2, . b.. 

ki j kj ji 

where b.. is the j th element of the vector of constants {b.}. 

The arbitrary normalization on M is taken as 


W T [M] [♦] = [I] 


(D-7) 


where [I] is a unit matrix. 

Substituting the approximate mode shape and eigenvalue into the modified t 
dynamic equilibrium equations will yield residual forces. For the i^ mode, the 
residual forces are 


= C k m ]|*iI -<x, + «) M )*i( 


(D-8) 


.th 


where A . is the approximate i eigenvalue ,and | Jthe corresponding eigenvector. 


D-2 


Deflections caused by the residual forces are obtained from the equation 


M = [K m ] - 1 {F, ( • 

Substituting Equation D-8 into D-9, we obtain 

f“,( =)♦[! -(\ + «) [Km ]' 1 M 

Substituting Equation D-6 into D-10 results in 

{»,[= W|k,|-lV') [k m ] ' X [m] [♦] fb t }. 


Using Equation D-5 reduces D-ll to 


Kh [♦] [CO -<V«> 


Let 


[E] = [I] -(X i + <) [A] 


-1 


\ 


X i + < 

1 - 


V‘. 


and substitute Equation D-13 into D-12 

{u,h [*] [K] ib t } . 

The modal mass associated with { u. \ , using Equation D-14, is 
)u,[ T [M] {u 1 )Hb,[ T [E] T M T [M] [♦] [E]{b,[. 


Making use of Equation D-7 and D-13, Equation D-15 reduces to 


{u 1 ( T [M] {a t (= ib t J T 

/ 



Where j = 1, 2 , , but the i index always refers to the same approximate 
eigenvalue in Equation D-16. 


(D-9) 


(D-10) 


(D-ll) 


(D-12) 


(D-13) 


(D-14) 


(D-15) 


(D-16) 


D-3 


Equation D-16 can be stated as 


l“,l T [«] 


" ("W. 


(D-17) 


Since the terms in the summation are squared and yield positive quantities, it 


is possible to select the smallest value 


- m 


, which usually occurs 


when j = i (exceptions occur when more than one X ^ closely surround X j) 




Therefore, Equation D-17 can be stated as 


Kf T M <u il- T~rrj s ) <b ) 


/• 


The approximate modal mass is [m] 

Making use of Equations D-6 and D-7, 

)*il T [M] [*] T [M] [♦] )b,( 


(D-18) 


(D-19) 


(D-20) 


D-4 



Substitute Equation D-20 into D-19 to obtain the desired ERROR BOUND equation 


(M T [M] {a,| 

l5,l T [m] {*,) 



V‘\ 
V* / 


(D-21) 


Change Equation D-21 into matrix form so that the error bounds for more than 
one mode can be obtained simultaneously by matrix operations 


([+] T [M] [+])" 1 [U] T [M] [U]> 



(D-22) 


where 

[♦ ] = Approximate mode shapes for any number of modes desired. They are 
obtained from an analysis using any reduction scheme. Do not include 
rigid body modes, as they will cause singularities 

[u] = [K m ] 1 [f] = Deflections caused by residual forces. K M is defined 
in Equation D-3 

[F]= [K ] [#] - [m] [«*>] [A] = Residual forces 





+ € 

\ 

\ 


= Diagonal matrix of approximate eigenvalues, incre- 
mented by € , obtained from analysis. The rows and 
columns must correspond to the columns of 


Let 


[c] =([*] T [m] [♦Jr 1 [u] T [M] [u] 

substitute Equation D-23 into D-22 


(D-23) 


[C] > 


V 

V 


(D-24) 


The resulting [c ] matrix will not be a diagonal matrix. The off-diagonal or 
cross mode terms are discarded since they were not used in deriving the error 
bound Equation D-21. Therefore, the i*"* 1 error bound equation from Equation D-24 
would be 




V* j 


where C.. is the i^ diagonal term corresponding to the i^ mode. 
Let 


€ 


= a 


i . 


Using Equation D-26 in the reduction of D-25 yields 


(D-25) 


(D-26) 


^i 

’r’V« 1 yc7i>. 


(D-27) 


Keep in mind that in computing the error bound for each mode, . changes and 

can be obtained from Equation D-26. In most cases, « can be taken small enough 

so that a. can be set to zero. 

1 


D-6 


Changing the eigenvalues to frequencies in Equation D-27 gives 



2 


>(W5L) 

( 1+ °i v'S) 


or 



which is the desired frequency error bound for the i^ mode. 


(D-28) 


D-7 
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LARGE DEFLECTION BEHAVIOR OF INITIALLY IMPERFECT RECTANGULAR 
PLATES UNDER UNIFORM END-SHORTENING AND EXTENSION 

by 

H. Armen, Jr. and E. Saleme 

Figure E-l illustrates the problem under consideration; the edges of a simply 
supported rectangular plate are loaded so as to cause a uniform shortening or ex- 
tension. The lateral motion of the plate in the x direction is prevented by rigid 
supports. 



T-i44 Fig. E-1 Geometry and Loading of Rectangular Plate 


E-l 


A gross measure of the stiffness of the above panel is represented by the load 
p(x) associated with a prescribed uniform end -shortening ft . For the case of small 
deformations with no initial imperfections, a uniformly distributed load p(x) will re- 
sult in a uniform end-shortening, ft . The load-deformation relation can be written 
as 

P(lbs)=Kft (E-l) 

where 

f b 

P=_ b J p(x) dx 

ft = ae Q , e Q is the end-shortening strain 
and 

K= _2Et_^ £ 

(1 - * ) 

The foregoing expression cannot be expected to remain valid for thin-skin 
panels having relatively large (of the order of the thickness) initial imperfections. 
The stiffness influence coefficient, K, for this case will vary not only with the as- 
pect ratio, load distribution (a uniformly distributed load will not necessarily result 
in a uniform end-shortening), and boundary condition, but also with the applied load 
level and the amount and shape of the initial imperfection. The method used to de- 
termine this coefficient for the large deflection, initial imperfection case, is based 
on a procedure presented by Timoshenko and Gere,* where assumed displacement 
functions are used to determine the total potential energy of plate. The energy func- 
tion is then minimized with respect to the set of generalized displacement coeffi- 
cients resulting in a set of simultaneous nonlinear algebraic equations. The solu- 
tion of these equations for the generalized displacement coefficients is obtained and 
used to determine the load corresponding to a prescribed end-shortening state. 


♦Timoshenko, S. P. and Gere, J. M. , Theory of Elastic Stability , McGraw-Hill, 
pp. 411-418, 1961 


E-2 


El METHOD OF ANALYSIS 

Using the coordinate system shown in Fig. E-l, an assumed expression for 
the initial imperfection is given by the following expression that satisfies the bound- 
ary conditions 

w o = w 0 cos ar cos Hr (E - 2 > 

where W q is the magnitude of the initial imperfection at the center of the panel. An 
illustration of this assumed shape for Section A-A is shown in Fig. E-2. A similar 
shape exists along the y-axis. 



z 


T-145 Fig. E-2 Initial Imperfection Along Section A-A 


The application of uniform end-shortening along the edges y = ±a results in a 
deformation state that can be described by the displacement components in the x, y, 
and z-directions as u, v, and w^, respectively. As suggested by Timoshenko and 
Gere, these components of displacement are written in the following form: 

u = Uj sin ^-cos^L (E-3a) 

v = Vj sin cos ey (E-3b) 

and 

W 1 = W 1 cos cos la” (E-3c) 


E-3 


where U^, Vj, and are generalized displacement coefficients to be determined. 
The form of the above equations for the in-plane displacements u and v were chosen 
to satisfy the boundary conditions; u = 0 along x = ±b and y = ±a, v = 5 (due to ey) 
along y = ±a. The total lateral displacement may be written as 

w=w 0 +w i = (W + W i> 008 -W cos (E_4) 

where, from Equation E-3c, it is assumed that the lateral displacement due to an 
end-shortening is of the same functional form (shape) as the initial imperfection. 
This assumption is reasonable for a uniform end-shortening condition only. 


The nonlinear strain displacement relations may be written in terms of the 


bending and stretching components as follows: 


b s 

e = e + e = 

XXX 


,2 
d W „ 


dx 


+ || +1 /2 


'<?w\ 2 ( dw o\ 2 

JFJ -(ar-/ J 


(E-5a) 



(E-5b) 





,2 

W 1 <?u «?v 

— + — + — + 

dx dy dy dx 


dw dw o 

dx dy dx 



(E-5c) 


The first term e* 3 in the foregoing equations represents the bending components 

g 

of the strain, and the remainder e is associated with the middle surface strains. 

The total strain energy (stretching and bending of the sheet) is given by 
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For a prescribed uniform end-shortening condition the constants U^, V^, and 
Wj are foimd by minimizing the total strain energy. The conditions associated with 
this minimization are 
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The three conditions represented by Equations E-8a through E-8c may be 


written in the following form 
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Equations E-9a and E-9b can be used to obtain expressions for and in terms 
of W 2 . Substituting these expressions into Equation E-9c, we obtain the following 
cubic equation in terms of 

wj + pW 1 - qW Q = 0 (E-10) 

where 
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If we define Q = (p/3) + (q/2) , then the real root of Equation E-10, for 


Q > 0, is 


w i ~y y ~2~ 

For the case where Q £ 0 then 
where 
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A solution for and subsequent determination of and provides suffi- 
cient information to determine the strain and corresponding stress state of the 
middle plane of the plate. If we define a = vx/2b and 0 = 7ry/2a then we can 
write the strains. Equations E-5a through E-5c, in terms of the displacements 
Equations E-3a through E-3c as 
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The load along the edges y = ±a can be determined from 
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Fig. F-2 Effective Fuselage Width Calculations at Sta 64 


F-3 




























A 

z 

A x z 

z - z 

A(z-z)» 

•ooy 

y 

Ay 2 

•ooz 


0 

0.100 

14.6 

1.46 

1.96 

9.92 

— 

12.5 

mm 



0 

0.072 

0 

0 

-4.64 

1.55 

— 

12.5 

Bo 



0 

0.069 

0 

0 

-4.64 

1.49 

- 

0 

0 



© 

0.31 2e 

0 

0 

-4.64 

6.72e 

— 

6.25 

12.19e 

4.06 1 1 -( 1 -e) 3 ) 


0 

0.292e 

7.3 

2.13 

2.66 

2.07e 

5.1 8[ 1 -( 1 -e) 3 ] 

12.5 

45.6e 



Z 

0.241 


1.46 


13.0 



26.9 




0.604e 


+2.13e 

1 

+8.8e 



+57. 8e 



e 

2A 

7 

SAxz 


SAlz-T) 1 

^Oooy 


SAy 2 

2 Iqoz 

•zz 

0.3 

0.422 

4.98 

2.10 


15.3 

3.4 

18.7 

44.2 

2.7 

46.9 

0.5 

0.543 


2.53 


17.4 

4.5 

21.9 

55.7 

3.6 

59.3 

0.7 

0.663 

4.45 

2.95 


18.4 

5.0 

23.4 

67.4 

4.0 

71.4 

1.0 

0.845 

4.25 

3.59 


21.8 

5.2 

28.0 

84.7 

4.1 

88.8 


Determine equivalent skin cap areas with Neutral Axis assumed unchanged from value circled above 
Let a 2 = a! + a 3 

1) (a, + 2a 3 + 0.141) 4.64 = (a, + 0.1) 9.96 

2) (a, + 2a 3 + 0.141) (4.64) 2 + (a, + 0.1) (9.96) 2 = l yy 

(a, +0.1) 9.96 x 14.6 = l yy 

'w 

31 " 145.3 _0 - 1 

/9.96 \ 

2a 3 =(a, +0.1) ^4^4 -l) -0.041 

a 3 = 0.575 (a, +0.1) -0.021 


e 

ai 

a 3 

a 2 

SA 1 

•zz‘ 

<|< 

'zz' 

'zT 

2A 

.845 

lyy 

~w 

'zz 

88.8 

0.3 

0.029 

0.053 

0.082 

0.405 

44.2 

0.96 

0.94 

0.50 

0.67 

0.53 

0.5 

0.051 

0.066 

0.117 

0.475 

53.1 

0.87 

0.89 

0.64 

0.78 

0.67 

0.7 

0.061 

0.072 

0.133 

0.507 

57.2 

0.77 

0.80 

0.78 

0.84 

0.80 


2A‘ = (a, +a 2 +a 3 + .241) 

'zz' = ( 1 2.5) 2 (a, + a 2 + .172) 


T-147 


Fig. F-3 Effective Fuselage Width Calculations at Sta 78 
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Determine equivalent skin cap areas with Neutral Axis assumed unchanged from value circled above 
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Fig. F-4 Effective Fuselage Width Calculations at Sta 116 
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Determine equivalent skin cap areas with Neutral Axis assumed unchanged from value circled above 


Let a 2 = a t + a 3 
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Fig. F-5 Effective Fuselage Width Calculations at Sta 125.5 
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Determine Equivalent skin cap areas with Neutral Axis assumed unchanged from value circled above 


Let a 4 = a, & a 2 = a 3 
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Fig. F-6 Effective Fuselage Width Calculations at 141.75 
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NOTE: Analysis will be set up for e = .5. For Re-analysis with a different 
e, multiply the Modulus of Elasticity (E) for the equivalent caps 
which are on 4 property cards, by the appropriate area ratio from 
above curve. 


Fig. F-8 Effective Width Modification Factors 
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Fig. F-13 Effective Skin Width Ratio for Wing 
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Fig. G-2 Static Load-Defl. Curves — Run 18 Fin Ballast Applied Load (+ Px) 
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Fig. G-7 Static Fin Rotation Relative to Fuselage Run 31 — Fin Ballast Applied Load (-P) 
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Fig. 6-8 Static Torsion-Rotation Curves Run 30 — Applied Nose Torsion 
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Fig. G-9 Static Torsion —'Rotation Curves Run 30 — Applied Nose Torsion 
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Fig. G-10 Static Torsion-Rotation Curves Run 30 — Applied Nose Torsion 







i. G-11 Comparison of Test Run 30 With Analyses Nose Torsion Case 
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Fig. G-12 Static Load-Defl. Curves — Run 26 Payload Applied Load at Sta 117.5 (— P z ) 
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Fig. G-13 Static Load-Defl. Curves — Run 26 Payload Applied Load at Sta 117.5 (— P z ) 
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Fig. G-14 Static Load-Defl. Curves — Run 26 Payload Applied Load at Sta 117.5 (P^) 
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T-173 Fig. H-3 Mid-Fuselage - Sta 1 16.0 Grid 905 
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Fig. H-4 Mid-Fuselage - Sta 116.0, Grid 905 









5/25/74 MAX-DEF. 


f 



H-6 


T-175 Fig. H-5 Wing Tip - Grid 3017 
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Fig. H-6 Wing Tip - Grid 3017 
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T 178 Fig. H-8 Fin Ballast - Grid 4400 




Fig. H-9 OMS Ballast - Grid 2200 
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T-182 Fig. H-11 Mid Pay load- Stall 7. 5, Grid 4886 
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Fig. H-13 Mid-Fuselage - Sta 116.0, WL 51.5 
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Fig. H-14 Mid-Fuselage - Sta 116.0, WL 51.5 
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Fig. H-15 Nose Torsion - Sta 46.15, Grid 120 
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Fig. H-16 Nose Torsion - Sta 46.75, Grid 120 
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Fig. H-17 Mid Torsion - Wing Tip, Grid 3017 
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Fig. H-18 Mid Torsion - Wing Tip. Grid 3017 
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T-189 Fig. H-1 9 Fin Ballast -Grid 4400 
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T-190 Fig. H-20 Fin Ballast - Grid 4400 
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